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XXII. BAND ERSTES HEFT 1954 


Uber die Lésung ebener Elastizitiitsaufgaben in komplexer 
und hyperkomplexer Darstellung. 


Von K. Stahl. 


1, Einleitung. Bei der Behandlung ebener Elastizitaitsaufgaben sind zwei 
Verfahren bisher weniger beachtet worden, obwohl sie gegeniiber der iiblichen Methode mit 
Hilfe der Airyschen Spannungsfunktion den Vorteil einer wesentlich kiirzeren und iiber- 
sichtlicheren Darstellung bieten. 

G. Kolossoff+leitet die Spannungskomponenten und Verschiebungen aus zwei Funk- 
tionen einer komplexen Veranderlichen ab. Darauf aufbauend hat N. Muschelisvili 2 
einen sehr allgemeinen Lésungsweg fiir die verschiedenen Randwertaufgaben iiber die kon- 
forme Abbildung auf den Einheitskreis angegeben. In den letzten Jahren ist eine ganze Reihe 
bemerkenswerter Aufgaben in komplexer Darstellung gelést worden ®. A. E. Green und 
W. Zerna * verwenden die komplexen Spannungsfunktionen in einer umfassenden Darstellung 
der Elastizitatstheorie. Unabhangig von N. Muschelisvili sind dieselben Funktionen auch von 
A, C. Stevenson® entwickelt und verwendet worden. 

Die zweite Methode, die auf L. Sobrero® zuriickgeht, gestattet, die Spannungen als Ele- 
mente einer Funktioneiner hyperkompl!lexen Verdnderlichen anzugeben. Hiermit 
ist neuerdings von K. Schmidt’ eine Reihe von Aufgaben des ebenen Spannungszustandes be- 
handelt worden. 

Die vorliegende Arbeit § weist die enge Verwandtschaft dieser beiden Methoden nach und 
gibt Transformationsformeln an, die ein Ergebnis in der einen Darstellung leicht in die andere 
iiberfiihren lassen. Anhand der durchgefiihrten Beispiele lassen sich Anwendungsmidglich- 
keiten und Vorteile der beiden Methoden miteinander vergleichen. 


2. Die komplexen und hyperkomplexen Spannungsfunktionen. Die Probleme des ebenen 
Spannungs- und Formanderungszustandes in einem homogenen, elastisch-isotropen Medium 
lassen sich nach Airy auf die Aufgabe zuriickfiihren, eine Funktion U (x, y) zu bestimmen, die 
bei Abwesenheit von Massenkraften in dem vorgegebenen Gebiet der Differentialgleichung 


AAG = 0 (1) 


geniigt und gewisse Randbedingungen erfillt. Die Spannungskomponenten, in tiblicher Weise 
mit O,, dy und t,, bezeichnet, ergeben sich dann zu 


_ eu _ eu Seu 5 
mie Ge? Ht 2 aety" o 

Durch 
U(x, y) = V(z,2) (3) 


1 G. Kolossoff, Z. Math. Phys. 62 (1914), S. 384. 

2 N. I. Muscheligvili, Math. Ann. 107 (1933), S. 282; Z. angew. Math. Mech. 13 (1933), S. 264; 
Einige Grundaufgaben der mathematischen Elastizitatstheorie, 3. Aufl. Moskau 1949. ; 

3° B. L. Mincberg, Prikl. Math. Mekh. 12 (1948), S. 415 (russ.); — D. V. Vainberg, Prikl. Math. 
Mekh. 13 (1949), S. 151 (russ.); — G. N. Polozhii, Prikl. Math, Mekh. 13 (1949), S. 297 (russ.); Dokl. 
Akad. Nauk. SSSR 73 (1950), S. 49 (russ.); — S. E. Birmann, Dokl. Akad. Nauk. SSSR 62 (1948), 
S. 187 (russ.); — C. A. M. Gray, Quart. Journ. Mech. appl. Math. 4 (1951), 5, 444. 

4 4. E.Green u. W. Zerna, Theoretical Elasticity (noch nicht veréffentlicht). 

5 A. C, Stevenson, Philos. Mag. 34 (1943), S. 766; Proc. Roy. Soc. A 184 (1945), 5.129 u, 218. 

6 I. Sobrero, Memorie della Reale Academia d’Italia 6 (1934/35), S. 1; Theorie der ebenen Elastizitat, 
Hamburger Math. Einzelschr., Heft 17, Leipzig 1934. 

7 K. Schmidt, Ing.-Arch. 19 (1951), S. 324. ; 

8 Dissertation Hannover 1953 (D 89); Hauptberichter Prof. Dr.-Ing. 0. Flachsbart, erster Berichter: 
Dozent Dr.-Ing. W. Zerna, zweiter Berichter: Prof. Dr.-Ing. E. Pestel. Die Arbeit entstand am 
Festigkeitslaboratorium der TH Hannover (Leiter Prof. Dr.-Ing. 0. Flachsbart). 
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werde eine Funktion der beiden komplexen Variablen z = x +-iy und 2 = x — iy definiert, 
die nur reelle Werte annehmen kann. Fiir diese Variablen gilt 


ees ete 2 
Ox Oz | Oz’ oy (5, 5) (4) 
ie ead leet)? 
Oz Ox dy’ Oz Ox Oy 
und daher ar 
AU=4,7—. 


Aus (1) folgt, daB AU eine harmonische Funktion ist und somit auch als Realteil einer Funktion 
der Verdnderlichen z etwa wie folgt 1 geschrieben werden kann: 
2 “, 
AU =4 27 = Re [4 O()] =2[9() + O@)I. (5) 
Hier ist ®(z) eine Funktion von z und @ (z) die dazu konjugiert komplexe Funktion. Durch 
Integration von (5) beziiglich z ergibt sich mit der Bezeichnung (z) = [@ (z) dz 


OV Ss Bp 

42 2 [p() +2 OG) +9@)]. (6) 
wobei y(zZ) eine Funktion von z darstellt, die als Integrationsfunktion hinzutritt. Eine Inte- 
gration von (6) beziiglich z liefert mit y.(z) = fe) dz 

AV =2 [2 plz) +29 (2) +2) + Y2(2)] = 4 Re [2 pz) + y2(2)]. 

Hierbei muB y,(z), die zu p,(z) konjugiert komplexe Funktion, als Integrationsfunktion hinzu- 
treten, da V(z,z) voraussetzungsgem48 nur reeller Werte fahig ist. Mit Hilfe von (2), (3), (4) 
und mit der Bezeichnung Y(z) = y’(z) lassen sich nun leicht die folgenden Spannungskombi- 


nationen herleiten: 
OV = 
Ox 1 Oy = 40 — = 2[G(z) + O(@)], (7a) 


uae 


= 2[F O' (2) + P(a)]. (7b) 


Oy — Gy fj 2UTey = 4 
Aus diesen Kombinationen kénnen die Spannungen selbst in einfacher Weise gewonnen 
werden. 
Die Verschiebungen in Richtung der x- und y-Achse seien mit u und v bezeichnet. Im 
Falle des ebenen Formanderungszustandes geben bekanntlich die Gleichungen 2 


7) 
26 = 0, — (6, + dy), 
dv 
A ie Mt Uh) (8) 


die Verkniipfung der Verschiebungen mit den Spannungen an. Dabei ist mit G der Schub- 
slats und mit y die Querzahl bezeichnet. Zur Abkiirzung sei die Verdrehung eines Scheiben- 
elementes 


gesetzt. Die Gleichungen (8) sind gleichwertig mit 


a(u + iv) 1 
2 0 (Oy Oe 2 tel 


Ouray sD = 
26" = 5 (1 — 2») (6, +0,)— 2618, 


ae ; j é 
ie gies wie auch im folgenden werden die von N. I. Muscheli3villi eingefiihrten Bezeichnungen be- 
zt. 


2 Vel. z. B. C. B. Biezeno und R. Grammel, Technische Dynamik Bd. I, S, 108, 2. Aufl. Berlin 1953. 
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wie durch (4) nachzuweisen ist. Durch Vergleich mit (7) folgt 


O(utiv) av 
Po rar ramen ie ak ve) 
O(ut+iv) _ eV SG 
2G — = (1— 29) 25 — 2618. (9b) 
Die Integration von (9a) beziiglich Z ergibt 
26(u +i) =— 2745), (10) 


wobei ¢ (2) eine Funktion von z darstellt, die als Integrationsfunktion hinzutritt. Eine Dif- 
ferentiation von (10) nach 2z liefert 


d(u+ iv) CAV arco 
eens 5,55 69S) 
woraus sich durch Vergleich mit (9b) 
~ 2 Cie 
#2) = (4—49) 2618 (11) 
ergibt. Mit Hilfe von (5) kann nun die Funktion % (z) bestimmt werden: 
P'(2) =(4—4) Bz), G(2) =(4—4 7) 9). (12) 
Durch Einsetzen von (6) und (12) in (10) ergeben sich die Verschiebungen mit x = 3 — 4 zu 
26 (u +iv) =xp(z)—2 0) —H@) (13) 


fiir den Fall des ebenen Formanderungszustandes. Die Ausgangsgleichungen fiir den ebenen 
Spannungszustand haben die Form der Gleichungen (8), es ist lediglich y durch »/(1 +) zu 
ersetzen. Dementsprechend gilt Gleichung (13) fiir den ebenen Spannungszustand, wenn 
% = (3 — y)/(1 +7) gesetzt wird. 

Mit den hier gegebenen Formeln ist das ebene Problem der Elastizitatstheorie auf die 
Aufgabe zuriickgefiihrt, die in dem vorgegebenen Bereich anal ytischen Funktionen @(z) und 
W(z) zu bestimmen, die auf dem Rande gewisse vorgegebene Bedingungen erfiillen. 

So wie die Elemente einer Funktion der komplexen Verdnderlichen z = x +1 y auf Grund 
der Beziehung 1 +i =0 der Potentialgleichung geniigen, so erfiillen die Elemente einer 
Funktion der hyperkomplexen Variablen s = x +jy die Bipotentialgleichung, wenn fiir 7 


(1 +322 =0 (14) 
gefordert wird. Die imaginadre Einheit 7 ist dann durch 

: 1 ; : 

wee (3 + J*) (15) 


mit der hyperkomplexen Konstanten j verkniipft. 
Eine Funktion von s kann in der Form 


f(s) =a4+jb+j?ce+j%d (a,b, c,d reell) (16) 
geschrieben werden, da (14) erlaubt, alle héheren Potenzen von j durch die ersten vier auszu- 
driicken. Es kann nun gezeigt werden}, daB durch 

On =—C, Oy =4, Try Had (17) 


die Spannungskomponenten mit den Elementen einer hyperkomplexen Spannungsfunktion 
zu identifizieren sind. 

Damit ist andererseits das ebene Problem der Elastizitatstheorie auf die Aufgabe zuriick- 
gefiihrt, die in dem vorgegebenen Bereich analytische Funktion f(s) zu bestimmen, die auf 
dem Rande gewisse vorgegebene Bedingungen erfiillt. 


3. Allgemeine Eigenschaften der Spannungsfunktionen. a) Transformations- 
formeln. Nach T. Levi-Civita kann jede hyperkomplexe Zahl und damit auch jeder Wert, 
den die Spannungsfunktion f(s) annimmt in der Form 


f(s) =a +01 +5") 8 (18) 
1 K, Schmidt, Ing.-Arch, 19 (1951), S. 326. 
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durch zwei komplexe Zahlen « =o’ +i’ und Bp =f’ +i f” dargestellt werden. Die An- 
wendung von (15) und (14) fiihrt dann zu 


; , . 3 4/7 9 4/7 2 Zs 4/ o 4s 
fl) ad tp pple ee 5 ppp pe (19) 
in der Form von (16), und mit (17) ergeben sich daraus die einfachen Beziehungen 
o, toy = Re [a], (20a) 
Oy — Oz 2itzy =4 +28 (20b) 


fiir die oben betrachteten Spannungskombinationen. Im folgenden sei vorausgesetzt, daB 
sich die hyperkomplexe Spannungsfunktion in eine Laurentsche Reihe entwickeln 1aBt, also als 


fos 2 arts > ae (21) 


r=— © r= — 


geschrieben werden kann, wobei die hyperkomplexen Konstanten nach T. Levi-Civita durch 
die komplexen Konstanten a, und f, ausgedriickt sind. Diese Voraussetzung ist in den prak- 
tischen Anwendungen so gut wie immer erfiillt. Durch 


s=-[1+0+/4(2 1) (22) 


sind hyperkomplexe und komplexe Variable miteinander verkniipft und die Anwendung des 
binomischen Satzes liefert hieraus 


Beep aRta tei | 


fiir die Potenzen von s. Damit wird (21) zu 


-+ 00 ++ 9° + 9° 
f(s)= >So #+04+7) 63) De RE >], 


Der Vergleich mit (18) zeigt, daB « eine Funktion von z allein ist, und daher kann aus (20a) 
und (7a) geschlossen werden 


+ oo 
0 2 = AO (a), 


r= — oo 
Kine Differentiation dieser Gleichung liefert 
+ © 


2D) Pope a2 AD (Zz), 
Damit ist nun unter Benutzung von (18), (20b) und (7b) analog dem Vorhergehenden zu 
schlieBen : 


fo. 3] 


+ 
> Bw = Pz) — 2 @(z) +2 O'(z). 


Diese Ergebnisse erlauben es, durch 

f(s) 4 OG) = (Le) 2 OU) ee OG) (23a) 
einen Zusammenhang zwischen der hyperkomplexen und den komplexen Spannungsfunktionen 
herzustellen, Werden unter G(s), O’(s) und Y(s) diejenigen Funktionen verstanden, die sich 
aus P(z), P’(z) und ¥(z) ergeben, wenn z durch s ersetzt wird, so kann die Gleichung (21) mit 
den obigen Beziehungen auch in der Form 


f(s) = 4 O(s) + (1 +3°) [P(s) — 2 @(s) + s (s)] (23b) 
geschrieben werden. 
Die Gleichung (23a) zeigt besonders anschaulich die enge V. i 
g ge Verwandtschaft der beiden 
Methoden und Gleichung (23b) gibt die Méglichkeit, ein Ergebnis in der einen Darstellung 
leicht in die andere zu transformieren. Von Nutzcn ist auBerdem die Beziehung 
F(s) = J f(s) ds = 49) +(1 +7) [p() —39(@) +2 ©@)], (24) 


die sich aus (23b) durch Integration und Einsetzen von (22) ergibt. 
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b) Umlaufrelationen, Es sei AB ein Stiick einer Kurve, die im Innern des von 
dem betrachteten Kérper eingenommenen Gebietes verliuft (Abb. 1), und dX, dY seien die 
in x- und y-Richtung von der rechten Seite her auf das Bo- 
genelement dz wirkenden Krafte. Hier ist, wie auch weiterhin, 
unter Kraft der auf die Langeneinheit der Dicke entfallende 
Anteil der Gesamtkraft zu verstehen. Dann verlangt das 
Gleichgewicht 


dX +idY = (6, dy —t,y dx) +i(t, dy — Oy dx). 
Diese Gleichung geht unter Beachtung von (2) und (4) 


iiber in 


Abb. 1. 
dX ~- CGY = (5 Fe | — ap a(s) is Zum Gleichgewicht am Bogenelement. 
oy Ox Oz 
Das positive drehende Moment dieser Krafte bezogen auf den Koordinatennull punkt ist dann 
ae Ae Re |a(23  — a 


Ist C eine ganz im Innern des Bereiches liegende geschlossene Kurve, und sind P, und P, die 
Resultierenden der innerhalb von Cin x- und y-Richtung wirkenden duBeren Krafte, M deren 
Moment bezogen auf den Koordinatennullpunkt, dann ist mit P = P, +i Py unter Ver- 
wendung von (6) 


P = il|pe) +2 ©) +H lc. (25a) 
M = fte|z z ©) +2 9) -— ple (25b) 
fiir einmaliges Umlaufen von C im positiven Sinne. 
Eine entsprechende Forme] fiir die hyperkomplexe Funktion 
IEG) =a By jas. tjto.. Ej* P 
kann aus (24) mit Hilfe von (19) und (25) oder auf direktem Wege erhalten werden. 


In mehrfach zusammenhiangenden Bereichen sind neben den beiden Gleichgewichts- 
bedingungen (25) drei Eindeutigkeitsrelationen! zu beachten. Die Forderung nach Ein- 
deutigkeit der Verschiebungen fiihrt mit (13) zu 


xp (2) — + ®@) — p@)lc =0. (26) 
Fir die hyperkomplexe Funktion ergeben sich aus (24) in der Form von (19) und durch Ver- 
gleich mit (13) die Gleichungen 


[= 14° 0 +A] FO) = Re eo) — 2 BH) —PO] 4° =26u4-- 
ee) Re Sip) — OBS va 2 Cojo 


deren erstes bzw. zweites Element die Verschiebung u bzw. v bis auf eine Integrationskonstante 
angibt. Die Verschiebungen sind eindeutig, wenn bei einmaligem Durchlaufen der Kurve C 
das erste bzw. zweite Element der obigen Gleichungen zu seinem Ausgangswert zuriickkehrt. 


Gleichung (11) in Verbindung mit (12) zeigt, daB die Verdrehung eines Scheibenelementes @ 
bis auf einen konstanten Faktor durch den Imaginarteil der Funktion @(z) bestimmt ist. 
Daher ist @ eindeutig, wenn @(z) eindeutig ist. 


c) Verhalten gegeniitber einer Parallelverschiebung des Ko- 
ordinatensystemes. Es sei ein Koordinatensystem («,, y,) gegeniiber dem (x, y)- 
System um x bzw. yp in x- und y-Richtung verschoben. Die Spannungen scien im (x, y)- 
System gegeben durch @(z), Y/(z) und f(s) und im (x, y;)-System durch )P(z,), ov (2) und 
of (s1), wobei 

By =F — 2%, % = Xp UY; S$] =S— 8p, 8 = % +I Y%o- 


1 Vgl. C. B. Biezeno und R. Grammel, Technische Dynamik Bd. I, S$. 110, 2. Aufl. Berlin 1953. 
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Auf die Spannungskombinationen (7) kann die Verschiebung des Koordinatensystemes keinen 
EinfluB haben, es muB also gelten 
D(z) = oP (a), (27a) 
z D'(2) + P (2) = 21 oP'(%) +P (%)- 


Hieraus folgt mit der Ableitung von (27a) 
(2) = oF (4) — 2 oP (4). (27b) 
Fiir die hyperkomplexe Funktion gilt die einfache Beziehung 
f(s) =of (s1)- (28) 
Wird die Transformationsgleichung (23b) auf das verschobene Koordinatensystem bezogen, 
so ergibt sich mit (28) 

f(s) = 4, ®(s,) +(1 Sei) [oY (s;) — 2 pP (sy) +4 oP (s,)], . (29) 
eine Formel, die dann von Nutzen ist, wenn die gegebene Funktion in der Form einer Reihen- 
entwicklung an der Stelle z = 2 bzw. s = Sy vorliegt. 

d) Die-Spannungen in krummilinig ore 
thogonalen Koordinaten. In der (x, y)-Ebene sei 
ein System von krummlinig orthogonalen Koordinaten 
durch die Gleichung 

z=a(t) mit €=€+i7 
definiert. Die Linien € = konst. und 7 = konst. bilden da- 
bei das neue Koordinatennetz. In einem Punkte 2 seien o;,, 
oy und t,y die Spannungskomponenten im (x, y)-System. 
In diesem Punkte sei das (& 7)-System, gebildet aus den 
Tangenten an die Koordinatenkurven, um den Winkel 6 
gegeniiber dem (x, y)-System gedreht (Abb. 2). Dann gilt 


bekanntlich, wenn o;, o, und 7;, die Spannungskomponenten im (£, 7)-System bezeichnen, 


Abb. 2. 
Krummlinig orthogonale Koordinaten. 


1 

as = (0. + ay) —= (Gy — ox) c0825 + try sin 25, 
J 

Gy => (oz + oy) +> (Gy — ox) cos 20 —Tzysin2 0, 
a ; 

te > (oy — ox) sin 2.6 + Try cos 26, 


oder in anderer Form 
Oz 0; = Oy + Oy ? 
O, — 6; + 21%, =e? *9 (a, —o, +2itzy). 
Aus der Definition der neuen Koordinaten folgt 


rid 7 : 30 
@’(¢) Se 
und mit (7) sind daher die Spannungskomponenten durch 


0: +o, =2[G(z) + O@)], 
0 +2i ty, =< 2 Oe) + PC] (31) 


mit ee co(¢) in den allgemeinen Koordinaten £ und 7 ausgedriickt. 
Die hyperkomplexe Funktion, deren Elemente die Spannungskomponenten im (é, )- 
System angeben, hat mit den Bezeichnungen von (18) die Gestalt 


0, +J°*'— fo; +f te, =a +(1 lg i a =a (32) 


wie sich durch (20) und (31) nachweisen 1aBt. Aus Gleichuug (23a) geht hervor, daB in einer 
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paves umiplexen Funktion von s in der Form von T. Levi-Civita der komplexe Faktor von 
1 +-j?) als Summe einer Funktion von z und einer mit z multiplizierten Funktion von z ge- 


schrieben werden kann. Der Faktor von (1 + j?) in (32) wird sich jedoch wegen q'(¢) im all- 


gemeinen nicht als Summe einer Funktion von € und einer mit ¢ multiplizierten Funktion von C 
schreiben lassen. Daher diirfte (32) nur in speziellen Fallen als Funktion der Variablen E+jn 
darzustellen sein. 


4. ae mathematischer Hilfssatz. Das Verfahren yon N. MuscheliSvili in seiner neuesten 
Form stiitzt sich auf einige weniger bekannte Siatze aus der Funktionentheorie, von denen 
hier nur das fiir die vorliegende Arbeit notwendige angefiihrt werden soll. 


Es sei L eine einfache geschlossene, glatte Kontur, versehen mit einem Umlaufsinn, die die 
komplexe Ebene in zwei Gebiete, St und S-, teilt. Mit S+ sei dabei das links von L liegende 
Gebiet bezeichnet. Auf L sei eine Funktion @(t) gegeben (t bezeichnet die Punkte z auf L), 
die der Hélder-Bedingung 

|® (t.) — O(t,)| < Alt, 4)" mit A>O0 und p> 0 
fix zwei Punkte t, und t, auf L geniigt. Dann ist nach einem Satz von J. Plemelj 2 durch 
_ 1 fod 
a fice 
L 

eine Funktion @(z) bestimmt, die sowohl in S*+ als auch in S~ analytisch ist und im Unend- 
lichen verschwindet. AuBerdem strebt @(z) bei einer Annaherung von z an L von links bzw. 
von rechts bestimmten Grenzwerten zu, die mit ®+(t) bzw. ©-(t) bezeichnet seien und fiir 
die gilt 


@* (t) — D-(t) = Dit). (33) 
Eine Funktion von diesem Verhalten werde abschnittsweise analytisch genannt. Ist die Auf- 
gabe gestellt, die abschnittsweise analytische Funktion @(z) zu bestimmen, deren Randwerte 
auf L die Bedingung (33) erfillen und die fiir groBe |z| von der Form 
k 
1 
P (z) Soe al 


ist °, so gibt in diesem Fall 


k 
aye a! P(t) dt a 
0) =z4, | 2: + Pan 2 (34) 
L n=0 
die eindeutige Lésung dieser Aufgabe, denn die Differenz zweier Lésungen wiirde auf der 
Kontur L keinen Sprung aufweisen, also in der ganzen Ebene analytisch sein und im Unend- 
lichen verschwinden. Das ist aber nur méglich, wenn diese Differenz iiberall identisch ver- 


schwindet. 


5. Anwendungen. Bei den folgenden Anwendungen werden wir im allgemeinen von den 
komplexen Spannungsfunktionen ausgehen und die hyperkomplexe Funktion durch eine 
Transformation nach (23b) erhalten. Die Ergebnisse in der 
komplexen Schreibweise sind zum Teil bereits von A. E. Green x 
und W. Zerna* gegeben worden. 


a) Allgemeine Spannungsfunktion fir 4 
die unendliche Ebene. dn derunendlichausgedehnten 
Ebene greife im Koordinatenursprung die Kraft P = P,- 1% Py Py, 
mit den Komponenten P, und P, und das positiv drehende 
Moment M an (Abb. 3). Die Summe aus Gleichung (25a) und ~S 


der mit i multiplizierten Gleichung (26) ergibt 
Abb, 3. Kraft und Moment in unend- 
P licher Ebene. 


lp@le =— te 5- 


1 N. I. Muschelisvili, Singulare Integralgleichungen. 2. Aufl. Moskau 1949. Ubersetzt von I. R. M. 
Radok u. W. G. Woolnough, Melbourne 1949. 

2 J. Plemelj, Monatsh. Math. Phys. 19 (1908), S. 205. 

3 Q(1/z) = Ordnungsfunktion von 1/z bezeichnet Glieder von der GréBenordnung 1/2. 

4 A. E. Green u. W. Zerna, Theoretical Elasticity (noch nicht veréffentlicht). 
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Damit ist aus (25b) und der mit z multiplizierten Gleichung (26) 
nee, 
p2(z)lc a af =| yee 


fiir eine den Nullpunkt in positivem Sinne umschlieBende Kurve C abzuleiten. Diese Be- 
dingungen werden durch die Funktionen 


RP aa EF iM, 
ee sree ee yo Ua Nora Sey oe ee a 
erfiillt. Durch Differentiation werden hieraus die Spannungsfunktionen 
P. 1 seeeer dyes eik ion 35 
O18) oe roan ig aol ae Ed Corer af % 2. 2 (35) 


gewonnen, die in der ganzen Ebene mit Ausnahme des Punktes z = 0 analytisch sind und im 
Unendlichen verschwinden und daher die Lésung der Aufgabe darstellen. 

Greifen die Krafte P,(r = 1, 2,...n) und Momente M, an den Stellen (%,, yo,) der (x, y)- 
Ebene an und herrscht auBerdem im Unendlichen ein vorgeschriebener Spannungszustand, 
dann folgen mit (27) und mit den Bezeichnungen 2, = % — %o,. 29, = Xor +4 Yor die Spannungs- 
funktionen @,(z), ¥%(z) fiir diesen allgemeinen Fall zu 


n 


— P 1 
r I be 
a 2a(%+1) a | Z 
eat 


Py (2) 
(36) 

. e “ { x Py 1 iM 1 Zor P, 1 | / 

Yo (2) oa er er ae Orsi ie 


Die Konstanten J’ und J” sind durch den Spannungszustand im Unendlichen bestimmt, denn 
wie (30) und (31) zeigen, ist 


Re] => (M +N), I’ => (N,— MN) e-2*%, 


wenn JV, und N, die Hauptspannungen und 6 den Winkel bezeichnen, den N, mit der positiven 
x-Achse einschlieBt. Ein imaginarer Anteil von I" wiirde lediglich eine starre Verdrehung des 
gesamten Systems bedeuten, die keinen EinfluB auf die Spannungen hat. Zum Vergleich sei 
hier die entsprechende hyperkomplexe Funktion angefihrt, die mit (29) und den Bezeich- 
nungen Ss, = S — Sg,, Sor = Xr +J Yor folgt: 


ea acl 
fo(s) SD race Sr +T4P 


Pah alivergoe i asters we i ; | 
add sted) 1 lee Sr 2m Ss? | 2(% +1) 1) L adel 3 


(37) 


In dem Sonderfall einer einzigen Kraft im Koordinatenursprung wird diese Funktion zu 


f(s) =sa5- pps (- 0-4) +P] +01. +1-J8 +4) +P0—#) Ph, 


wenn die Kraft in ihren Komponenten dargestellt ist. Dieses Ergebnis stimmt fiir den Fall 
des ebenen Spannungszustandes mit dem von K. Schmidt 1 angegebenen iiberein. 


b) Lésung der Randwertaufgaben fiir den Kreishbereich. Der 
elastische Kérper, eine unendliche Scheibe mit Kreisloch oder eine Kreisscheibe, nehme das 
Gebiet St ein. Die Kontur L habe den Radius 1; der Umlaufsinn sei so festgelegt, daB St 
links von L zu liegen kommt. Durch 


eae ere?, 6 =Inr ig 


werden diesem Bereiche angepaBte Koordinaten eingefiithrt und damit nach (30) die Spannungs- 
komponenten in bezug auf diese Koordinaten zu 


0, +6) = 2[@(z) + b(z)], 
1 K. Schmidt, Ing.-Arch. 19 (ORIN, sh Beil. 
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aot its = 2=[z @'(z) + Y(z)] 
erhalten. Hieraus folgt 


Or tits = B(x) + (2) —2 Oz) 2 VQ). (38) 

Die Ableitung der Verschiebungskomponenten (13) nach dem Winkel 3 ergibt 
O(u +iv ‘ oe Sh ea aN 
26 ts [x OG)— BG) +7 HG) +E VQ). (39) 


Die Veranderliche z auf der Kontur L sei mit t bezeichnet. Dann gilt hier, da LZ durch den 
Einheitskreis gebildet wird, # = 1/t und damit kénnen die Spannurgen und die Ableitung der 
Verschiebungen auf dem Rande L durch 


6, tity |, = (1) 4 Ae) tes ate) (40a) 


t t t 


: sala US aaa I Seay 
nite OO o(-) =o ie cS ¥ (>| (40b) 
ausgedriickt werden. 


Im folgenden sei vorausgesetzt, daf auBere Krafte nur am Rande der Scheibe angreifen. 
Dann sind die Funktionen @(z) und Y(z) in St analytisch. Unter @(1/z), @'(1/z) und W(1/z) 
seien die Funktionen verstanden, die sich aus @(z), @’(z) und Y(z) ergeben, wenn Z durch 1/z 
ersetzt wird. Die Bedeutung von @(z) in dem Gebiet S~ sei durch die Definition 


@(2) = o(7) + : aes ! Ps) far 2 in St (41) 


0(u+ iv) 


2G 55 


2 
& z 


festgelegt. Dann ist ®(z) unter den gegebenen Voraussetzungen eine abschnittsweise analy- 
tische Funktion (mit eventueller Ausnahme des Punktes 2 = 0) und der zweite Teil der Glei- 
chungen (40a) und (40b) ist der Grenzwert von (41) fiir z > t. Mit den Bezeichnungen von (33) 
ist daher 


0, titre |x = Ot (t) — D- (i), (42a) 
22 tty) itt ot + o-()], (42b) 
L 


und es ist ersichtlich, daB @(z) sich durch die Teile des Randes, an denen keine 4ufBeren Span- 
nungen (¢, und T,») wirken, analytisch fortsetzt. 
Liegt zin St, so liegt 1/z in S~. Dieser Wert in (41) eingefiihrt ergibt 


Me) = — Oz) +7 0 @) + 2G — fir 2 in St, (43) 
und der Ubergang zu den konjugiert komplexen Werten liefert hieraus 
1 1 Saal Ike ey 
(2) = + &(2) + Ola) + 6'(), (44) 


womit die Funktion Y(z) durch @(z) ausgedriickt ist. Daher geniigt die Kenntnis von @(z) 
in S+ und S~ zur Bestimmung der Spannungen und Verschiebungen. 

Handelt es sich um die unendliche Scheibe mit Kreisloch, so werden 
die Spannungsfunktionen fiir groBe |z| mit (35) iibereinstimmen, also von der Gestalt 


aan 2 Li le Ps. 1 45 
TOD Sa Pape he: 0%) Ao: =aaen +3) 2 


sein, wenn P die Resultierende der Randbelastungen darstellt. Dann ist nach (41) fiir kleine 


|2| @(z) von der Form 


op Lee J 
~ Qa(e +1) @ pe Eee 


(2) (46) 


Im Falle der Kreisscheibe sind die Spannungsfunktionen fiir |z|< 1 analytisch, 
haben also die Gestalt 


D(z) =A, +4, 2+4,2..., V(z) = B+ By2+ B,2..., (47) 
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woraus sich ergibt, daB fiir grofe |z| die Funktion ®(z) von der Form 
ae tl 
Ce (3) (48) 


sein mu. 
Liegt die erste Randwertaufgabe vor, so ist auf dem Rande L 
6, titel, = N(t) +i T(t) 
vorgegeben. Dabei ist N(t) die nach auBen gerichtete Normal- und T(t) die Tangentialbe- 
lastung in Richtung des Umlaufsinnes von L, bezogen auf die Langeneinheit. Es sei voraus- 
gesetzt, daB N(t) und T(t) der Hélder-Bedingung geniigen. Dann ist nach (34) unter Beriick- 


sichtigung von (46) 
i N LT 1p 1 
®() =57; [ (t) + 4 () a x — (49) 
L 


i 2% (% + 1) 


die eindeutige Lésung im Falle der unendlichen Scheibe mit Kreisloch und unter Beriick- 
sichtigung von (48) 


oie t—z 


@ (2) = ie ceeege (50) 
L 


die eindeutige Lésung fir die Kreisscheibe. Ein Vergleich von (48) mit (47) liefert zusatzlich 
die Bedingungen 


Fe NO ef | 
eo (51) 
0 =lim{O@) +4}1= 7, (VO +iTOa | 
oder z J 
oe 22 
fa | NO ad = Ke [4], [toa =0, 
: (52) 


TIN® +iT@]e? do =0, 


die den Realteil der Konstante A) bestimmen und besagen, da die resultierende Kraft und 
das resultierende Moment der Randbelastung verschwinden miissen. Der Imaginirteil von A, 
ist als starre Drehung hier bedeutungslos. 


Bei der zweiten Randwertaufgabe sind die Verschiebungen am Rande 
u--iv |, = g(t) +i g0(t) 


vorgegeben. Ihre Ableitungen nach ¢ seien mit g; und g, bezeichnet. Fir @(z) besteht dann 
nach (42b) die Randbedingung 


“Or (t) + Ot) =26[g (*) tig, I]. 
Wird eine Hilfsfunktion K(z) durch 
fiir z in St 
Ki) 53 
i ae ey) 
definiert, so mu das Produkt K(z) D(z) die Bedingungen 
KY (2) O* (2) — K(z) O-(z) = 2 Cl gi () +85 (0) 
erfiillen. Unter der Voraussetzung, daB gi (t) und g(t) die Hélder-Bedingung erfillen, ist 


_ 1 2G lanl) + igs (0)] xP 

K @ ee 1 §o 

eee aif ie? “Fae FY) Of 
L 


fiir die unendliche Scheibe mit Kreisloch und 


] us ig’ = 
K (2) (2) =; | 2 6 [gi(t) + igs (t)] oh (55) 


t—2z 
L 
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fir die Kreisscheibe die eindeutige Lésung. Die zu (51) analogen Integralbedingungen besagen 
hier, daB die auf dem Rande vorgegebenenVerschiebungen keinen Sprung aufweisen diirfen 
und geben mit g, + ig) = (g, +ig,) e—** die Bedingung 

2a 

ae I . 
% Ap — Ay = 5h [26( + i gy) dd 

0 
zur Bestimmung von A); dabei sind g, und gy die Randverschiebungen in radialer und tan- 
gentialer Richtung. 

Die gemischte Randwertaufgabe liegt vor, wenn auf einem Teil des Randes 
L’ die Spannungen, auf dem iibrigen Rand L’’ die Verschiebungen vorgeschrieben sind. Hier 
sei der Kiirze halber vorausgesetzt. daB L’ einen Bogen zwischen den Punkten t, und t, 
einnimmt. Die Erweiterung auf den Fall, da® L’ aus n Abschnitten besteht, ist ohne weiteres 
miglich 1, doch von geringerem praktischen Interesse. 

Nach dem Vorausgehenden muf dann @(z) auf L’ und L’” die Randbedingungen 


@* (t) — O-(t)|,, = N(t) +i T(t), 
x O*(t) + O-(t)|z. = 2 [gi (t) +ig(s)] 
erfillen, die kiirzer geschrieben lauten 
@* (t) — k O-(t) = h(t) 
1 auf L’ ( N(t) +i T(t) auf L’ 
mit k | ; und h(t) = } (56) 
—— auf L | — [gi () +g; ()] auf L”. 
Wird (56) mit den Randwerten H(t) einer Hilfsfunktion H(z) multipliziert, die in St und S~ 
analytisch ist mit der Eigenschaft, daB 
k H*(t)— H-(t) =0 
ist, so wird 
H*(t) OF (t) — H-(t) O-(t) = h(t) H(t) 
die Randbedingung, die die Funktion H(z) @(z) erfiillen mu8B. Die geforderte Eigenschaft be- 
sitzt die Funktion 
ee i ae 
H(z) = (z—t,)2 ‘ el tee. ” mit y =e ; 
wenn als Funktionswerte die Werte des Zweiges genommen werden, der in der langs L’’ auf- 
geschnittenen Ebene analytisch ist mit der Higenschaft, da8 lim H(z) =z wird. Unter der 


Z—> 0 


Voraussetzung, da H*+(t) h(t) die Hélder-Bedingung erfiilt, ist die Lésung durch 
1 | h(t)H* (t) ap 12 _ P- H(0)-% (57) 


H(z) O(2) =>] =e (x-+1)2% ' 2a(x+1)z 


L 


fiir die unendliche Ebene mit Kreisloch und durch 
;: 1 h(t) H+) = 
Hs) ®@) =75; | ENN hy (58) 

| is 
fir die Kreisscheibe eindeutig gegeben. Die Konstante A, ist durch die Bedingung ®(0) = A, 
bestimmt. 

Greifen Einzelkrafte und Momente im Innern der Scheibe an und sind die Rand- 
spannungen bzw. Verschiebungen gleich Null, so kann die Spannungsverteilung durch die 
Funktionen 


D(z) = Do(z) +O, (2), PY (2) = M2) + Fi (2) (59) 


beschrieben werden. Dabei sind ®,(z), Y%(z) die Spannungsfunktionen (36), die sich ergeben, 
wenn die vorgegebenen Krafte und Momente in der unendlichen Ebene angreifen wiirden; 


1 N. I. Muschelisvili, Einige Grundaufgaben der mathematischen Elastizitatstheorie S, 463. 3. Aufl. 
Moskau 1949. 
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M,(2) und Y(z) sind in S* analytische Funktionen, die so zu bestimmen sind, daf ihre Rand- 
spannungen bzw. Verschiebungen entgegengesetzt gleich denjenigen von ®,(z), Y%(z) werden. 
Es ist also fiir den Anteil @,(z), Y%4(z) eine Aufgabe analog dem .Vorhergehenden zu lésen, was 
mit den angegebenen Formeln in einfacher Weise méglich ist. 

Ist das Verschwinden der Spannungen auf dem Rande gefordert, so mu nach (40a) 
und (42a) @,(z), deren Bedeutung in S~ durch (41) definiert ist, die Randwerte 


ts = e509 =— 000-85) 82) HERE 


annehmen. Es ist demnach 
1 


rg fal he eal ie 
1 — #5) — (7) + t (7) re} (7) 
- dt (60) 


Ding U t—z 


®, (2) = 


die cindeutige Lésung fiir die unendliche Scheibe mit Kreisloch, da die Resultierende der Rand- 
spannungen verschwindet. Die Lésung erfolgt nach dem Residuensatz und liefert 


®,(2) = b,(-) tee oi(~) Ee ¥, (7) 4 @,(0) far x in St, ” 
P, (2) = Dy(z) + ©, (0) | fiir z in S- 
Die Funktion ¥4,(z) bestimmt sich dann nach (44) zu 
A! 1 1 , 1 = 
Ys(@) = ®5(Z) +3 Pa(F) — > MO) + 100) + BO). (62) 


Die Lésung fiir die Kreisscheibe unterscheidet sich hiervon nur in der additiven Konstante 
und zwar ist ®)(0) durch — A, zu ersetzen, wobei der Realteil von Ay aus 
A, =lim @,(z) =lim 2 Y(z) — Ay 
z+0 %—> oo 
erhalten wird. 


Ist das Verschwinden der Verschiebung auf dem Rande gefordert, so muB die 
Funktion @,(z) nach (40b) und (42b) die Randwerte 


2 P:() + B() = — x B,(0) + ®(-) (7) 7 ¥o(+) 


annehmen. Die Lisung erfolgt analog dem Vorhergehenden und ergibt 


9 =2[0()—40,(2)—-2%(2)— a0), 


& & 


Fs(2) = — 3 Go() +32 ol) +25 1) — > Bie) 


(63) 


fiir die unendliche Scheibe mit Kreisloch, wenn die resultierende Kraft, die auf dem Rande L 
itbertragen wird, verschwindet. Die Lésung fiir die Kreisscheibe unterscheidet sich hiervon 


nur darin, da — A, an die Stelle von @,(0) tritt, wobei A, aus 


Ay = ®,(0) =lim 


Z-* CO 


(1 2B (| 440 
\# Gs Po(*), pres 
erhalten wird. 

Ist auf einem Teil des Randes L’ das Verschwinden der Spannungen und auf dem 
itbrigen Rand L” das Verschwinden der Verschiebungen gefordert, so kann die 
Lésung in der Form 


H(z) ®,(z) = 1 re H*(t) dt 


| Qt t—z 
L 


iF Py (t) ,(-) : ®5(-) | : ¥, (7) creme LY 
UW liteaeel Ie eee al a 
fe Po (t) + | deel ; Bole leaes ¥(7)| aateds 


—§ 


mit h(t) = 
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Mg die unendliche Scheibe mit Kreisloch angegeben werden, wenn die Resultierende der auf L’’ 
ttbertragenen Krafte verschwindet. Die Lésung fiir die Kreisscheibe ergibt sich hieraus durch 
Addition von — A, z. Die Konstante Ay ist durch ©, (0) = Ay bestimmt. 
A pues eine Kombination der hier gegebenen Lésungen ist es immer méglich, eine Aufgabe 
fiir den Kreisbereich, bei der auf dem Rande Spannungen oder Verschiebungen und auBerdem 
im Innern angreifende Krafte und Momente vorgegeben sind, zu behandeln. 

Als Beispiel seien die Spannungsfunktionen fiir eine Kreisscheibe mit freien Rindern an- 


gegeben, in der an den Stellen (ages Yen iS dy en. 8 n) die Krafte P, und Momente M, an- 
greifen: . 


= 1 : tan P, #% Pr P, 
D (z) moat | or ; ae = 2 (zor Zor) 


= or 25 S) 
os 1 . My(% + 1) 
Pigg x ! ge eae 
Hebe — 24) <5 EY Py, + EEN, 
Tp) oe pee Zor Pr  iMr(w+1) , PF 
2 3(#% +1) \ Zr 22 m2 ee 
rT Ls iS 
P, , 2 Pp, , © P, Tr 
re (For Zgr) Zo iS (20 r Zor) : 7a 
r PO Tey “r 
_ iMr(e +1) | 2¢ Mee +1) 
ee rer: 


Darin ist zur Abkiirzung z = z— 2,, 2, = 1/zor = X%0r +1Yor gesetzt. Es zeigt sich, daB 
die Randbedingungen dadurch erfiillt werden, daB zu jeder Singularitaét an einer Stelle 
(%or, Yor) eine zweite an der Stelle (xo,, yo,) hinzutritt. Dabei ergibt sich die Lage der zweiten 
Stelle durch ,,Spiegelung‘ der ersten am unbelastéten Kreisrand. Hierin kann eine Analogie 
zum Schwarzschen Spiegelungsprinzip gesehen werden. 

Ein Ergebnis wie das obige kann nun leicht mit Hilfe von (29) in die hyperkomplexe 
Schreibweise transformiert werden. In diesem Falle wird mit den Bezeichnungen 
Sp SS See Sie = Xo r +7 Yor 


nm 


4, = EE x Py P, , = / 1 <= 
Oliscorrrmery oh) Ie) aT 


m(% + 1) <col Sr se pees 


Pe Mnles Os, Cag) Se _#Pr  iMr(e +1) | 3%Pr , Pr 
an 


2s" { T 2a(e+ sou easy, 32 ou a 
4 P, , P, , x P, —/ _ Dy P , 
= =—(z z % Zz | = (G2 Zz = = ( Zy,)2 
ee ( Or or) $73 ( Or or) l 32 ( Or or) ansor oe ( Or or) 
ee, = 4i M,(z# 4-1) i M,(x% + 1) 27% Mr(# + 1) 
+ 2 Pr (zo, — %r) + (% +1) P, 20, S2 Ue 52 Bode : 
SOT een r eo 


In dem Sonderfall, daB zwei diametral gerichtete Krafte auf einem Durchmesser der Scheibe 
angreifen, stimmt diese Formel mit dem von K. Schmidt + erhaltenen Ergebnis iiberein. 


c) Lésungen fir die Halbebene. Die gerade Linie, die die Begrenzung der 
Halbebene bildet, kann als ein Kreis mit dem Radius Uncudlich aufgefaBt werden. In den oben 
angegebenen Lésungen fiir den Kreisbereich ist daher jeweils die entsprechende Lésung fiir 
die Halbebene enthalten. Insbesondere werden die Gleichungen fiir den Fall, daB der elastische 
Kérper die obere Halbebene einnimmt, aus den entsprechenden Gleichungen fiir die unend- 
liche Ebene mit Kreisloch dadurch gewonnen, da das Koordinatensystem in den oberen 
Randpunkt des Kreises verschoben und die Langeneinheit so gewahlt wird, daB der Radius 
des Kreises gegen Unendlich strebt. Es werde dabei vorausgesetzt, daB alle duBeren Krafte 
in der Halbebene in einer endlichen Umgebung des Koordinatenursprunges angreifen. Dann 
wird durch diesen Grenziibergang aus (41) die Definition 


@(z) = — @ (z) — z @’ (z) — Wz) fir zin S- (65) 


1 K. Schmidt, Ing.-Arch. 19 (1951), S. 336. 
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fir die Funktion @(z) in der unteren Halbebene S~ gewonnen. Hier sind unter D(z), D(z) und 


Wz) die Funktionen zu verstehen, die sich aus @(z), D(z) und W(z) ergeben, wenn z durch z 
ersetzt wird. Liegt zin der oberen Halbebene St, so liegt z in S-, daher gilt ebenfalls 


@(z) =— O(z2)—z O'(z)— Ve) fiir z in St, ' 
woraus durch Ubergang zu den konjugiert komplexen Werten folgt 
W (2) = — D(z) —z O'(z)— G(z) _— fiir z in St. (66) 


Diese Gleichung erlaubt, die Funktion Y/(z) durch @(z) auszudriicken. Es geniigt also auch hier 
die Kenntnis von ®(z) in St und S~ zur Bestimmung der Spannungen und Verschiebungen. 


Aus den Gleichungen (49), (54) und (57) werden mit h(t) nach (56) die Lésungen 


(2) =5n5 [a 4h, } 
ii 
Tee. cis igs 
L 
NIG) es) eee oe 
L 


der ersten, zweiten und der gemischten Randwertaufgabe fir die 
Halbebene erhalten. Dabei ist L die von — o bis + co zu durchlaufende reelle Achse 
wahrend die ibrigen Bezeichnungen beibehalten sind. 


Greifen Einzelkrafte und Momente im Innern der Halbebene an und ist auf dem Rande L 
das Verschwinden der Spannungen bzw. der Verschiebungen vorgeschrieben, so werden aus 
(61), (62) und (63) die Lésungen 

P, (2) = — Bo(z) — 2 Do(z) — Molz), Pilz) = Po(z) — = Pilz) + z Bo(z) 


fir die erste Randwertaufgabe, beziehungsweise | 
eee Su) ‘ia = , 
PD, (2) = —- [Po (2) +2 (2) +¥(2], Y(z) =x Dy(z) — ®,(z) — z Bi (z) 


fiir die zweite Randwrrtaufgabe gewonnen. Ist auf einem Teil des Randes L’ 
das Verschwinden der Spannungen und auf dem Rest des Randes L’’ das Verschwinden der 
Verschiebungen gefordert, so wird die Lésung aus (64) zu 


H(z) ®,(2) = 1 ee H+ (t) ap 


200% t—z 


1D, (1) Dy (1) == t. On (1) =) eae 


mit A(t) =) _ ®,(t) 4+— [B,(t) tO; (t) + Y(t)] auf L’”’ 


erhalten, wenn die resultierende Kraft, die auf L’’ itbertragen wird, verschwindet. 

Als Beispiel seien die Spannungsfunktionen angegeben, die sich fiir eine in der oberen 
Halbebene erstreckende Scheibe ergeben, deren Rand unverschieblich festgehalten wird und 
in der an den Stellen (%,, yor) (r = 1, 2,...n) die Krafte P, und Momente M, angreifen: 


1 < P, 
Oy 


— 1 < Pr Zor Pr i Mr(x + 1) x Py Pr Zor Pe , 
BO) saga) 2 ar i al ea earatale : —- 72 (20+ — Zor) 


z 2 12, 
ip a ze a, z u z 


P, 


12 
% By 


| 
| 


(207 — Zor) 4 


P, i M,(x + 1) 
a % 2 x 


72, I 73 
% By BY 


iMy(%#+1) , mace 


. Uf Uh —_ ° . ° . . 
Hier bedeutet z = z— %, und 2, =%, = x9, +i yor. Durch (29) wird hieraus leicht die 
entsprechende hyperkomplexe Spannungsfunktion mit den Bezeichnungen s, = s — so,, 
Sor = Xr +J Yor gewonnen: 
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peaks Soil whet ee i My (x 
OSes + bg (Hoe — aoe) $e 


_ Sr 37 eg 8 % s°3 J 
; n — Bas - 
+ 1S on joe — EM TD iP teh See ees c 
2a (x + 1) oy Pan s2 Pd ae (20, — Zor) 
a 3iM.(z+1) , 2iM 
eee, po ae r r(% +1) 7 
=i 4 oe (Zor Zr) oe 372 -+ Be 373 (zor — ‘on ° 


Bei gradlinigen Begrenzungen gestaltet sich die Lésung der sogenannten dritten 
Ran d wertaufga b e, bei der auf dem Rande die Tangentialspannung und die Normal- 
peers vorgeschrieben sind, besonders einfach. Aus (7) und (13) leiten sich die Glei- 
chungen 


Ep tates O(@) +20'(@) + VQ), 
2G F**) _ 4 (2) — 6@)— + OG) — ¥@) 


ab. Auf dem Rande L gilt hier z =z =t und daher kénnen die Spannungen und die Ab- 
leitung der Verschiebungen auf dem Rande auch durch 


Oy —itey|z = N(t) +iT() = Bt) + D(t) +t O'(t) + V(t), 

2 G20) = 26 (8h) +igi)] =x PH) — GB —1 Hy — FH 
angegeben werden. Der Imaginirteil dieser Gleichungen liefert die Beziehungen 

2% T(t) =t O'(t)—1t O'(t) + P(t) — P(t), 
2i[ T(t) +26 g(t] = (x +1) [O() — OO). 
Dieses Ergebnis legt es nahe, die Fortsetzung der Funktionen @(z) und P(2) in S~ durch 
D(z) = O(2) fir zin S-, W(z) = (2) fir z in S- (68) 

zu definieren. Dann setzen sich D(z) und Y(z) durch die Teile des Randes, auf denen Tan- 


gentialspannungen und Normalverschiebungen verschwinden, analytisch fort. Fiir O(z) und 
W(z) ergeben sich dann die Randwertaufgaben 
Ont ’ 

@* (1) — OW) =TGITH +26 HOI, 


Wt (t) — W(t) = — 2i T(t) + t[O'(t) — O'(0)]. 


Wird vorausgesetzt, daB T(t) und go(t) sowie deren erste Ableitungen der Hélder-Bedingung 
geniigen, so ist die Lésung durch 


1 ( 2i(T() + 26830) 
2) =s5| OPN 2 at, 
L 
P (2) = (eee oe a 
~~ Dae t—z 
L 


eindeutig gegeben. 

Greifen Einzelkrifte und Momente im Innern der Scheibe an und ist auf dem Rande das 
Verschwinden der Tangentialspannung und der Normalverschiebung gefordert, so kann durch 
eine Rechnung analog der Vorausgegangenen die Lésung 

D(z) = Dp(2) + Pole), PY @) = Pole) + Yol@) (69) 
gefunden werden. 

Soll dieses Ergebnis in hyperkomplexer Form dargestellt werden, so ist es zweckmaBig, 
durch f,(s) eine Funktion zu definieren, die aus f(s) dadurch entsteht, daB iiberall j durch — j 
ersetzt wird. Wie (15) zeigt, geht dadurch gleichzeitig 1 in — i tiber. Auerdem sei fo(s) die- 
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jenige Funktion, die aus fo(S) entsteht, wenn s durch s ersetzt wird. Mit dieser Bezeichnungs- 
weise ist die Lésung der obigen Aufgabe durch 


f(s) =fols) +.o(s) (70) 
gegeben, wie durch (23b) nachzupriifen ist. 


d) Lésungen fiir-Bereiche, die sich durch rata vale Funk- 
tionen auf den Kreis abbilden lassen. Die elastische Scheibe sei von einer 
einfach geschlossenen, glatten Kontur L begrenzt. Durch die rationale Funktion 


z = w(C) (71) 


werde der von der Scheibe eingenommene Bereich auf das Gebiet S+ der ¢-Ebene konform ab- 
gebildet, das von dem Einheitskreis A begrenzt ist. Dabei sei die Abbildung so vorgenommen, 
daGB im Falle cines endlichen Bereiches, dieser auf das Innere, ein unendlicher Bereich auf das 
AuBere von A abgebildet wird. Der Umlaufsinn von A wird so festgesetzt, daB das Gebiet S* 
links von A zu liegen kommt. Die Funktion w(¢) hat dann die Form 


w(6) =cl +a,(6), (72) 


wobei «,() eine in St analytische Funktion und c eine komplexe Konstante darstellt. Unter 
@() und ¥(€) seien diejenigen Funktionen von ¢€ verstanden, die sich aus @(z) und W(z) er- 
geben, wenn z durch w(C) ersetzt wird. Dann ist 


@'(C) = w'(C) O'(z). 
Die Funktion w(¢) mit € = ge’ definiert gleichzeitig in der z-Ebene ein krummliniges System 
der Koordinaten @ und © mit der Eigenschaft, daB die Kurve 9 =1 mit der Kontur L zu- 
sammenfallt. Die Spannungskomponenten in Bezug auf diese Koordinaten sind nach (31) 
Og +09 =2[(9(C) + @(C)], 
ow! (C)o Oi 


: Cilla ¢) 
Vo 2 oo = es 
59 — oe +2itee = 2 =a) iy + YO), 


woraus 


0 +i mo = O) + OO) — Sp (OO OO +0'@) YO) 
zu erhalten ist. Die Ableitung der Verschiebungen (13) nach O ergibt 


2 Uti — ie wt) (x OO)— O@ +2 oO #@ +H O FON. 


Ks erweist sich hier als zweckmafig, die Funktion w’(¢) @(¢) als Hinheit zu betrachten. Die 
Fortsetzung dieser Funktion in S~ sei durch 


wl) ®G) =— ai) ®(=) +z[o0) Be) Ha! (=) ¥(z)| ftir fim S~ (78) 


definiert. Dabei sind unter @ (1/6), ®'(1/¢), W(1/¢) wud @(1/C) diejenigen Funktionen zu 


verstehen, die sich aus ®(), ©’ (¢), Y() und w(¢) ergeben, wenn € durch 1/€ ersetzt wird. 
Diese Definition ist so gewahlt, daB sich die Funktion @(¢) durch die Teile der Kontur A, die 
Bilder der unbelasteten Stiicke des Randes L sind, analytisch fortsetzt. Liegt ¢€ in St, so 


liegt Wg in S~, Dieser Wert als Veranderliche in (73) eingesetzt ergibt nach Ubergang zu den 
konjugiert komplexen Werten 


a ()¥O) =s6' (=)[e© | ®(=)| a(z)®'() fiir Cin St, (74) 


Diese Gleichung erlaubt, die Funktion Y(C) durch @(¢) auszudriicken. Es geniigt daher die 
Kenntnis der Funktion @(¢) in S+ und S~ zur Bestimmung der Spannungen und Verschie- 


bungen. Die Werte von € auf dem Kreise A seien mit 4 bezeichnet. Dann ist 2 — 1/a und 
damit kénnen die Werte der Spannungen und Verschiebungen am Rande durch die Rand- 
werte der Funktion w’‘(¢) O(C) ausgedriickt werden: 


[ar'(2) P(A) — [w'(A) OA) = (A) [N(@) +i TA), 
#4 [w'(A) OAT + [w'(a) O(a) = 2 E[gi(a) tigs(ay]. | i‘ 
Hier bedeuten g; und g die Ableitung der Verschiebungen g, und gy nach J. 
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Greifen dufere Krafte nur am Rande des elastischen Kérpers an, so sind @(C) und Y(¢) in 


S* analytisch und daher wegen (72) auch w'(£) @(6) und w’'(£) W(C). Im Falle eines unendlichen 
Bereichs ist nach (45) fiir grofe |¢| 


: - 1 1 
00) 0) =—saaay FI (=), 
; = Pz 1 1 
HPO =a ey et tole): 
wobei P die Resultierende der auf L angreifenden auBeren Krifte darstellt. Die rationale 
Funktion w’(¢) habe an den n Stellen €, (r = 1,2,...n), die in S~ liegen, Pole von der maxi- 
malen Ordnung m,. Wird w‘(¢) O(€) in das Gebiet S~ fortgesetzt, so zeigt (73), daB w’(f) B(C) 
ebenfalls Pole an den Stellen ¢, von der maximalen Ordnung m, aufweist. Hinzu kommt im 
Falle eines unendlichen Bereiches an der Stelle ¢ = 0 derPol 


00) © =a = +90), 


der sich durch Anwendung von (73) auf (76) ergibt. Damit ist die Lage aller Pole der Funk- 
tion w’(¢) O(¢) und deren maximale Ordnung bekannt. Eine beliebige Funktion @(f) wiirde 
jedoch, wie aus (74) hervorgeht, bei der Funktion w’(¢) Y(¢) an den Stellen 6) = Wie die in St 
liegen, Pole von der maximalen Ordnung m, hervorrufen, was den Voraussetzungen wider- 
spricht. Die Bedingung, daB w’‘(¢) Y(C) in St+ analytisch bleibt, liefert fiir jede Stelle Cf ein 
System von m, linearen Gleichungen, die die m, ersten Koeffizienten der Potenzreihen- 
entwicklung der Funktion @(€) an der Stelle €; mit den Koeffizienten der Hauptteile des 
Poles der Funktion w’(¢) O(¢) an der Stelle ¢, in Beziehung setzt. Diese n Gleichungssysteme 
seien kurz mit 


(76) 


ea) eet Ot a) (77) 
angedeutet. Besonders einfach werden diese Bedingungen wenn w(¢) die Form einer Laurent- 
schen Reihenentwicklung an der Stelle € = 0 hat, da dann Pole nur an den Stellen €= 0 und 


¢ = oauftreten. 


Es sind nun nach (34) bei den iiblichen Voraussetzungen iiber die Stetigkeit der vorge- 
gebenen Randwerte 


' 1 N(A) + iT(A)] w’(A P 1 
ee ge ERO) eeeceay a 
A 
; fa 2G) +120 P 1 
K() a0) BC) ae, | EP ah + RO aa ey 
A 
; 1 of h(d) H+(A Px H(0) 1 P 78 
HG) 0'@) OC) =a; [AGE + RO+ aaah eae ED |S 
A 


a! @)EN(A) +8 TQ) aut A’ 
mie) 28 igi (a) tig (auf A” 


x 
die Liésungen der ersten, zweiten und der gemischten Randwertaufgabe fiir einen unendlichen 
Bereich. Dabei sind K(¢) und H(¢) die zu K(z) und H(z) analogen Hilfsfunktionen in der C- 
Ebene, A’ und A” die Bilder von L’ und L” und R(€) ist eine rationale Funktion, die nur an 
den Stellen ¢, Pole der maximalen Ordnung m, besitzt und im Unendlichen verschwindet. Die 


Lésung (78) liefert mit (77) ein System von = m, linearen Gleichungen fiir die Koeffizienten 
der Hauptteile der Pole von R(¢). Nach Lisung dieses Gleichungssystemes ist R(C) voll- 
standig bestimmt. 

Im Falle eines endlichen Bereiches ist P — 0 zu setzen; dafiir treten bestimmte additive 
Konstanten und bei der gemischten Randwertaufgabe ein Pol ersten Grades an der Stelle € =oo0 
hinzu. Diese sind durch (73) mit den Koeffizienten der Potenzreihenentwicklung von ®(¢) 
an der Stelle £ = 0 verkniipft und durch Untersuchung des Verhaltens der Lésung (78) im 


Unendlichen zu bestimmen. 
2 
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Es wurde bereits gezeigt, da® sich die Aufgaben, bei denen im Innern des Kérpers an- 
ereifende Krafte vorgegeben sind, auf den hier behandelten Fall, da® Krafte nur am Rande 
angreifen, zuriickfihren 1a8t. Es kann daher auf eine besondere Behandlung an dieser Stelle 
verzichtet werden. 

e) Angeniherte Liésung fiir Bereiche, die von einer belie bi- 
een einfach geschlossenen Kontur begrenzt sind. Zu jedem Bereich 
S, der von einer beliebigen einfach geschlossenen Kontur begrenzt ist, existiert eine Funktion 


z= o(C) =cf + (C), 
die diesen Bereich auf das AuBere beziehungsweise auf das Innere des Einheitskreises konform 
abbildet. Die Funktion «,(¢) ist in dem Bildbereich von S analytisch und kann daher durch 
die Entwicklung 


(Oy Oe ee 
fiir einen unendlichen Bereich und durch 
Cy (5) = €2 C? + 050% ++ ° 


fiir einen endlichen Bereich dargestellt werden. Wird diese Entwicklung nach einer be- 
stimmten Anzahl n von Gliedern abgebrochen, so wird durch die Funktionen 


Onl GC) ==0'C Pe Ss eo a 


Wn(C) = C6 + ¢, ¢? ¢, 8 es Cn Ge 

ein Bereich S, auf das AuBere bzw. Innere des Einheitkreises abgebildet, der sich von S um so 
weniger unterscheidet, je gréBer n gewahlt wird. Die Funktion w,(C) ist rational, daher ist die 
Lésung der entsprechenden Aufgabe fiir den Bereich S, mit den oben angegebenen Formeln 
méglich. Diese Lésung wird sich von der gesuchten um so weniger unterscheiden, je weniger 
sich S, von S unterscheidet. Durch Wahl eines geniigend groBen n ist daher die gesuchte 
Lésung mit beliebiger Genauigkeit zu erreichen. Als Beispiel fiir ein derartiges Vorgehen sei 
die Arbeit von C. A. M. Gray} genannt, die eine spezielle Aufgabe fiir einen quadratischen 
Bereich lost. ; 


f) Naherungslisungen fiir Bereiche, die von mehrere nn Kieu 
turen begrenzt sind. Ist der betrachtete Bereich von einer endlichen Anzahl n von 
Konturen begrenzt, wie zum Beispiel die unendliche Ebene mit n Offnungen, so kann die 
Lésung durch ein alternierendes Verfahren angenahert erreicht werden. Es werden dabei 
zunidchst diejenigen Spannungsfunktionen bestimmt, die die vorgegebenen Randbedingungen 
jeweils an einer Kontur befriedigen wiirden, wenn die tibrigen nicht vorhanden waren. Die 
Summe dieser n Lésungen stellt die erste Naherung dar. Dabei hat jedoch jeder der n Anteile 
an den iibrigen n — 1 Konturen Randwerte zur Folge, die mit den Randbedingungen nicht ver- 
traglich sind. Die nachst bessere Naherung wird durch Addition derjenigen Spannungs- 
funktionen erhalten, die die stérenden Randwerte jeweils an einer Kontur zu Null machen 
wiirden, wenn die iibrigen n — 1 Konturen nicht vorhanden waren. Auch diese Lisung weist 
Abweichungen der Randwerte von den Sollwerten an allen Konturen auf, die jedoch geringer 
sind als die bei der vorhergehenden Naherung aufgetretenen, da die Spannungen, die infolge 
einer Anderung an einer Kontur auftreten, mit der Entfernung abnehmen und daher auf den 
ibrigen Konturen Abweichungen zur Folge haben, die kleiner sind als die vorgenommene 
Anderung. Es ist einzusehen, da dieses Verfahren nach einer geniigend groBen Zahl von 
Schritten zu einer Lésung von gewiinschter Genauigkeit fiihren wird. Die Konvergenz des 
Verfahrens ist um so besser, je gréBer der Abstand der einzelnen Konturen im Verhiltnis zu 
ihrer Ausdehnung ist. Der Rechenaufwand steigt mit der Anzahl der Konturen jedoch er- 
heblich an. 

Als einfaches Beispiel sei hier cin Parallelstreifen behandelt, der von den Geraden y =O 
und y = a begrenzt wird. Auf diesen Randern sei das Verschwinden der Tangentialspannung 
und der Normalverschiebung gefordert (dritte Randwertaufgabe). Im Innern des Streifens 
greifen an den Stellen (29,, ¥o,) (r = 1, 2,...n) die Krafte P, und Momente M. an. Da zur 
Aufstellung der Lésung ein haufiges Verschieben des Koordinatensystems erforderlich ist, er- 


1 C, A. M. Gray, Quart. Journ. Mech. appl. Math. 4 (1951), S. 444. 
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eis sich hier die hyperkomplexe Spannungsfunktion als vorteilhafter fiir die Aufstellung der 
Ssung. Die durch die vorgegebene Belastung im Streifen hervorgerufenen Singularitaten 


werden durch die Funkti i i wenn ; : 
ec ics see — beschrieben. Die erste Naherung ist nach (70) durch Spiegelung 


f(s) =fals) +Fo(s) + Sols — 2j a) 
zu erhalten. Kine Fortsetzung dieses Verfahrens bis ins Unendliche liefert dann die Lésung 


der gestellten Aufgabe: 


+ co 


Fs) = 2X thls +2 mja) + f(s +2 mj a} 


n=— 


oder mit (37) vollstandig ausgeschrieben: 


+ co n — 
a 1 —— 4 P, +-(1 +- 3?) (% PB, -- 3 P,) — 4 P, + (1+ j?) (% Pr +3 Pr) 
Zt) | Sr +2mja a s’ +2mja 


(79) 


_iMw 1 1 
2 (sr + 2 mja)? er imjar| 


unter Verwendung der Bezeichnungen -s, = 5 —5o,,° 3), = %r +i Yor und s&s, =s— so,, 
Sor = %r—J Yor- 

Das Ergebnis stellt eine einfach periodische Funktion dar, die sich auch wie folgt zusammen- 
fassen 1aBt: 


1 { “ = IU TELS 
© P 2 ae r 
f(s) Paw 1) I 4 ecm sad) ear oto Oe 2aj 
eee a7 P23 Py ae | 
_iM 2% 4 a2 j? 4 a2 j? 
= 2 pera Ie ; 7S) | 
rai} sin : pore | [ee 
(*.i) (Fes) 
Die (79) entsprechenden komplexen Spannungsfunktionen sind 
en ee < set P, 
SO 22 paces ieee falas | 
-+ oo n = 
ot | 1 | x Pr x P; 
#(2) a aiee | time . z +2ima 
(Zor-+2ima)P,  (@,+2ima) Pe iMie +1) , i M(x +1) 
(zr + 2ima)? (2, + 2ima)?* (2% +2ima)® | (2 + 2ima)? 


unter Verwendung der Bezeichnungen z, = 2— %,, %r = %or +tYor und ay == Z— Zor, 
en Kye iy¥or- Hier ist Y(z) keine einfach periodische Funktion, so daB eine Zusammen- 
fassung nach (80) nicht méglich ist. Die Ausdriicke fiir die Spannungen (7a) und (7b) werden 
jedoch, wie auch die hyperkomplexe Spannungsfunktion erwarten |aBt, einfach periodische 
Funktionen von z und z. 


6. Zusammenfassung und Schlufbetrachtung. Im ersten Teil der Arbeit werden die kom- 
plexen und hyperkomplexen Spannungsfunktionen als Mittel zur Lésung ebener Elastizitats- 
aufgaben eingefiihrt. In beiden Darstellungsarten werden die Bedingungen angegeben, die in 
mehrfach zusammenhangenden Bereichen zusitzlich erfiillt werden miissen, das Verhalten 
der Spannungsfunktionen gegeniiber einer Verschiebung des Koordinatensystems wird 
untersucht, und die Spannungskomponenten werden in krummlinig orthogonalen Koordinaten 
ausgedriickt. 

Es wird gezeigt, da8 sich in einfacher Weise ein in komplexer Darstellung vorliegendes Er- 
gebnis in die hyperkomplexe Form itberfithren 1a6t und umgekehrt, wobei die Gestalt des Er- 

'gebnisses nur wenig verdndert wird. Daraus kann geschlossen werden, da8B sich die beiden 


2* 
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Methoden nicht grundsitzlich unterscheiden und da auch der Arbeitsaufwand bei einer 
numerischen Auswertung entsprechender Formeln etwa der gleiche sein wird. 

Der zweite Teil behandelt die Anwendungsmiéglichkeiten. In komplexer Darstellung werden 
Methoden aufgezeigt, die ebene Elastizititsaufgaben von groBer Allgemeinheit, auch fir Be- 
reiche, die von mehreren beliebigen Konturen begrenzt sind, exakt oder in beliebig genauer 
Naherung zu lésen gestatten. Liésungen in hyperkomplexer Darstellung sind bisher nur fiir 
Bereiche, die von Geraden oder Kreisen begrenzt sind, auf denen bestimmte Randwerte zu 
Null werden, bekannt. In diesen Fallen wurde die hyperkomplexe Lésung der komplexen 
gegenitbergestellt, um die Ahnlichkeit der beiden Darstellungsarten erkennen zu lassen. 

Es wurde gefunden, da in manchen Fallen, insbesondere bei den alternierenden Lésungs- 
methoden fiir Bereiche, die von mehreren Kreisen oder Geraden begrenzt sind (wie zum Bei- 
spiel der Parallelstreifen) die hyperkomplexe Funktion fiir die Aufstellung der Lisung besser 
geeignet ist. 

In Ergainzung der hier gegebenen Anwendungsméglichkeiten sei darauf hingewiesen, daB 
beide Verfahren auch zur Liésung von Plattenaufgaben benutzt worden sind und daB dariiber 
hinaus bei entsprechender Anpassung der Formeln auch orthogonal anisitrope Scheiben oder 
Platten behandelt werden kénnen. Auf diesen Gebieten werden noch manche interessante 
Ergebnisse zu erwarten sein. 

Nach Abschlu8 dieser Arbeit wurde dem Verfasser ein Aufsatz von R. Legendre! bekannt, 
der sich mit dem Zusammenhang zwischen der hyperkomplexen und einer sogenannten 
doppelt komplexen Spannungsfunktion befaS8t und dhnliche interessante Zusammenhdange 
aufdeckt. Das angefithrte Beispiel (Kreisscheibe mit zwei im Innern aber nicht auf einem 
Durchmesser angreifenden Kraften) ist jedoch unvollstandig, da der Einflu8 der Querzahl 
auBer acht gelassen wurde. 


1 R. Legendre, La Recherche Aeronautique 27 (1952), S. 3. 


(Eingegangen am 18. April 1953.) 
Anschrift des Verfassers: Dr.-Ing. Kurt Stahl, (17a) Mannheim-Feudenheim, OdenwaldstraBe 12. 


XXII. Band 1954. Truckenbrodt: Ein Quadraturverfahren zur Berechnung der Reibungsschicht usw. 2] 


Ein Quadraturverfahren zur Berechnung der Reibungsschicht 
an axial angestrémten rotierenden Drehkérpern. * 


Von E. Truckenbrodt 


1, Einleitung. Mit der Berechnung der Reibungsschicht an Rotationskérpern, die sich 
um ihre Achse drehen und gleichzeitig in Richtung der Drehachse angestrémt werden, haben 
sich in neuerer Zeit H. Schlichting und E. Truckenbrodt befaBt. Zunichst wurde von beiden 
Verfassern der Fall der angeblasenen rotierenden Scheibe in laminarer Strémung behandelt }. 
Die mittels eines Naherungsverfahrens gefundenen Ergebnisse wurden inzwischen durch 
exakte Rechnungen von A. N. Tifford und Sheng To Chu? gut bestatigt. Mittlerweile wurde 
das Verfahren auch auf den Fall turbulenter Strémung iibertragen®. Das wesentliche Ergebnis 
der genannten Arbeiten besteht in der Berechnung des Drehmomentes der angeblasenen 
rotierenden Scheibe. Dabei hat sich gezeigt, da® dieses stark abhingig ist vom Verhaltnis 
der Umfangsgeschwindigkeit der Scheibe zur Anstrémungsgeschwindigkeit. Das Problem der 
angeblasenen rotierenden Scheibe liegt auch vor im Staupunkt des vorn stumpfen angeblasenen 
rotierenden Drehkérpers. Dieses sieht man sofort ein, wenn man bedenkt, da® man die Um- 
gebung des Staupunktes durch seine senkrecht zur Anstrémung stehende Ebene ersetzen kann. 


Der allgemeine Fall eines angeblasenen rotierenden Drehkérpers in laminarer Strémung ist 
von H. Schlichting* bearbeitet worden. Diese Arbeit stellt eine entsprechende Erweiterung 
des fiir den Fall ebener Strémung bekannten Verfahrens von K. Pohlhausen® und H. Holstein 
und T. Bohlen® dar. Die Lésung wird zuriickgefihrt auf zwei miteinander gekoppelte simul- 
tane Differentialgleichungen. Aus den Ergebnissen ist wieder die starke Abhangigkeit des 
Drehmomentes vom Verhaltnis der Umfangsgeschwindigkeit zur Anstrémungsgeschwindigkeit 
zu erwahnen. 


Der noch fehlende Fall des angeblasenen sich drehenden Rotationskérpers in turbulenter 
Strémung ist Gegenstand der vorliegenden Untersuchung. Die hierbei gewonnenen Formeln 
ergeben bei sinnvoller Ubertragung auf den Fall laminarer Strémung auch hierfiir recht 
brauchbare Naherungsformeln. 


Experimentelle Untersuchungen zu dem behandelten Problemkreis wurden durchgefihrt 


von C. Wieselsberger’ und S. Luthander und A. Rydberg ®. 


Ein unserer Aufgabe ahnliches Problem stellt der schiebende Fliigel unendlicher Spann- 
weitenerstreckung dar. Seine theoretische Berechnung inshesondere bei turbulenter Strémung 
wurde erstmalig von A. D. Young und T. B. Booth® durchgefiihrt. In Abb. 1 haben wir unser 
Problem des angeblasenen rotierenden Drehkorpers demjenigen des schiebenden Fliigels un- 
endlicher Spannweite gegeniibergestellt. Fiir den Teil des Drehkérpers, der nahezu zylindrisch 
ist, stellt die Oberseite des Fliigels die abgewickelte Oberfliche des Drehkérpers (gestrichelte 


* Die vorliegende Arbeit stellt den zweiten Teil einer von der Fakultat fiir Maschinenwesen der Tech- 
nischen Hochschule Braunschweig angenommenen Habilitationsschrift mit dem Titel ,,Die Str6mung um 
rotierende Drehkérper bei axialer Anstroémung‘ dar. (Berichter Prof. Dr. H. Schlichting, Mitberichter 
Prof. Dr. H. Blenk.) Der erste Teil mit dem Titel ,,Die turbulente Strémung an einer angeblasenen 
rotierenden Scheibe‘‘ wird in der,,Zeitschrift fiir angewandte Mathematik und Mechanik“ erscheinen. 

Diese Untersuchungen wurden mit Unterstiitzung der Deutschen Forschungsgemeinschaft im 
Institut fiir Str6mungsmechanik der T. H. Braunschweig ausgefibrt. 

1 H. Schlichting u. E. Truckenbrodt, Z. angew. Math. Mech., 32 (1952), 5.97; vel. auch Journ. 
Aeron. Sc. 18 (1951), S. 638. 

2 A. N. Tifford u. Sheng To Chu, Journ. Aeron. Se. 19 (1952), 5. 284. 

8 E. Truckenbrodt, Die Strémung an einer angeblasenen rotierenden Scheibe bei turbulenter Stré- 
mung; erscheint demnachst in der Z. angew. Math. Mech. 

4 H. Schlichting, Ing.-Arch. 21 (1953), 5, 227. 

5 K. Pohlhausen, Z. angew. Math. Mech., 1 (1921), 8. 252. ; 

6 H. Holstein u. T. Bohlen, Ein einfaches Verfahren zur Berechnung laminarer Reibungsschichten, 
die dem Naherungsansatz von K. Pohlhausen geniigen, Lilienthal-Bericht 5. 10 (1940), 5. 5. 

7 C, Wieselsberger, Physik. Z. 28 (1927), 5. 84. 

8 §, Luthander u. A. Rydberg, Physik. Z. 36 (1935), S. 552. 

® A.D. Young u. T. B. Booth, Aeronaut. Quarterly, 3 (1951/52), Sy PALATE 
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Fliche) dar. Die Geschwindigkeitskomponente in Spannweitenrichtung V,, entspricht der 
Umfangsgeschwindigkeit vy = @ R, und die Geschwindigkeitskomponente in Sehnenrichtung 
U der Geschwindigkeitskomponente in meridionaler Richtung U am Drehkérper. 


_ 2. Die Grundgleichungen. Die Ausgangs- 
gleichungen fiir unser Verfahren bilden die 
Prandtlschen Grenzschichtgleichungen. Diese 
gelten unter der Annabme, da die Umdrehungs- 
geschwindigkeiten nicht allzu groB im Ver- 
gleich zur Anstrémungsgeschwindigkeit sind. 


Abb. 1. Gegeniiberstellung eines angestrémten sich drehenden Abb. 2. Erlauterungsskizze: 
Rotationskérpers (a) mit einem schiebenden Kérper unend- Geschwindigkeiten x: u, U; y: v, » = oR; 2: w. 
licher Spannweite (6). Schubspannungen *: T,3 y: Ty: 


Lediglich fiir den Sonderfall der angeblasenen rotierenden Scheibe 1a8t sich mit den Grenz- 
schichtmethoden der gesamte Bereich von der angeblasenen stillstehenden Scheibe bis zur 
nicht angeblasenen rotierenden Scheibe erfassen. Ist die Reibungsschichtdicke 6 klein im 
Vergleich zum Kérperhalbmesser, dann gilt fiir den angeblasenen rotierenden Drehkérper 


in meridionaler Richtung 


du du dU vdR 1 0x 
ies ee eta 
Uapeiailias Oe Srdx somo: ” 
und in azimutaler Richtung 
dv dv wovdR il Ge 
= | 2, 
Pee as Reuss (2) 
Die Kontinuitatsgleichung lautet 
0(u R) O(wR) 
Ox i: Oz us (3) 
Die Randbedingungen sind 
2302005 0 Syn, w= 0s ee 0 ee (4) 


Die Bedeutung der einzelnen GréBen geht aus Abb. 2 hervor. Es verlaufen x in meridionaler 
Richtung, y in azimutaler Richtung (Umfangsrichtung) und z senkrecht zu diesen beiden 
Richtungen. R(x) stellt die Abhingigkeit des Kérperhalbmessers von der Lauflange x dar. 
Die Geschwindigkeitskomponenten innerhalb der Reibungsschicht in Richtung der drei 
Achsen x, y und z werden mit u, v und w bezeichnet. U(x) ist die auBerhalb der Reibungs- 
schicht (z > 6) vorgegebene Potentialgeschwindigkeit. Die GroBe v) (x) =a R gibt die von 
« abhangige Umfangsgeschwindigkeit an, wobei w die Winkelgeschwindigkeit ist, mit der der 
Kéorper um seine Langsachse dreht. Die Bezeichnungen 1, und Ty stellen die Schubspannungs- 
komponenten in «- und y-Richtung dar, und 9 ist die Dichte des strémenden Mediums. 

Durch Integrationen der Gleichungen (1) und (2) itber den Wandabstand von z — 0 bis 
z = 06 (6 = Grenzschichtdicke) und Beriicksichtigung der Kontinuitatsgleichung (3) ergeben 
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sich zwei Impulsgleichungen in meridionaler und azimutaler Richtung. Die von H. Schlich- 
ting * abgeleiteten Impulsgleichungen gelten sowohl fiir 
die turbulente Strémung, 
Schreibweise 


gen die laminare Strimung als auch fiir 
so daB wir sie sofort ttbernehmen kinnen. Sie lauten in etwas anderer 


fiir die meridionale Richtung 


z 2 
Ua ja” oe =) Dy dR to 


UR dx U de '\U) R dx ~ 9 (5) 
und fiir die azimutale Richtung 
. ad 1 
bs af ae y 
mp Ridx\UR Ixy) sane te (6) 
Hierin bezeichnen 
3 
u u 
i, : (1 : dz (7) 
f 
3 
t= er dz. Impulsverlustdicken, (8) 
F 0 
@. 
uv 
8 
s* = (1 " a) dz  Verdrangungsdicke, (10) 
0 
y R 
a a= a Drehparameter (11) 


und Ty und ty, die Wandschubspannungskomponenten in meridionaler und azimutaler 
Richtung. 

Die Formeln (5) und (6) stellen die Grundgleichungen unserer Aufgabe dar. Um sie an- 
wenden zu kénnen, sind jedoch noch einige Angaben iiber die Wandschubspannungskomponen- 
ten erforderlich. Hieriiber wird im nachsten Abschnitt berichtet werden. 


3. Die Ansatze fiir die Wandschubspannungskomponenten. Wahrend fiir die laminare 
Strémung ein fester Zusammenhang zwischen den Wandschubspannungskomponenten und 
den Geschwindigkeitsgradienten an der Wand in der Form 


Ou Ov 


besteht, liegen die Verhaltnisse bei der turbulenten Strémung nicht so einfach. Wir sind hierfiir 
wie bisher bei allen turbulenten Strémungen auf Naherungsansatze angewiesen. Da es sich bei 
unserem Problem um eine dreidimensionale Reibungsschicht handelt, deren Verhalten im 
Einzelnen nicht ganz einfach zu iibersehen ist, gehen wir bei der Betrachtung so vor, daB wir 
zunachst den nahezu zylindrischen Teil des Drehkérpers, dR/dx ~ 0, behandeln. Wie wir 
bereits in der Einleitung zeigten, verhalt sich dieser Teil ahnlich wie ein schiebender unendlich 
langer Fliigel, so da8 wir die hierfiir bekannten Ergebnisse ibernehmen kénnen. Daran an- 
schlieBend untersuchen wir den Fall, daB die Oberflache des Drehkérpers kegelig ist, dR/dx--0. 
Die stairkste Anderung des Kérperhalbmessers mit der Lauflainge ergibt sich im Staupunkt, 
wo dR/dx = 1 wird. Wie wir schon in der Einleitung erwahnten, ist die Lésung hierfiir bereits 
durch die Untersuchungen an der angeblasenen rotierenden Scheibe bekannt. Wir werden die 
Wandschubspannungskomponenten in meridionaler und azimutaler Richtung im folgenden 
getrennt voneinander behandeln. 


a) Die Wandschubspannungskomponente in meridionaler 


Richtung. Die Untersuchungen von A. D. Young und T. B. Booth® sowie N. Rott und 


1 Siehe FuBnote 4 von S. 21. 
2 Siehe FuBnote 9 von S. 21, 
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L. F. Crabtree} am schiebenden unendlich langen Fliigel haben ergeben, da die Wandspan- 
nungskomponente in Sehnenrichtung von der Querstrémung in Spannweitenrichtung unab- 
hingig ist. Auf unsere Verhiltnisse ithertragen bedeutet das, daB die Wandschubspannungs- 
komponente in meridionaler Richtung fiir den Teil des Drehkérpers, der nahezu zylindrisch 
ist, von der Drehung des Kérpers nicht abhangt: 

Tx Tx =i 

ae = (a), (R =konst.). (13) 
Unter dieser Annahme ergibt sich aber, wie wir sofort aus Gleichung (5), in der das letzte Glied 
auf der linken Seite verschwindet, ersehen, daB® die Impulsverlustdicke #, und die Ver- 
dringungsdicke 6* durch die Drehung unverandert bleiben. 


Ist die Oberflache des Drehkérpers jetzt kegelig, so kommt noch die Wirkung der Zentri- 
fugalkraft auf die Reibungsschicht hinzu. Nehmen wir auch hierfiir weiterhin an, da8B sich die 
Impulsverlustdicke und die Verdrangungsdicke nicht mit der Drehung aindern?, dann folgt 
aus der Impulsgleichung (5) fiir die Wandschubspannungskomponente in meridionaler Richtung 
sofort 


Tx Tx U\2 Oy dR 
8 Gale (a) | (R+ konst.). (14) 


DaB dieser Ansatz auch die Verhaltnisse im Staupunkt, wo dR/dx = 1 ist, gut wiedergibt, 
kann man durch Vergleich mit der Wandschubspannungskomponente in radialer Richtung 
bei der angeblasenen rotierenden Scheibe zeigen. Fir maBige Umfangsgeschwindigkeiten 
im Vergleich zur Anstrémungsgeschwindigkeit kann das aus der Zentrifugalwirkung in (14) 
herriihrende Glied vernachlassigt werden. 


b) Die Wandschubspannungskomponenteinazimutaler Rich- 
tung. Von Young-Booth und Rott-Crabiree werden keine Angaben iiber einen Ansatz fir 
die Wandschubspannungskomponente in Spannweitenrichtung des schiebenden unendlich 
langen Fligels gemacht, so da8 wir hierauf nicht zuiiickgreifen kénnen. Fir den nahezu 
zylindrischen Teil des Drehkérpers wollen wir setzen R=konst. und U = konst. Aus den 
Impulsgleichungen (5) und (6) folgt dann 

dx Tx Ixy Vy Tyo 


PE a ENE d Set 
dz (pU" ) oe dx U ov 


(R =konst.) (15) 


Wenn wir jetzt annehmen, daB zwischen den Geschwindigkeitsverteilungen in meridionaler 
und azimutaler Richtung Ahnlichkeit besteht, derart, daB man v/v) = 1 — u/U setzen kann, 
dann ergibt sich aus (7) und (9) Gleichheit fiir die Impulsverlustdicken #,, = 0,. Unter 
Beachtung dieses Ergebnisses erhalten wir die folgende Beziehung zwischen den beiden Wand- 
schubspannungskomponenten 


Tyo = — Gtx (R =konst.). (16) 


Als Ergebnis haben wir also gefunden, daB sich die Wandschubspannungskomponenten in 
azimutaler und radialer Richtung verhalten wie die Umfangsgeschwindigkeit zur Anstrémungs- 
geschwindigkeit. 


Fiir die Wandschubspannungskomponente 1) machen wir ausgehend von den Gesetz- 
maBigkeiten der langsangestrémten ebenen Platte den bekannten Ansatz 
BLE OG i 
aE eo o 
Noy 
worin fiir den laminaren Fall n = 1 und « = 0,220 und fiir den turbulenten Fall n — 1/6 und 
«% = 0,0065 zu setzen ist. Diese Beziehung setzen wir in (16) ein und erhalten, wenn wir noch 
%, = 0, beriicksichtigen, fiir die Wandschubspannungskomponente in Umfangsrichtung 


eee 
ov eee : (18) 
v 


1 N. Rott u. L. F. Crabtree, Journ. Aeron. Sc. 19, (1952), S. 553. 


as Mer wird fiir die laminare Strémung durch die Rechnungen von H. Schlichting sehr 
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Diese Gleichung gilt entsprechend den oben gemachten Angaben zunichst fir den Fall, 
daf R=konst. und U=konst. sind. Wir wollen sie sinngemaf auf den Fall itbertragen, daB 
die Kérperoberflache kegelig ist, R+= konst. und U-+: konst., indem wir fiir R und U die von x 
abhangigen Werte einsetzen. Eine Begriindung hierzu sei noch im folgenden gegeben: Nach 
der Impulsgleichung (6) ist die Wandschubspannungskomponente abhangig von der Reynolds- 
schen Zahl U #,,/y und von der Lauflange x. Die von x abhangigen GréBen sind die Verteilung 
des Kérperhalbmessers R(x) und die Geschwindigkeitsverteilung U(x) in meridionaler Rich- 
tung. Wenn wir annehmen, daf zwischen diesen beiden GréBen ein eindeutiger Zusammenhang 
besteht, dann kénnen wir auch schreiben 
Tyo U%sy U 
gt e(9 mp” 
Diese Abhangigkeit wird von (18) richtig wiedergegeben, so da8 wir sie auch fiir den allgemeinen 
Fall als giiltig ansehen wollen. Da8 (18) auch fiir die Staupunktstrémung brauchbar ist, zeigt 
wieder ein Vergleich mit den Ergebnissen an der angeblasenen rotierenden Scheibe. 

Den Wert a in (17) bzw. (18) kénnen wir in Beziehung setzen zum Widerstandsheiwert CF 
einer mit der Geschwindigkeit U,, langsangestroémten ebenen Platte der Lange 1. Man findet 


dabei den Ausdruck 
va men (th 


der sich aus dem folgenden allgemein geschriebenen Widerstandsgesetz ergibt: 


n 


(19) 


ary lin 
= es } (21) 
Im laminaren Fall ist n = 1 und c = 1,328. Im turbulenten Fall gilt nach Falkner n = 1/6 


und c = 0,0306. Will man mit der alteren Formel von Prandtl rechnen, dann ist n = 1/4 und 
c = 0,074. 


4, Die Berechnung der Impulsverlustdicken. Wir setzen jetzt die Ergebnisse des vorigen 
Abschnittes in die Impulsgleichungen (5) und (6) ein und leiten daraus einfache Quadratur- 
formeln zur Berechnung der Impulsverlustdicken #, und 3, , ab. Wegen der Ansatze fiir die 
Wandschubspannungskomponenten Tt, und ty, nach (14) und (18) werden die beiden Impuls- 
gleichungen voneinander unabhangig, so daB wir sie im folgenden getrennt voneinander be- 
handeln kénnen }. 


a) Die Impulsverlustdicke #,. Setzen wir (14) in (5) ein, dann erhalten wir 
fir die Impulsgleichung in meridionaler Richtung 


1 d Ox dU Txo 

UR FN Bea U dx ot? laa (22) 
Unter den gemachten Annahmen fiir die Wandschubspannungskomponente in meridionaler 
Richtung folgt, daB die Impulsgleichung in meridionaler Richtung unabhangig von der Drehung 
des Kérpers wird. Zur Lisung von (22) kinnen wir daher unmittelbar die Formeln anwenden, 
die der Berechnung der Impulsverlustdicke # des nicht drehenden Kérpers dienen. Wir be- 
nutzen hier die Quadraturformel von E. Truckenbrodt * 


x/I 
c,\l+n U \3+2n7/R\1l+n x 
ler(e) flee) ir) 
; ! 


= = a/ : (23) 


1 
lt+n 


Hierin bedeuten U,, die Anstrémungsgeschwindigkeit des Kérpers, | seine Lange (oder eine 
sonst anderswie definierte Lange) und cy den Reibungsheiwert einer mit der Geschwindigkeit 
U,, angestrémten ebenen Platte der Linge J. Bei laminarer Strémung ist n = 1 und bei 
turbulenter Strémung n = 1/6 zu setzen. 


1 In dem Verfahren von H. Schlichting, FuBnote 4 von S. 21, fiir den laminaren Fall bleibt die Kopplung 
beider Impulsgleichungen dadurch erhalten, daB die Reibungsschichtdicken in meridionaler Richtung Ox 
und in azimutaler Richtung dy voneinander verschieden angenommen werden und zu den beiden Unbe- 
kannten des Problems werden. 

2 EB. Truckenbrodt, Ing.-Arch. 20 (1952), S. 211; vgl. auch Journ. Aeron. Sc. 19 (1952), S. 428. 
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Die Integrationskonstante Cj folgt aus dem laminaren Wert fiir die Impulsverlustdicke 0, 
im Umschlagspunkt x = x,: iS 
n 


* U\3R ox) *” ef in Uae Nee oe ae (24) 
welbalcube. pauliestial | 4p 
: 0 


Bei voll-laminarer bzw. voll-turbulenter Strémung ist mit x, = 0 auch Cj = 0. 


b) Die Impulsverlustdicke #,y. Aus (6) erhalten wir mit (18) 


GT Roh ys UE 25 

7, (UR Dx y) = (een UR’. ( ) 
v 

Wie man leicht zeigen kann, laBt sich diese Differentialgleichung geschlossen lésen. Man macht 

hierzu zweckmabig die Substitution 


y= U R30+%) ea Bey - 
Als Ergebnis der Integration erhalten wir 


K,+(1+n)af UR 


U Oxy ug ee % “ 26 
( y ze U R38 (+) ; (26) 


Die GréBe (1+ n)a driicken wir jetzt nach (20) durch den Reibungsbeiwert cy einer mit der 
Geschwindigkeit U,, angestrémten ebenen Platte der Linge J aus. Dann ergibt sich nach #, y 
aufgelést 

1 


CASS [ HURT RASC ee 
a 3 Sats Se eve Mtge 
at (2 an (ah eg 


Duy x,/l 


l U /R\3 
Geer 


Bei laminarer Strémung ist n = 1, bei turbulenter Strémung n = 1/6 zu setzen. 
Die Integrationskonstante Kf bestimmt sich wieder aus dem laminaren Wert fiir die 
Impulsverlustdicke #,,. im Umschlagspunkt x = xj. 


(27) 


ot al 1/2]1-+n 
se (RNG Gey es loa U /(R\6 , x 
= Lo. (7) ree ; Wray se 
0 


Bei voll-laminarer bzw. voll-turbulenter Strémung ist mit x, = 0 auch Kj = 0. 


Bemerkenswert an der Formel (27) ist, daB die Impulsverlustdicke #, , unabhangig von der 
Drehgeschwindigkeit ist. Auch dieses Ergebnis wurde fiir den laminaren Fall von H. Schlich- 
ting gefunden. 

Vergleicht man die beiden Formeln fiir die Impulsverlustdicken (23) und (27), so stellt man 
fest, daB bei der Impulsverlustdicke in meridionaler Richtung #,, die Geschwindigkeitsvertei- 
lung starker eingeht als bei der Impulsverlustdicke },,. Beziiglich der Verteilung des Kérper- 
halbmessers liegen die Verhaltnisse gerade umgekehrt. 


5. Die Abschatzung der Lage des Ablésungspunktes, Durch die Drehbewegung wird die 
Lage des Ablésungspunktes verindert werden. Eine genauere Berechnung ist jedoch zur Zeit 
noch nicht méglich, da selbst im Fall der ruhenden Wand das Kriterium fiir die Ablésung nur 
angenadhert bekannt ist. Wir miissen uns deshalb damit begniigen, die Lage des Ablésungs- 
punktes grob abzuschatzen. Zu einer genaueren Ermittlung miiBte man die Verteilung der 
Schubspannungen iiber die Reibungsschichtdicke kennen, um daraus dhnlich wie fiir den Fall 
des nicht drehenden Kérpers eine Gleichung fiir den das Geschwindigkeitsprofil in der 
Reibungsschicht kennzeichnenden Formparameter herzuleiten. 


Wir wollen ahnlich wieim laminaren Fall vorgehen. Dort wird als charakteristischer Parameter 


7, U0e, Tl dU) Vay eta 
af y OG ae R d| (29) 
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eingefiihrt. Erreicht x*—x4 einen bestimmten Wert (nach dem Pohlhausen-Verfahren 
x4 =— 0,1567), dann tritt Ablésung ein. 


Man kann sich iitberlegen, da im turbulenten Fall ein entsprechender Parameter folgender- 


mafen lauten miifte: 
_ [U@s\» . [1 dU , /m\21 dR 
L ( v Bs E dx +(¥] R ae" (30) 


Fir den Fall ohne Drehung (v) = 0) geht I’ mit n = 1/4 in den von Buri! verwendeten 
Parameter tiber. Wie im laminaren Fall nehmen wir nun an, daB Ablésung eintritt, wenn J" 
einen bestimmten Zahlenwert erreicht hat. Diesen zunadchst noch unbekannten Wert jie 
bestimmen wir aus den Verhaltnissen des nicht drehenden Kérpers zu 


Ud\" bx dU 
(reg. 0 


Hierin bedeutet x4, die Lage des Ablisungspunktes fiir v) = 0. Fiir die Impulsverlustdicken 
%, nehmen wir aus dem im vorigen Abschnitt Gesagten an, daB sie von der Drehgeschwindig- 
keit unabhangig sind. 

Die Lage des Ablésungspunktes x4) ermitteln wir z.B. nach dem Verfahren von E. 
Truckenbrodt?. Es soll an dieser Stelle auf Einzelheiten einer solchen Rechnung nicht naher 
eingegangen werden. 

Die Bestimmungsgleichung fiir die Lage des Ablésungspunktes lautet jetzt nach (30) und 


(31) 


n 


i! 
Ox / U \l4n/ 1 dU\14n 
l te) G i) xAQ 


n 1 
U\ttepl dU , /m\? 1 dR]i+= 

(az) la page (B) R aa chy 
Nach dieser Beziehung berechnet man fiir verschiedene Werte x4 die Impulsverlustdicken 
@,/1 und bringt diese Kurve zum Schnitt mit der Kurve fiir 3/1 nach (23). Der Schnittpunkt 
liefert dann den gesuchten Abszissenwert x4. Man vergleiche hierzu das Beispiel der Kugel, 
Abb. 4. | 

Da im Druckanstiegsgebiet, wo Ablésung auftreten kann, sowohl dU/dx als auch dR/dx 
gewohnlich negativ sind, folgt, daB im allgemeinen 0,(%4) < %x(%49) wird. Das bedeutet 
aber auch, daB die Ablésung beim drehenden Korper friiher eintritt als beim nicht drehenden 
Kérper. Auch dieses Ergebnis wird fiir den laminaren Fall von H. Schlichting bestatigt. 


(32) 


6. Die Quadraturformel zur Berechnung des Drehmomentes. Die Berechnung des durch 
die Drehbewegung entstehenden Drehmomentes |4Bt sich in einfacher Weise durchfiihren. Es 
braucht hierbei lediglich tiber die in Umfangsrichtung wirkenden Schubspannungskomponenten 
Tyo integriert zu werden. 

Das gesamte Drehmoment erhalt man zu 


M =— 2m [ ty, Redx. (33) 
0 


Es ist x, die Bogenlange des Korpers in meridionaler Richtung vom Staupunkt aus gemessen. 
Lost die Reibungsschicht ab, so ist der Punkt, bis zu dem die Integration auszufihren ist, 
nicht mehr klar definiert. Wie man aber aus den durchgerechneten Beispielen ersieht, spielt 
die GréBe x, keine entscheidende Rolle, so daB man keinen grofBen Fehler macht, wenn man bis 
zum Kérperheck die Integration erstreckt. 

Die Schubspannungskomponente ty 9 entnehmen wir der Impulsgleichung fiir die azi- 
mutale Richtung (6) mit » =@ R zu 


d 
Tyo R? =-—oo7 (UR Oiser)s (34) 


1 —, Buri, Eine Berechnungsgrundlage fiir die turbulente Grenzschicht bei beschleunigter und ver- 
zégerter Stromung, Dissertation, Ziirich (1931); vgl. H. Schlichting, Grenzschicht-Theorie, 5. 422, Karls- 
ruhe (1951). 

2 Siehe FufSnote 2 von S. 25. 
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Diesen Wert setzen wir in (33) ein und finden nach Ausfiihren der Integration 


M=2200(U R98, ,)(c=2,) - (35) 
Einen Momentenbeiwert definieren wir folgendermaSen: 
fi) Gs aioe eR (36) 
& v2 Ril 


Hierin ist R,, der gréBte Kérperhalbmesser (Radius des Hauptspantquerschnittes), Vi, = @ Rn 
die zugehérige Umfangsgeschwindigkeit und 1 eine beliebig definierte Kérperlange. Die Defini- 
tion fiir cy geht fiir die rotierende Scheibe mit 1 = R,, in die frither angegebene Formel 
iiber. 
In die Momentenformel (35) setzen wir jetzt die Formel (27) fiir die Impulsverlustdicke 0, y 
ein und finden fiir den Momentenbeiwert nach (36) 
1 


ail ee 
Ue é\lgn | U /R\80+a) , 412° 2) 
iar x Eta ole lra| uy Ct 
jl 


%y 


Die Integrationskonstante ergibt sich aus dem laminaren Wert fiir die Impulsverlustdicke 0, y 
an der Stelle x, zu 


x,/1 om 
joe U /[ R\3 Oxy I+n a Cry "U a ee Viziie: (38) 
t= [ea (e) Flew [2 Af Balke) ST) J 
Uber die Widerstandsheiwerte cy der ebenen Platte als Funktionen der Reynoldsschen Zahl 
U.l/v wurde bereits berichtet, [S. 25].- 
Aus (37) folgt nun sofort, daB sich der Momentenbeiwert bei kon stant gehaltener 
Anstrémungsgeschwindigkeit U, umgekehrt proportional der Umfangs- 


geschwindigkeit andert. Weiterhin findet man aber unter Beachtung von (36), daB das Dreh- 
moment linear mit der Umfangsgeschwindigkeit anwachst. 


M ~ Vm (Us = konst.) . (39) 


Bei voll-laminarer bzw. voll-turbulenter Strémung verschwindet mit x, = 0 auch K{ = 0 in 
(37). Der Drehmomentenbeiwert verhalt sich dann folgendermafen: 


at n 


Us (USING len Uso\1+2 [Vin Rm\ 147 
ace cs) a 


Bei konstant gehaltener Umfangsgeschwindigkeit folgt dann fir 
den Momentenbeiwert sowie das Moment selbst: 
1 


M~U‘**  —(V,,= konst.) . (41) 


Bei laminarer Strémung (n = 1) nimmt das Moment mit der Wurzel aus der Anstrémungs- 
geschwindigkeit zu. Bei turbulenter Strémung (nm = 1/6) wachst es nahezu linear mit der 
Anstrémungsgeschwindigkeit. Das Ergebnis im laminaren Fall ist wieder in Ubereinstimmung 
mit den Angaben von H. Schlichting. 

Auf eine Berechnung des Strémungswiderstandes (Reibungs- und Druckwiderstand) wird 
verzichtet, da hierfiir die Entwicklung des Nachlaufes hinter dem drehenden Rotationskérper 
bekannt sein miiBte. Hierfiir fehlen aber noch jegliche Angaben sowohl theoretischer wie auch 
experimenteller Art. Eine weitere wichtige Tatsache ist aber die, daB gerade diejenigen Kérper, 
wie z. B. Geschosse, von denen man den Widerstand kennen méchte, am hinteren Ende scharf 
abgeschnitten sind, wodurch sich ein stark abgeléstes Totwasser hinter dem Kérper aus- 
bildet. Dadurch wird aber der Widerstand entscheidend bestimmt, wie auch aus einer Nahe- 
rungsbetrachtung von C. Wieselsberger1 hervorgeht. 


7. Beispiele. a) Rotierende Scheibe. Wie bereits in Ziff.1 erwaihnt wurde, 
liegen fiir die angeblasene rotierende Scheibe sowohl im laminaren wie auch turbulenten Fall 
theoretische Rechnungen vor. 


1 Siehe FuBnote 7 von S. 21, 
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Potentialtheoretisch verhalt sich die angeblasene rotierende Scheibe wie die rotations- 
symmetrische Staupunktsstrémung. Hierfiir gilt in bekannter Weise 


Cah bo se IK (he (42) 


Setzen wir in (37) x; = 0, x, = R, und Kj = 0, so liefert die ausgefiihrte Integration unter 
Beachtung des Widerstandsgesetzes (21) 


1 n 
l+n 2 T 1a 
om aly a 
vD Vv 
mit 
Uae ies 
Pee od Fete By eas (44) 


Setzen wir die Zahlenwerte fiir c und n ein, dann erhalt man fiir den laminaren Fall A = AO 5a 
wahrend sich fiir die asymptotische Lésung der Wert A = 3,17 ergibt!. Fiir den turbulenten 
Fall wollen wir mit der Prandtlschen Widerstandsformel, d.h. mit n = 1/4 rechnen, um mit 
dem friiheren Ergebnis? vergleichen zu kénnen. Wir erhalten A=0,115, wahrend sich fiir die 
asymptotische Formel A =0,113 ergibt. Unter der asymptotischen Lésung verstehen wir den 
Fall, bei dem die Umfangsgeschwindigkeit im Vergleich zur Anstrémungsgeschwindigkeit 
klein ist. Dieses gilt etwa bis zu Werten von a/w > 1 fir laminare Strémungen und a/w > 2 
fiir turbulente Strémungen. 


Unsere Untersuchung hat uns nun gezeigt, da8 unsere Momentenformel (37) die angeblasene 
rotierende Scheibe recht befriedigend wiedergibt, wenn wir uns auf maGige Verhaltnisse von 
Umfangsgeschwindigkeit zu Anstrémungsgeschwindigkeit beschranken. 


b) Rotierende Kugel. Als zweites Beispiel wollen wir die angeblasene rotierende 
Kugel wahlen und hierfiir eine vollstandige Rechnung nach dem beschriebenen Naherungsver- 
fahren durchfiithren. Wir wollen uns dabei auf die Falle voll-laminarer und voll-turbulenter 
Strémung beschranken. 

Ist R,, der Kugelhalbmesser, dann gilt fiir die vom Staupunkt aus gemessene Lauflange 


x=R,gp (0<9 <n), (45) 
und fiir die Verteilung des senkrecht zur Anstrémungsrichtung gemessenen Kérperhalbmessers 
Ri=.Ry, sin.@ . (46) 


Als Bezugslange 1 in unseren bisherigen Formeln wollen wir setzen 1 = R,,. Die potential- 
theoretische Geschwindigkeitsverteilung erhalt man in bekannter Weise zu 


U =+ U, sing. (47) 


Fir die Impulsverlustdicken #, und #, , nach (23) und (27) kann man sofort zeigen, dafs sie sich 
wegen der Proportionalitat von U ~ R nicht voneinander unterscheiden. 
Bei voll-laminarer und voll-turbulenter Strémung setzen wir x, = 0 sowie Cj] = Kj = 0. 


Wir erhalten dann 
1 


5 P USts 
| [sin 3" ag] 


a en De wis Oxy 2 eo Cr = 
Rm Rn Rae 3 2 sin’ @ 


n 


(48) 


Im laminaren Fall 1aBt sich die Integration mit n = 1 geschlossen ausfiihren. Auf die Wieder- 
gabe dieser etwas umstandlichen Formel wollen wir verzichten. Das Ergebnis der Rechnung 
ist in Tabelle 1 angegeben und in Abb. 3 dargestellt worden. Die Ubereinstimmung mit der 
genaueren Kurve fiir 3,/R, von H. Schlichting, die wir miteingetragen haben, ist sehr gut °. 


1 Der Lésung in der in FuBnote 1 von S. 21 genannten Arbeit liegt das Pohlhausen-Polynom fiir die 
Geschwindigkeitsverteilung zugrunde. Hierfiir ist ¢ = 1,372. In (44) eingesetzt wiirde sich dann ergeben 
PA = 3,00. 

2 Siehe FuBnote 3 von S. 21. - o 

8 Beachten wir das in der obigen Fufinote Gesagte, dann miifiten, um den Vergleich wouseag ae 
durchzufiihren, unsere Werte noch mit 1,372/1,328 = 1,033 multipliziert werden, was die gute Uberein- 
stimmung noch verbessern wiirde. 
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Fir den turbulenten Fall haben wir die Integration graphisch ausgefithrt und das Ergebnis 
in Tabelle 1 sowie Abb. 4 dargestellt. Nach dem Verfahren von E. Truckenbrodt 1 haben wir 
fiir die nicht drehende Kugel noch das Reibungsschichtdickenverhialtnis ASO e; perce 
(Tabelle 1). Wir haben dabei die Reynoldsschen Zahlen Re = U,, Rn/y = 10° und 10 gewahlt. 
Nimmt man an, daB® die Reibungsschicht sich ablést, wenn H =1,8 wird, dann gilt 


Ay = 126,5° fiir Re = 10° und bei 49 = 132,0° fir Re = 10’. 


7 wo 
Vie Foe 
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06 G08 


OS GOS 
04 GUY 
03 903 Un 
UE 
G2 | i 0} — 10 
— — Re 2 nach tt Schlichting 
Me 45 
a7 { oe got 
P _ Ai) fe | 
0 20° 40° 60° 80° 0° «ea? 0 20° 40° 60° 80° 700° 700° 109° 
Abb. 3. Die Impulsverlustdicken #, und Oxy fiir die Abb. 4. Die Impulsverlustdicken #, und Oxy fiir die Kugel bei 
Kugel bei laminarer Strémung; turbulenter Strémung; 
6=6,=%,,, Ro= UoR |, Wm, = ORp- P= 6, = 8,y, Re=Uokm/% Vm = ORm- 


JI] Ablésungs beginn. 


Tabelle 1. Die Impulsverlustdicke 0 = 0, = Oxy bei laminarer und turbulenter Strémung sowie der Form- 
parameter H bei turbulenter Strémung (vy = 0) fiir die Kugel; Re = U., Rm. 


| laminar | turbulent 
x R U —.  ¢ 7— @ H 
ae pats eee te, / i Se ii a 
oO Rn Rm c= IEE OME Re = 10° 107 
Om 0,000 0,000 0,192 0 
OS 0,174 0,260 
20° 0,342 0,513 
30° 0,500 0,750 0,202 0,0020 
40° 0,643 0,964 0,212 0,0027 
50° 0,766 1,149 0,226 0,0035 
60° 0,866 1,299 0,244 0,0044 
70° 0,940 1,410 0,271 0,0054 
80° 0,985 1,477 0,309 0,0068 
90° 1,000 1,500 0,366 0,0088 1,40 32 
100° 0,985 1,477 0,455 0,0117 1,4] iso} 
OS 0,940 1,410 0,606 0,0166 1,46 1,36 
1208 0,866 1,299 0,0254 1,58 1,44 
130° 0,766 1,149 0,0441 2,09 1,66 
140° 0,643 0,964 0,0917 


1 Siehe FuBnote 2 von S. 25. 
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Um jetzt die Verainderung der Lage des Ablésungspunktes mit der Drehgeschwindigkeit 


abschatzen zu koénnen, gehen wir von Gleichung (32) aus, die sich fiir den vorliegenden Fall der 
Kugel folgendermaBen schreibt: 


Ox Pe sinl—” PA lfn 
Rn (4) B( cos p4 (49) 
mit 1 A 
= cos P4 ltn Ux 
B Eee : 4 V. 2 R x (pa 0) 
—n a m m 
sin ao|+t 9 ta 


Die Kurven #,/R, (~4) nach (49) haben wir fiir den turbulenten Fall mit den Kurven 


Ox] Rm (p) nach (48) zum Schnitt gebracht, vel. Abb. 4, und daraus die Lagen der Ablisungs- 
punkte ermittelt. Das Ergebnis unserer Un- 


ae 1] 
tersuchung haben wir in Tabelle 2 mitgeteilt. - 
Mitangegeben sind auch die laminaren Werte Messing nach Luthander Ee LS UNA 
von H. Schlichting. Abb. 5 zeigt, wie in beiden = aa —~ 


Fallen laminarer und turbulenter Strémung 130" 
der Ablésungspunkt mit wachsender Dreh- 
geschwindigkeit nach vorn wandert, ein Ergeb- 
nis, das durch die Messungen von Luthander 


und Rydberg + bestatigt wird. 


turbulent 


120°|— 


Tabelle 2. Die Lage des Ablésungspunktes y4 bei 
laminarer und turbulenter Strémung fiir die angebla- 
 sene rotierende Kugel; Re = U,, Rm/v, Vm = Rm. 


{NS 


110° | t 
ie. farstnae ache turbulent |. _____|//aminar(nach 1. Schlichting) 
a. H. Schlichting Re = 10° 107 as 
0 LD © 126,5° 132502 ba che 
0,25 108,0° 126522 USGS 
0,50 HOP oe 125,4° 130,8° | | 
ae eee eee 
195 ’ 120,3° 126,1° Abb. 5. Die i ee paneer ane fiir die ange- 
1,50 (RES ARS Re =Uw Ry /v, m=OR,- 


Den Drehmomentenbeiwert bei voll-laminarer bzw. voll-turbulenter Strémung errechnen 
wir nach (37) mit x, = 0 und Kf = 0 und unter Beachtung von (21) zu 


aa Re Ifn 
U.. Rm ny alles || : 
CM ey ee = 2m 0( ==") a va[z | aetna (50) 
Cae e v 
“9 Rin 


In Abb. 6 und 7 haben wir tiber dem Winkel y, die Werte / Re a cy fiir den laminaren Fall 


und Ret cy fir den turbulenten Fall aufgetragen. Dabei ist fiir die Reynoldssche Zahl 


ire = U7, R,|v gesetzt worden. Es zeigt sich in beiden Fallen, da sich der Momentenbeiwert 
nur noch wenig andert, ob man bis zum Abliésungspunkt oder bis zum Kérperheck die Inte- 
gration erstreckt. Im laminaren Fall, Abb. 6, haben wir noch die Werte von H. Schlichting 
eingetragen, die mit unseren Werten recht gut iibereinstimmen®, Fiir p, =z 1aBt sich die 


Integration nach (50) geschlossen ausfithren und liefert 
1 


: (ey " its 
ore 3 Vx : 
Re't Co Cy = 276-6 Ag 1/6523 a z (51) 
erm 


1 Siehe FuBnote 8 von 5S. 21. ; Noh ; ; 
2 Die Ubereinstimmung wiirde wieder noch besser sein, wenn wir die Unterschiede in den Ansatzen 
fiir die Widerstandsbeiwerte (vgl. FuBnote 3 auf S. 29) beriicksichtigen wiirden. 
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Es bedeutet [’ die Eulersche Gammafunktion. Aus (51) ergibt sich fir den laminaren Fall 


9,76 U : 
= 2 laminar (51a) 
CM /Re We, ( ) 
und entsprechend fiir den turbulenten Fall 
_ 0,299 Us ars 
ty = i = Vn (turbulent). (51b) 


' Re=107, = 
| 7U. 
| /Re oP ow ta a 
ed ohhh 0 ae der Ablosun 
em ya | 5 Ven a5 9 
pool 
80 ie ele + G2 
a Bereich der Ablosung 
60 = ea 
0S gf 
Yeo 
60 
0 ie 
90° SO EE BO 6 16 0 120° 149° 6° 
Abb. 6. Der Drehmomentenbeiwert der Kugel bei laminarer Abb. 7. Der Drehmomentenbeiwert der Kugel bei turbulenter 
Strémung (g, = Winkel, bis zu dem die Integration ausgefiihrt Strémung (py, = Winkel, bis zu dem die Integration ausgefiihrt 
worden ist); worden ist); 
Q@ 72 pd 
Cy= M/2 V2 RS ; Re =UxRy/vi Vin = oOR,- Cy = M/S Vn RSs Re =U Ry [v3 Vi, = ORm- 


c) Rotierende Rotationsellipsoide. Fir die langsangestrémten sich 
drehenden Rotationsellipsoide wollen wir im folgenden die Drehmomentenbeiwerte nach (37) 
fiir die Falle voll-laminarer und voll-turbulenter Strémung berechnen. 

Fiir die Geschwindigkeitsverteilung eines in Richtung der groBen Achse L mit der Ge- 
schwindigkeit U,, angestrémten Rotationsellipsoid gilt in bekannter Weise 


Wy Um sin @ 


SSS et ee EE CY 52 
Ux Ux //sin® p+k? cos? p ee 

Hierbei gilt fiir die Kérperkontur die Parameterdarstellung 
X=Zcosy und R=Rysing, (53) 


wobei X von der Ellipsenmitte aus gemessen wird. U,, ist die maximale Kérpergeschwindig- 
keit. Sie ist eine Funktion des Achsenverhdltnisses 


Qin D ; 
p=- (54) 
und berechnet sich zu 
On 2, 
in ey (55) 


Zahlenwerte fiir aq sind von A. F. Zahm 1 angegeben worden. Fiir die Anderung der Lauflange x 
mit dem Polarwinkel ¢ gilt noch 


i = 4 Yin? p+ F cos® (56) 


1 A. F. Zahm, Flow and force equations for a body revolving in a fluid, NACA Report Nr. 323 (1928). 
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Die angegebenen Beziehungen fiihren wir in Gleichung (37) ein, in der wir x, = 0, Ki = 0 und 
l = L setzen. Die Formel, die wir finden, ist dadurch besonders einfach, weil sich die Wurzel 
im Nenner von (52) gegen die Wurzel von (56) heraushebt. Es folgt: 


1 
7 


M tet am 
tu = PAG Aa at | aut d : 57 


Wie wir bereits beim Beispiel der Kugel sahen, lat sich das Integral geschlossen lésen. Wir 
finden, wenn wir fiir cynoch die Beziehung (21) einsetzen, 
1 


me “J T 5+ 3n |i+n 
Ret"? eg =2n0\~ : 58 
eee oti SS oh A (58) 
dos 
J 
Fiir die Reynoldssche Zahl gilt jetzt Re = U,,L/y statt Re = U,, R,,/y wie beim Beispiel der 
Kugel. Zahlenmafig ergibt sich fiir den laminaren Fall 
Ane ae p= (laminar) (58a) 
und fiir den turbulenten Fall 
0.117 fae (Oe 
i 7 
Re wo) Vin 
Die Momentenbeiwerte sind auBer von der Reynoldsschen Zah] Re und vom Geschwindig- 
keitsverhaltnis U./Vin (Vm =@Rmn) noch vom Achsenverhiltnis k abhangig, das die GriBe 
U,,/ U,, nach (55) bestimmt. Fiir einige Werte : 
von k haben wir U,,/U,, berechnet und dar-  Tabelle 3. Die Drehmomentenbeiwerte bei laminarer 
aus die Momentenbeiwerte ermittelt. Das und turbulenter Strémung fiir DESL rotierende 
Ergebnis ist in Tabelle 3 angegeben. Da- Roteousellipseiden: 


(turbulent). (58b) 


co 


nach wachst der Drehmomentenbeiwert im laminar turbulent 
turbulenten Fall schneller mit dem Achsen- ipa iad fo a 

o : ; : : ‘ k pee, VRe C Rear 
verhaltnis des Rotationsellipsoides als im Ses Ue °M Uco °M 


laminaren Fall. In Abb. 8 haben wir den 


Drehmomentenbeiwert in Abhangigkeit von 0 1,000 | ae eee 
der Reynoldsschen Zahl fiir die beiden Falle ee een | Bee 0187 
laminarer und turbulenter Strémung fiir die 0,75 1.354 6,55 0,151 
verschiedenen Rotationsellipsoide darge- 1,0 1,500 6,91 0,165 
stellt. 1 Siehe Fufnote 1 von S. 34. 
| 
k=0 tI vi 
1 : < i 7 
tbe 
1 turbulent 
Th 
aed 2 y 6 8 we 2 ¢. 6 8 0 2 eS a Be a 


Abb. 8. Die Drehmomentenbeiwerte fir Rotationsellipsoide bei laminarer und turbulenter Strémung; 
2 2 
Cy = M/ 2 Vn RmL; Re = Ux L/; Vin = OR, - 
Achsenverhiltnis: k = D/L ; k = 0; Faden; k= 1: Kugel 


34. Truckenbrodt: Ein Quadraturverfahren zur Berechnung der Reibungsschicht usw. _Ingenieur-Archiv 
inj AER MMRDA TREES SeeNAE Sie Pine 2 Es 


d) Rotierender Halbkérper. Als letztes Beispiel wollen wir jetzt einen rota- 
tionssymmetrischen Halbkérper behandeln, der durch Uberlagerung einer rdumlichen Quell- 
strémung mit einer Translationsstrémung entsteht. 

Legen wir den Koordinatenanfangspunkt in den Staupunkt, dann gilt fiir die geometri- 
schen Daten des Kérpers die folgende Parameterdarstellung : 


ee ee ee (59) 


es | ee ce Ee e3) 60 
ere pees 9 ||) rie Sa (60) 


Es sind F und E die unvollstandigen elliptischen Integrale erster und zweiter Gattung fiir den 
Modul « = 60°?. 
Bezeichnen wir mit L die Abszisse auf der Achse des Kérpers vom Standpunkt aus ge- 


messen, dann ist 
L I cos Pp 
ETE ( A eo) 


cos — 
2 


Es bedeutet R,, den Kérperhalbmesser im Unendlichen L—» oo, Die Quelle befindet sich im 
Punkt L => R,,, von wo aus auch die Winkel m gemessen werden. 


Die Geschwindigkeitsverteilung ergibt sich zu 


U hee Oy oc 
Teg Pig (62) 
Die maximale Geschwindigkeit betragt 
ai 116 63 
Un. 3 ( 
und befindet sich an der Stelle cos p = — 1/3, d. h. L = 0,83 R,,. 


Die geometrischen Daten sowie die Geschwindigkeitsverteilung sind in Tabelle 4 mitgeteilt. 
Die Momentenbeiwerte wurden fiir die Falle voll-laminarer und voll-turbulenter Strémung 
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Abb. 9. Die Drehmomentenbeiwerte fiir den rotierenden Halbkérper bei laminarer Strémung; 
Cy = M/S VRS Re=Ux Ry, |; Vin = OR, « 
2 Der Fallk =1 (rotierende Kugel) wurde bereits als zweites Beispiel behandelt. Die Unterschiede 
Re Fs ER fiir a prep mone pte DBM ert (laminar 6,91 bzw. 9,76; turbulent 0,165 bzw. 0 299) 

ruhren her von den verschiedenen BezugsgréBen, Gl. (50) bzw. Gl. (5 ad pnt 
TEE HOP ELL ches ese (50) bzw (57), und der verschiedenen Definition 

ce 2 de g/2 
2 Es bedeutet F(”, P 4/3 
(3 3) ae wud Hea a [ /i—Ssim oa0. 
0 0 
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nach (37) mit x, = 0, Kt = 0,1 = R,, und x, = L berechnet. 


graphischen Darstellungen. 
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Abb. 10. Die Drehmomentenbeiwerte fiir den rotierenden Halbkérper bei turbulenter Strémung; 
Q 2 3 
cM =m/< Vin Rms Re =UcoRm/ys Vim = ORm 
Dabei gelten die folgenden Beziehungen: 
M (Wes, 1 
‘y= 3 Gs See, Wp = Oia. 


8. Zusammenfassung. Als 
Abschlu8B eines umfangrei- 
chen Programmes zur Berech- 
nung der Strémung an ange- 
blasenen rotierenden Rotati- 
onskérpern wird der bisher 
noch fehlende Fall turbulenter 
Strémung behandelt. 

Fir die Berechnung der 
Impulsverlustdicken sowie 
des Drehmomentes konnten 
einfache Quadraturformeln 
erhalten werden, die fiir den 
laminaren Fall ebenfalls eine 
recht brauchbare Naherung 
darstellen. In den Formeln 
kommen nur die Verteilung 
des Kérperhalbmessers und 
die Geschwindigkeitsvertei- 
lung um den Kérper vor. Es 
wird eine Abschatzung zur Be- 


+ View Ban 


Die Abb. 9 und 10 zeigen die 


Tabelle 4. Die geometrischen Daten und die Geschwindigkeitsverteilung 


Oe 
10° 
20° 
30° 
40° 
50° 
60° 
OE 
80° 
90° 

100° 
L105 
120° 
130° 
140° 
150° 
160° 
TiO 


des rotationssymmetrischen Halbkérpers. 


x L R U 
| Rn Rn Rn ce 
A 
0 0 0 0 

0,087 0,006 0,087 0,174 
0,176 0,023 0,174 0,343 
0,266 0,052 0,259 0,504 
0,358 0,092 0,342 0,653 
0,455 0,145 | 0,423 0,787 
0,557 0,211 0,500 0,901 
0.665 0,291 0,574 0,996 
0,783 0,387 0,643 1,068 
0,914 0,500 0,707 1,118 
1,061 0.635 0,766 1,146 
L233 0,798 0,819 1,155 
1,440 1,000 0,866 1,146 
1,703 1,260 0,906 ES) 
2,065 1,620 0,940 1,092 
2,618 Dales 0,966 1,059 
3,651 3,206 0,985 1,028 
6,595 6,149 0,996 1,007 
foe) oo 1,000 1,000 


180° 


stimmung der Lage des Ablésungspunktes infolge der Drehbewegung des Kérpers angegeben. 

Das Verfahren wurde erprobt an den Beispielen der rotierenden Scheibe, der Kugel, der 
Rotationsellipsoide und des rdumlichen Halbkérpers. In allen Fallen wurde der Momenten- 
beiwert bei laminarer und turbulenter Strémung berechnet. Allgemein kann man zeigen, daB 
sich das Drehmoment bei konstant gehaltener Anstrémungsgeschwindigkeit linear mit der 
Umfangsgeschwindigkeit andert. 


(Eingegangen am 23. Mai 1953.) 


Anschrift des Verfassers: Priv.-Doz. Dr.-Ing. E. Truckenbrodt, Braunschweig, Kérnerstrafe 12. 
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Zur Erzielung optimaler Wirkung bei Pendel-Schwingungstilgern. 
Von A. Slibar und K, Desoyer. 


1. Einleitung. Bekanntlich liefert die linearisierte Behandlung *) des an einer rotierenden 
Masse angelenkten Pendel-Schwingungstilgers die als Abstimmbedingung bekannte Relation 


@ =F of, (1) 
welche zur Erzielung giinstigster Tilgerwirkung einzuhalten ist. Darin bedeuten 92 die Kreis- 
fre quenz der erregenden Momente, L die Entfernung des Anlenkpunktes des Pendels vom Dreh- 
punkt der rotierenden Masse, w, die mittlere Winkelgeschwindigkeit der umlaufenden Massen, 
l die Lange des als mathematisches Pendel aufgefaBten Tilgers. Die Gleichung (1) nimmt die 
einfache Form 


2 == yp (1’) 


an, wenn die Kreisfrequenz Q der erregenden Momente als ein y-faches der mittleren Winkel- 
geschwindigkeit @, der rotierenden Masse eingefiithrt wird. Die linearisierte Rechnung liefert 
keine Beziehungen zur Bestimmung von optimalen Werten fiir die Pendelmasse, sondern die 
Aussage, diese méglichst groB zu machen. 

In einer vorangegangenen Arbeit * konnten die Verfasser zeigen, dafs sich bei nichtlineari- 
sierter Behandlung des Tilgers als mathematisches Pendel cine Bestimmungsgleichung fiir die 
optimale Pendelmasse auffinden la8t. Die Rechnung lieferte neben der Abstimmbedingung (1’) 
als Bestimmungsgleichung fiir optimale Tilgungsverhaltnisse die Beziehung 


mE\ . 394 _ (2) 
( 7 io 2(p? + 1)2° 


Darin bedeuten m die Masse des Pendels, J das Tragheitsmoment der rotierenden Masse, an die 
das Pendel angelenkt ist, L und J haben die bereits oben angegebene Bedeutung. Die Rechnung 
ergab weiters, da bei gewahltem L und J die doppelten Werte der aus (2) folgenden optimalen 
Massen zu unerwiinschten Resonanzzustanden fithren. 


Ist der Anlenkpunkt des Tilgers aus konstruktiven Griinden festgelegt, so sind zur Erzielung 
giinstigster Tilgungswirkung des Pendels, wenn dieses als mathematisches aufgefaBt wird, nur 
die zwei GréBen J und m frei verfiigbar. Legt man der Betrachtung jedoch ein physikalisches 
Pendel zugrunde, so ist der Tilger durch Festlegung von drei voneinander unabhangigen 
Gré8en, namlich der Masse m, des Abstandes s des Pendelschwerpunktes vom Anlenkpunkt 
und des Tragheitshalbmessers i, des Pendels um seinen Schwerpunkt, zu bestimmen. Die Be- 
trachtung von ausgefiihrten Tilgerformen ° zeigt, daB die rechnerische Behandlung von Tilgern 
als mathematische Pendel (mit Ausnahme eines Sonderfalles, vgl. Abschnitt 5 dieser Arbeit) in 
den meisten Fallen nur als eine grobe Naherung fiir das tatsachliche Verhalten des angelenkten 
Pendelkérpers angesehen werden kann. Zweifellos fiihrt aber gerade die Méglichkeit, drei von- 
einander unabhangige Parameter (m, s, i,) méglichst giinstig zu kombinieren, zu einer wesent- 
lich gréBeren Freiziigigkeit in der Konstruktion. Aus diesem Grunde soll im folgenden der als 
physikalisches Pendel aufgefaBte Tilger, auch unter Einbeziechung gré®erer Ausschlagwinkel, 
behandelt werden. 


Die vorliegende Berechnung der Tilgungswirkung eines Pendels an einer einzigen Dreh- 
masse ist vornehmlich fiir Sternmotoren anwendbar. Die Einwirkung von mehreren physika- 


1 0. Kraemer, Z. VDI 82 (1938), S. 1297, Z. VDI 83 (1939), S. 901, Motortechn, Z. 1 (1939), S..3: 
W. Schick, Ing.-Arch. 10 (1939), S. 303. 
A, Kimmel u. I. Lutzweiler, Ing.-Arch. 12 (1941), S. 100. 
R. Lambrich, Das Fliehkraftresonanzpendel, Jb. dtsch. Luftfahrtforsch. 1940, S. 1170. 
H. Kammerer, Die Eigenfrequenzkurven bei Drehschwingungssystemen mit Fliehkraftpendeln, 
Jb. dtsch, Luftfahrtf. 1940, S. 1180. 
2 K. Desoyer u. A. Slibar, Ing.-Arch. 22 (1953), S. 208. 
° z. B. H. Schrén, Die Dynamik der Verbrennungskraftmaschine. Wien, 1947, 
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lischen Pendeln auf eine einfach oder mehrfach zusammenhangende elastische Kette 


- : 
5 ‘ soll in 
einer nachfolgenden Arbeit untersucht werden. 


2. Die Bewegungsgleichungen. Der Rechnung liegt eine 
Anordnung gema8B Abb. 1 zugrunde. 

An einer Schwungmasse vom Tragheitsmoment J sei im 
Abstand L von der Drehachse ein physikalisches Pendel 
der Masse m angelenkt. Der Abstand zwischen dessen 
Schwerpunkt S und dem Anlenkpunkt A sei s. Der Trag- 
heitshalbmesser des Pendelkérpers beziiglich seines Schwer- 
punktes sei t,. Die Bedeutung der weiteren Bezeichnungen 


geht aus Abb. 1 hervor. 


Auf die Schwungmasse wirke ein auBeres Drehmoment 


M = M, -|- AM sin y» Q> (3) Abb. 1. Schema eines einfach angelenkten 
Pendel-Schwingungstilgers. 
von dem der konstante Betrag M, zur Aufrechterhaltung 


der konstanten mittleren Winkelgeschwindigkeit w, der Anordnung gegen die Arbeitswider- 
stande abgespaltet werde. 


In die Lagrangeschen Gleichungen 
d (oT oT 
dt (Fae) ~~ Ogn = 0 (4) 
fiihren wir als verallgemeinerte Koordinaten die beiden Winkel 
Cae, Ore? 


ein, die die Lage des Systems in jedem Augenblick kennzeichnen. Setzt man zur Bestimmung 
der verallgemeinerten Krafte Q, die Gleichheit der differentiellen Arbeiten 


Q,dp + Q,d8 =AMsinyvqp: dp (5) 
an, so ergibt sich 
0, =AMsinvy, QO, = 0. (5’) 
Die kinetische Energie T des gesamten Systems ist 
Tate 1 Mog 4 Ep + 8). (6) 


Aus den Koordinaten xs und ys des Pendelschwerpunktes in dem nicht mitrotierenden Koordi- 
natensystem der Abb. 1 


xs = L cos p + scos (y + #),") 
ys = Lsing +ssin(p + 9), | 
erhalt man das Quadrat der Absolutgeschwindigkeit vs des Pendelschwerpunktes: 


=a +78 =[(P +2 +2Lscos dP? +s Hh +292%~Pb+12LshGcosh. (8) 


(7) 


Die kinetische Energie T des gesamten Systems wird damit 


: > - ° mie. ° ; 
T HE 4 ML + 8421 scos 0] Gp? + F 28? + mse b+ mL sp Geos d+ = (p +9). (6’) 


Durch Einsetzen von T und dessen Ableitungen sowie von Q, bzw. Qy = 0 in (4) erhalt man 
die beiden Bewegungsgleichungen in der Form 


(J+ mL? + 2mL scos 8) p + m(s* +33) (G +0) —2mLsp dsin 9 — | 
—_mLs#sin 0 +mLs cos 8 =AMsinvg, 
(s? + 18) (p +9) + Ls (~ cos? + sin 7) = 0. | 


1 R.Grammel, K. Klotter und K. v. Sanden, Ing.-Arch, 7 (1936), 5. 439. 


W. Schick, Ing.-Arch. 10 (1939), S. 303. 
A. Kimmel und I. Lutzweiler, Ing.-Arch. 12 (1941), S. 100. 


(9) 
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3. Auswertung der Bewegungsgleichungen. Zerlegt man die Winkelgeschwindigkeit wm = 
in den konstanten Anteil @, und die Stérung Jw 


WO = + 4o 


dann ist wegen @ = 0 die GroBe @ durch Aw zu ersetzen. Fiihrt man die obige Zerlegung in 
die Gleichungen (9) ein und vernachlassigt wie iiblich Jw gegen w , so nimmt das System (9) 
die folgende Form an: 


(J + mL? 42m L's cos 8) Ady + m (i§ + 5%) (4a + 8) —2mL s wy O sin 6 — 


—mLs#sin 0 +mLs bcos # = AMsinv at, (9’) 
2 i ae 
: 8 (Ab + #8) +Aq@cos & + af sin ®? = 0. 
§ 
Fiihrt man die dimensionslosen GréBen 
s® + ig ; > mL? 4M _ A kde 10 
poe Ere ee sje Jo ’ Wo ( ) 


. . o y 
ein und kennzeichnet man die Ableitungen nach 7 durch einen Strich, so erhalt man aus (9’) 


+06 + Oo0i* +200 cos 8) _ 2608" sin? — Oo #2 sin B 


+ @Oa(A*+ cos #) 8’ = K sin yt, (11) 
Ae’ Ae OY + sin 
Wo — A* + cos O 


Beriicksichtigt man, da8 #, im Bogenma8 gemessen, praktisch stets kleiner als 1 bleibt und 
verwendet man demgema8 die Entwicklungen 


é gs 9 1 1 il Ba pe 
snPazd—Z, coodal—-z, Eee as 


unter Streichung der Glieder héherer als dritter Ordnung in # bzw. der Produkte von # und 
dessen Ableitungen, so geht die zweite Gleichung (11) in die Form iiber 
Aa’ 
ay 


= Fh" + kh + c¥ B34 ch" GB, (12) 


wobei zur Abkiirzung gesetzt wurde 
re z 1 . ay: ‘ a 


ere. oe fig 


1? mr + A* 2 Can (12a) 


64a) 2? Sot asi 
Die erste Gleichung (11) geht bei gleicher Entwicklung und Verwendung von (12) iiber in 
ath’ + akh + aF 08 + af" G + af 8? O + ak 0'O = K sin yt (13) 


mit den Abkirzungen 


1=[L+ 6+ Oo(2+4*)] f+O0(1+A*), af =[1+ O+ Oo(2+A*)] & — 00(2 +3) 
*=[1+ 0+ Oo(2+2*)] cf, a*— — Oo, (ie 
*=([1+ 0+ Oo(2+2*)] c¥ — Oo cf, a*= —2@0o. 


Die Gleichungen (12) und (13) sind formal gleich wie die fiir das mathematische Pendel giiltigen 
Gleichungen (12) und (13) der oben angefithrten Arbeit, nur die Bedeutung der Koeffizienten 
ist eine andere. Man iiberzeugt sich leicht, daB die Koeffizienten c, bis c} bzw. a; bis a; in die 
entsprechenden Grofen c, bis c, bzw. a, bis a, der dortigen Beziehungen iibergehen, wenn man 
auf die urspriinglichen Bedeutungen der eingefithrten Abkiirzungen zuriickgeht, die Lange s 
durch / ersetzt und weiters den Tragheitsradius igs gleich Null setzt. 


Wegen der formalen Ubereinstimmung der Gleichungen (12) und (13) mit den gleich be- 
zifferten Beziehungen des mathematischen Pendels kann die weitere Behandlung dieser beiden 
nichtlinearen Differentialgleichungen vollkommen analog durchgefiihrt werden. Der Ausdruck 
fiir die bezogene Winkelgeschwindigkeitsschwankung der Schwungmasse ergibt sich also ge- 


1 Siehe FuBnote 2 von S. 36. 
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maB Gleichung (19) der angefithrten Arbeit, wenn man darin die GréGen c, bis c, baw. a, bis a, 
durch die entsprechenden mit Stern versehenen GréfBen ersetzt : 


Aw | 1 fe —v? cf ss o,.#\ L3(a¥ —v? af)+4+ »? af] 3(c# — v2 cf) 
% ~  » \af— 0 af ew St) 4 (a*¥ —v? af) ee Gran ay) K+leosoa-+| 
[ex — (2 »)? cf] at 2 
+ Fey — Gy! oF] (af — oF ape S008 27 — te 
1 { [c¥ — (3 v)? c¥] [a* — v? (af + a*)] (cX —v? cf) 4 
By ATek — Bo)! oF] (oF vay? Tag —vPapyyy COS STH 


Da der Parameter K die gleiche Bedeutung wie in der oben angegebenen Arbeit hat, gelten fiir 
dessen Gré®enabschatzung die gleichen Gesichtspunkte. K ergibt sich an Hand von ausge- 
fiihrten Sternmotoren in der GréBenordnung von 10-5 bis 10-3 je nach Tangentialdruckdia- 
gramm, Kurbelwellenbauart und Drehzahl der Maschine. 

Die Gleichung (14) zeigt, daB sich das Glied erster Ordnung in K zu Null machen 1aBt. 
Setzt man fiir c¥, c¥, a¥, und a¥ gemaB (12a) und (13a) ein, so ergibt sich namlich 


cy — vce 1 — py? j* (15) 
at — vat [1+ O+ Oo2+a4)] —1*) + Oo(l +A*pP° 


Da der Nenner der rechten Seite von (15) nicht unendlich werden kann, mu8 fiir Verschwinden 
dieses Ausdrucks der Zahler zu Null gemacht werden. Die erste Abstimmbedingung fiir den als 
physikalisches Pendel gerechneten Tilger lautet somit 


T= 2A" = 0 (16) 

bzw. mit der Bedeutung von A* gemafB (10) 
Side 1 : 
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Abb. 2. Schichtenliniendiagramm fiir die Amplitude der verbleibenden bezogenen Winkelgeschwindigkeitsschwankung. 
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Nach Erfiillung dieser ersten Abstimmbedingung bleibt fir die bezogene Winkelgeschwindig- 
keitsschwankung der Ausdruck 


Aw y 3 Riggs Dero 

Se ee CORY E. 4 2 COR aT 

Tam aes h: 2 @o(v® + 1) jaca + 0} 60) an 
EEoeus K3 cos 3yt +--+ *<. 


8 v(v2 + 1)8 62 0? SF (9vt+208 +1) — 811 +8)| 


Das ausschlaggebende Glied in (17) stellt der mit K® multiplizierte Term dar. Es ist also danach 
zu streben, den Absolutbetrag von 
3 
T= = (18) 


4 y4 2? il 
2 Oa(v? + 1)2 [aa | ©)—Oc v > ae 


méglichst klein zu halten. Dies ist durch entsprechende Wahl der Parameter O = mL?/J und 
o = s/Linnerhalb der konstruktiven Schranken dieser beiden GréBen zu erreichen. 

Um zu einer fiir alle Werte von v verwendbaren graphischen Darstellung der GréfBe I’ als 
Funktion der beiden Parameter © und o zu gelangen, formen wir unter Verwendung der Ab- 
kiirzung 


Views cee oe ahs (18’) 
den Ausdruck (18) um in 
2412 5 3 
pe Ns O(c V) [311+ 6) — O(6P)] (18a) 


Die Beziehung (18a) wurde im Schichtendiagramm Abb. 2 und im Relief Abb. 3 dargestellt. 
Die an den Schichtenlinien eingetragenen Zahlen bedeuten jeweils den Wert der GréBe 

(»? + 1)2 

—F ie 


Sieht man von den Gliedern mit dritter und héherer Ordnung in K in (17) ab, so bleibt 
Aw 


Mo 


== 1 K* cos 290, t. (19) 


Unter dieser Voraussetzung stellen die an den Schichtenlinien der Abb. 2 und 3 eingetragenen 
Werte die Amplitude der GriBe (? +1? do 


VK?" dar. 


Der Schnitt des rdumlichen Diagrammes Abb. 3 mit den Ebenen 


c= es = Kosh. 
v 

liefert die bereits bekannten ebenen Diagramme ! fiir die Abhangigkeit der Amplitude von 
Aw/m,K? vom Parameter @ fiir das mathematische Pendel. 

Die Diagramme Abb. 2 und 3 zeigen, daB absolute Minima im Endlichen nicht existieren. 
Zu beachten ist, daf die Ordinaten in den Schaubildern unter der Voraussetzung der Einhaltung 
der ersten Abstimmbedingung (16) des physikalischen Pendels errechnet wurden. Die Be- 
ziehung (18a) zeigt, daB ein als Resonanz anzusprechender Effekt fiir den Fall 


F(0, 0) = Oo (1 + 0) — Oo VI] =0 (20) 
existiert. Das Gebiet F(0,0)> 0 ist im Schichtenliniendiagramm Abb. 2 mit A, das Gebiet 
F(0,c0) <0 mit B bezeichnet. Der Bereich A wird eingeschlossen von der 0 —, og V-Achse 
und der Resonanzlinie mit der Bezifferung -- oo, wahrend der Bereich B auf der anderen Seite 
der Resonanzlinie liegt. 


Geht man in (20) auf die urspriinglichen Bedeutungen der Parameter @ und o zuriick. so 
erhalt man 


F(m, L, ) = 7 BJ Lm L?)—mLsV). (21) 


1. Siehe FuBnote 2 von S, 36, 
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Abb. 3. Relief fiir die Amplitude der verbleibenden bezogenen Winkelgeschwindigkeitsschwankung. 


Die partiellen Ableitungen von F nach m, L und sliefern, gleich Null gesetzt, die im Schichten- 
liniendiagramm strichliert eingezeichneten Kurven, auf deren Bedeutung im folgenden einge- 
gangen wird. 


4, Zur Berechnung und Konstruktion. Je nach den konstruktiven Forderungen werden 
verschiedene Wege bei der Durchrechnung eines Tilgers zu beschreiten sein. In den meisten 
Fallen werden die GréBen J, 4M und @, als vorgegeben vorauszusetzen sein. Zulassig ist die 
freie Wahl von dreien der folgenden fiinf GréBen: (4m) zulassig, m, s, L und is. Uber die zu 
tilgende Ordnung vy muB — nach ausgefithrter Schwingungsrechnung fiir die tilgerlose Anord- 
nung — ein Entscheid getroffen werden. 

Je nach Auswahl der drei wahlbaren GréBen sind folgende Méglichkeiten zu unterscheiden: 

a) Sei z. B. die Pendelmasse m und der zulassige Hichstwert (4m), der Schwankung ge- 
geben. Dann ist vorerst die GroBe @ 


I’ = Amplitude von - er 
zu bestimmen. Entsprechend der zu tilgenden Ordnung y sind die zugehérigen Schichtenlinien 
im Diagramm Abb. 2 aufzusuchen. 

Jeder einzelne Punkt der beiden, dem gegebenen Wert von(4w),1 entsprechenden Schichten- 
linien liefert zum angenommenen Wert von m eine mégliche Kombination von s und L. Im 
allgemeinen wird sich die Bestimmung von s und L aus einem Punkt der Schichtenlinie im Be- 
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reich B aus konstruktiven Griinden nicht durchfiihren lassen; zur Konstruktion werden daher 
die Schichtenlinien im Bereich A zu verwenden sein. Nun kann noch entweder o oder L selbst 
gewihlt werden. Im ersteren Fall bestimmen sich auf der dem gewahlten o entsprechenden 
Parallelen zur @-Achse zwei Schnittpunkte mit der betreffenden Schichtenlinie. Damit liegen 
die beiden zugehorigen Werte fiir O fest. Aus den Werten fiir O folgen die Werte fiir s und L 
und ferner aus der ersten Abstimmbedingung (16’) die zugehérigen Werte des Tragheitshalb- 
messers ts. 

Wird im zweiten Fall L selbst aus der Konstruktion entnommen, so ist bei gewahltem m 
ein bestimmter Wert fiir 9 = mL?2/J gegeben. Auf der entsprechenden Parallelen zur o V-Achse 
ergeben sich zwei Schnittpunkte mit der gewahlten Schichtenlinie im BereichA und damit zwei 
mégliche Werte 0,5. Mit den so erhaltenen Werten 0,,. bestimmen sich beim vorgewahlten L 
die beiden Méglichkeiten fiir s. Die erste Abstimmbedingung (16’) liefert wieder das zugehorige 
Tgaia: 
Pa im zweiten der besprochenen Falle an Stelle einer vorgeschriebenen héchstzulassigen 
Winkel geschwindigkeitsschwankung die Forderung nach minimaler Schwankung zu erfiillen, 
so ist der gleiche Weg einzuschlagen und lediglich an Stelle der Schichtenlinie die strichliert 
eingezeichnete Kurve I aufzusuchen. Diese wurde durch partielle Differentiation von (21) 
nach s gewonnen. 


b) AuBer den vorgegebenen GréBen J, AM und a, seien (4),_1, L und is gewahlt; aus der 
ersten Abstimmbedingung (16’) folgen sofort zwei mégliche Werte fiir s: 


L EN2 : 
43 = 2 View) — id. (22) 


Dies zeigt, das bereits bei der Wahl von L und igs die Ungleichung 
L S27 ts (23) 


einzuhalten ist. Geht man mit den nun festliegenden Werten o,,. = $;,./L in das Schichten- 
liniendiagramm zu den Schnittpunkten mit der durch das zugelassene (4@),4 bestimmten 
Schichtenlinie, so ergeben sich die méglichen Werte fiir O, aus denen die zugehérigen optimalen 
Pendelmassen zu entnehmen sind. 

Wird im obigen Fall minimale Winkelgeschwindigkeitsschwankung gefordert, so ist an 
Stelle einer bestimmten Schichtenlinie die strichliert eingezeichnete Kurve II zu verwenden, 
welche aus (21) durch partielle Differentiation nach m gewonnen wurde. 


c) AuBer den vorgegebenen Gréfen J, 1M und wy seien (4@).1, m und s frei gewahlt; 
man geht mit einer Fehlannahme fiir 0 = s/Lin das Schichtenliniendiagramm ein und bestimmt 
mit den Schnittpunkten auf den dem vorgegebenen (A), entsprechenden Héhenlinien die 
zugehorigen Werte fiir @. Die sich daraus ergebenden L werden im allgemeinen der Beziehung 
OFehlannahme = $/L widersprechen. Es kann aber leicht ersehen werden, in welcher Richtung 
der Wert von o abzuandern ist. Nach wenigen Schritten wird sich die richtige GréRe fiir L 
ermittelnlassen. Kine zweite Méglichkeit der Bestimmung von L ware die rechnerische Ermitt- 
Jung aus (18), wenn man darin die GréBen O und o durch ihre urspriinglichen Bedeutungen 
ersetzt. Man erhalt 


=e Aw 
foe = Amplitude von any’ 
2m Le(t + 12 18(J 4 mL) — mg ee (oy K? ) 


ye 


Da diese Beziehung eine Gleichung dritten Grades in L darstellt, wird im allgemeinen der 
graphische Weg vorzuziehen sein. 

Mit dem gefundenen Wert von L und dem gewahlten s lat sich aus der ersten Abstimm- 
bedingung (16’) sofort is errechnen. 

Ist bei festem m und s eine minimale Winkelgeschwindigkeitsschwankung gefordert, so ist 
lediglich an Stelle einer bestimmten Schichtenlinie die in den Schaubildern strichliert ein- 
gezeichnete Kurve IJ zu verwenden, die durch partielle Differentiation von (21) nach L 
gewonnen wurde. 


d) Will man einen bereits in seinen Hauptabmessungen festgelegten Schwingungstilger, 
der die erste Abstimmbedingung (16’) erfillt, auf seine Wirksamkeit iiberpriifen, so ist einfach 
mit den bekannten Werten von @ und g in das Schichtenliniendiagramm einzugehen. 
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e) Sind konstruktive Schranken fiir m, L und s vorgegeben, so bestimmen diese ein einzu- 
haltendes Gebiet im Diagramm Abb. 2. In diesem Gebiet ist der Punkt, welcher der kleinsten 
Winkel geschwindigkeitsschwankung entspricht, aufzu- 
suchen. Er liefert die zugehérigen Kombinationen von 
m, Lund s. Aus der ersten Abstimmbedingung ist is zu 
errechnen. 


5. Der bifilar angelenkte Pendel-Tilger. Aus kon- 
struktiven Griinden werden Pendel-Tilger haufig nach 
dem in Abb. 4 dargestellten Schema ausgefiihrt. Bei 
dieser Art der Anlenkung ist die Relativbewegung der 
Pendelmasse gegeniiber der Hauptmasse eine reine 
Translation. Der Schwerpunkt S der Pendelmasse be- 
schreibt gegeniiber der rotierenden Hauptmasse eine 
Kreisbahn vom Halbmesser s um den Punkt A (Abb. 4). ae ee Fe ee ee 
Der Abstand OA werde mit L, das Tragheitsmoment aie ti BobWwingideidleccer ? mee ti 
der Hauptmasse mit J, und der Tragheitshalbmesser des 
Pendelkérpers um seinen Schwerpunkt S mit is, bezeichnet. Die Bedeutung der ubrigen 
Bezeichnungen geht aus Abb. 4 hervor. 


Die kinetische Energie T des gesamten Systems ist 


v2 


2 v s 
T = (J, + mis,) > Weer (24) 


Das Quadrat der Absolutgeschwindigkeit vg des Schwerpunktes des Pendelkérpers ergibt sich 
aus den Schwerpunktskoordinaten xs und ys in bezug auf das nicht mitrotierende Koordina- 


tensystem der Abb. 4: 


of = 28 +938 =[(L2 +82 +12Lscos dg? +22 4+22%¢9+2Lsp hcosd. (25) 
Durch Einsetzen in Gleichung (24) erhalt man 


“2 : : ; 
T= (J, + mi.) > + [1 +s? +2Lscos 3] 9? ++ vP+tHms~pd+mLsy dcosd. (24) 


Der Vergleich mit dem Ausdruck (6’) zeigt, daB der fiir das physikalische Pendel durchgefiihrte 
Rechnungsgang verwendbar ist, wenn man die dort eingefiihrte GréBe J durch den Ausdruck 
J, +m ig, und weiter den in der Beziehung (6’) verwendeten Trigheitshalbmesser is durch 
Null ersetzt. 

Die erste Abstimmbedingung (16’) des physikalischen Pendels geht damit iiber in die erste 
Abstimmbedingung 


eal (26) 


des mathematischen Pendels. Der Ausdruck nach (18), der in Hinblick auf die verbleibende 
Winkel geschwindigkeitsschwankung zum Minimum zu machen ist, geht dann tiber in 


P= sor7p a 1? © — 304] (27) 
mit 
eres et (28) 
Der Parameter K erhalt die Bedeutung 
AM 28’) 


(J; + m 13,) - 2° 


Abgesehen von der geanderten Bedeutung des Parameters 0 stimmt der Ausdruck J" mit dem 
Koeffizienten des mit K? multiplizierten Termes der Gleichung (20) der bereits zitierten Arbeit 


der Verfasser 1 iiberein. Die Abhangigkeit der Grébe I’ vom Parameter @ kann demnach aus 


1 Siehe Fufnote 2 von S. 36. 
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der Abb. 2 der angegebenen Arbeit entnommen werden. Die zweite Abstimmbedingung fiir 
bilifar angelenkte Tilger lautet analog (22) obiger Arbeit 
is Vins eee 

Oot = Ji + mopt 13, ~~ 2(y2 + 1)?” (29) 
die durch Differentiation des Ausdruckes J’ nach @ erhalten wird. Kine genaue Untersuchung 
der Abhangigkeit der GréBe I’ von den beiden prinzipiell voneinander unabhangigen Para- 
metern m und is, hat hier keine praktische Bedeutung, da erstens der Parameter is, fiir die 
Abstimmung nicht mit der gleichen Wirksamkeit eingeht wie bei einfach angelenkten physi- 
kalischen Pendeln und zweitens die im allgemeinen konstruktiv bedingte Form des Pendel- 
kérpers bei der bifilaren Ausfithrung nicht wesentlich variiert werden kann und daher mit 
bestimmter Pendelmasse m der zugehérige Wert is, im wesentlichen festgelegt ist. 


6. Zusammenfassung. Die Betrachtung ausgefiihrter Konstruktionen von Pendel-Schwin- 
gungstilgern zeigt, da® die rechnerische Behandlung von einfach angelenkten Pendelkérpern 
als mathematische Pendel nur grobe Naherungen liefern kann. In der vorliegenden Arbeit 
wird die Einwirkung und zweckentsprechende Dimensionierung eines an einer rotierenden 
Masse angelenkten Tilgers, der als physikalisches Pendel behandelt wird, unter Einbeziehung 
groBer Pendelausschlage untersucht. 

Bei Zugrundelegung eines mathematischen Pendels stehen an Verdnderlichen die Pendel- 
masse m, die Pendellange 1 und die Entfernung L des Anlenkpunktes des Pendels vom Wellen- 
mittel zur Verfiigung. Geht man jedoch vom physikalischen Pendel aus, so tritt der Tragheits- 
halbmesser is| des Pendelkérpers um seinen Schwerpunkt als weiterer Parameter hinzu. Die 
Méglichkeit, die vier voneinander unabhangigen Parameter m, L, s und is méglichst giinstig 
zu kombinieren, gibt eine wesentlich gréBere Freiziigigkeit in der Konstruktion, um geringste 
Winkel geschwindigkeitsschwankung der rotierenden Hauptmasse zu erzielen. Es zeigt sich, 
dafB die fiir das mathematische Pendel giiltige erste Abstimmbedingung (1’) auch fiir groBe 
Ausschlage in die Bedingung (16’) iibergeht. Die Amplitude der nach Erfiillung dieser Bedin- 
gung verbleibenden Winkelgeschwindigkeitsschwankung wird in einem Schichtenliniendia- 
gramm in Abhangigkeit von zwei dimensionslosen Parametern dargestellt. Das beigefiigte 
Relief (Abb. 3) zeigt, ebenso wie das Schichtenliniendiagramm Abb. 2, daB, wie bei der nicht- 
linearisierten Durchrechnung des mathematischen Pendels auch beim physikalischen Pendel 
teotz EKinhaltung der ersten Abstimmbedingung bei ungiinstiger Kombination der Parameter 
unerwiinschte Resonanz auch im konstruktiv erreichbaren Bereich auftreten kann. 

AnschlieBend werden die méglichen Wege zur praktischen Berechnung der KenngréfSen 
eines einfach angelenkten Pendel-Schwingungstilgers zusammengestellt. 

Die abschlieBend durchgefiihrte Behandlung des bifilar angelenkten Pendel-Tilgers ergibt 
weitgehende Ubereinstimmung mit den Ergebnissen der fiir das mathematische Pendel durch- 
gefiihrten nichtlinearisierten Rechnung. 

Die vorliegende Berechnung der Tilgungswirkung eines Pendels an einer rotierenden Haupt- 
masse ist vornehmlich fiir Sternmotore anwendbar, wobei der Ubersichtlichkeit halber vor- 
laufig vorausgesetzt wurde, das die Drehmasse keine Eigenfrequenzen besitze. Hine analoge 
Untersuchung iiber die Wirkung und zweckentsprechende Dimensionierung mehrerer an 
ein schwingungsfahiges Gebilde angelenkter Pendel-Tilger soll spater erfolgen. 


(Eingegangen am 26. Mai 1953.) 


Anschrift der Verfasser: Dozent Dr. A. Slibar und Dr. K. Desoyer, Wien IV., Karlspl. 13, Technische 
Hochschule. 
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Uber den Reziprozititssatz in der Dynamik der elastischen Korper. 
Von D. Graffi. 


Man betrachtet einen elastischen Kérper mit dem Volumen V und der Oberflache o. Nennen 
wir F dV und R do die Krafte, die auf das Kérperelement dV und auf das Oberflachenelement 
do wirken, und s die zugehérige Verschiebung eines beliebigen Punktes P des Kérpers. Ferner 
sei s’ die Verschiebung von P wenn statt F dV und Rdo die Krafte F’ dV und R’ do wirken. 
Die Krafte und die Verschiebungen mégen im allgemein variabel mit der Zeit sein. Ist 0 die 
Dichte in P, so hat man, wie Lord Rayleigh bemerkt hat!, zu jeder Zeit die Beziehung 


[Ss 2s —; Se si 5 mld OAS \ R’< 
| (Peds) av + | Rs'do=[(F—oSr)eav + [ Rida. (1) 
v 5 y a 


Die Gleichung (1) ist eine einfache Folgerung der Gleichungen der Dynamik der deformier- 
baren Kérper und des Hookeschen Gesetzes, das im Fall der isotropen Kérper, auf die wir uns 
beschranken, in folgender Form geschrieben werden kann: 


O47 =2Me 7; tAGj0, O=e, +ey +e (i,j =1,2, 3), (2) 
wo 0;; das Kroneckersche Symbol, o;; die Spannungen, ¢;; die Dehnungen, 4 und wu die Lamé- 
schen Konstanten sind. 

Die Gleichung (1) enthalt aber die im allgemeinen unbekannten Tragheitskrafte 0 ee ee 
Nur im statischen Fall verschwinden diese Krafte und (1) reduziert sich auf den bekannten 
Bettischen Reziprozitatssatz ?: 

Unser Ziel ist nun, eine Gleichung zu erhalten, die analog zu (1) ist, in der aber keine Trag- 


 heitskrafte auftreten. Der Einfachheit halber nehmen wir an, daB, fiir t = 0 die Verschiebun- 


gen s,s’ und ihre ersten Ableitungen nach t verschwinden. Dann wenden wir auf die Gleichun- 
gen der deformierbaren Kérper und auf das Hookesche Gesetz (2) die Laplace-Transformation 
an, das heift, wir multiplizieren jene Gleichungen mit e—P* (p komplexe Zahl mit negativem 
Realteil) und integrieren iiber t von Null bis Unendlich. Wir setzen voraus, da8 die so erhalte- 
nen Integrale konvergieren, und da bestimmte elementare Operationen erlaubt sind (wie 
Vertauschung der Ableitung nach den Raumkoordinaten mit der Integration nach t usw.), 
Voraussetzungen, die hier stets als giiltig angesehen werden diirfen. Es ist alsdann leicht, zu 
erkennen, daB die Gleichungen (2) und die anderen Gleichungen formal gleich bleiben, wenn 
nur s, Fusw., durch ihre Laplace-Transformierten {(s), &(2) usw. ersetzt werden. Daraus 
folgt, daB auch (1) gilt, wenn wir die dortigen GréBen durch ihre Laplace-Transformierten 
ersetzen. Erinnern wir uns daran, daB, unter unseren Voraussetzungen 


&(5r) =P 26) 
ist, so kommt . 
J &(F) Qs’) dV — p? f o &(S) 2G’) dV + f Q(R) BG’) do | 
ha V o 
= [ QF) 26) av — p* Je BR) LE) AV +f LR) LE) do. | 
Vv o 
Wendet man jetzt den Faltungssatz an, so erhalt man sofort 


ii dx [ F(t — 1) 3'(x) dV + ‘i dx [ R(t — 2) 5'(x) do | 
0 Vv 0 o 


(3) 


(4) 


_ far[Pa—natyav tf acf Ra—2s() a. 
0 Vv 0 o 


| Das ist der Reziprozitatssatz der Dynamik der elastischen Kérper, in welchem, im Gegensatz 
| zu Gleichung (1), keine Tragheitskrafte, sondern nur die gewohnlichen Krafte und die zugeh6- 


-rigen Verschiebungen erscheinen. 


E. H. Love, The mathematical theory of elasticity, Kap. VII Nr. 121, Cambridge 1934. 


1 A, 
2 4, E. H. Love, a. a. O. 
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Betrachten wir jetzt den speziellen Fall, in welchen die Kriafte als Produkt einer Ortsfunk- 
tion mit einer Zeitfunktion G (t), die fiir alle Krafte gleich ist, ausgedriickt werden kénnen, 
namlich 


Fi, P)=Gi)a, F(t, P)=G)@, Re, P)=Ci)b, Rit, P) =H, (5) 


wo @, a’, b, b’ nur Ortsfunktionen sind. Dann wird aus (4) 
t = t 7 
elo ar a s'(t) dV + [63(z) do| = (Cua) ae| Jw 5(r) dV+{ W §(x) do . (6) 
0 V C 0 V 6G 


Da diese Gleichung zu jeder Zeit erfillt sein muB, ist es leicht, zu beweisen, daB die Glieder, die 
in (6) mit G (t— 7) multipliziert sind, gleich sein miissen. Wenn wir diese Glieder mit G (t) 
multiplizieren und auf (5) achten, so haben wir (zu jeder Zeit) 


[Ps'dV+[Rs'do=[F’'sdV+][ RSdo. (7) 
Va o Ve o 


Diese Beziehung ist vollkommen analog zum Bettischen Satz in der Statik der elastischen 
Korper, natiirlich nur fiir die Ansatze (5) der Krafte und mit verschwindenden Anfangsbedin- 
gungen. 

An anderer Stelle habe ich die Ausdehnung der Gleichung (4) auf den Fall mit nicht ver- 
schwindenden Anfangsbedingungen gezeigt und andere Anwendungen angegeben }. 

Wir kehren jetzt zu (7) zuriick und nehmen an, daf das erste Kraftesystem sich auf eine 


einzige, im Punkt A konzentrierte Kraft H(A), das zweite Kraftesystem auf eine einzige, im 


Punkt B konzentrierte Kraft H’ (B) reduziert. Die Ansatze (5) sind bestimmt erfillt, wenn 
beide Krafte zu jeder Zeit den gleichen absoluten Betrag und eine von der Zeit unabhangige 
Angriffslinie haben (die auch verschieden fiir die zwei Krafte sein kann) und gleichzeitig ihre 
Richtung andern. Mit diesen Annahmen wird aus (7) 


H(A) s'(A) = H'(B) s(B), (8) 
wo 5(B) und s’(A) die Verriickungen in B und A infolge der Krafte H(A) und H’(B) sind. 
Die Beziehung (8), die einen bekannten Maxwellschen Satz auf die Dynamik der elastischen 
Korper erweitert, liefert die Komponente in der Richtung von H der Verschiebung, die in A 


infolge H’ entsteht, sobald die Komponente von s (B) in der Richtung von H’(B) bekannt ist ?. 

Endlich beweisen wir, daB die Gleichung (4) auch bei Reibungskraften noch gilt; genauer 
gesagt, wir nehmen an, da8 in den Beziehungen zwischen den Spannungen und den Dehnungen 
Ziahigkeitsglieder vorkommen, die proportional zur Ableitung der Dehnungen sind (Hypothese 
von Voigt), sowie Nachwirkungsglieder mit unveranderlicher Nachwirkung, die, gemaB der 
Theorie von Boltzmann und Volterra, die Erscheinungen der elastischen Hysteresis darstellen. 

Wenn wir fiir t < 0 den Korper im natiirlichen Zustand annehmen, so wird hiernach aus (3) 


Oi; (t) =2 “Ee; (t) +A 0: Q(t) 2 Ub dei (t) NV peo | 


Ot 0 
i ae (9) 
+ J pe 0) e5(2) de + 84 [ yt — 2) OC) dr, 


wo pv’ und A’ neue Konstanten und @ (i — 1), y (t — t) bekannte Funktionen von t und 7 sind. 
Wir wenden auf diese Gleichung die Laplace-Transformation an; nennen wir g(p) und h(p) 
die L-Transformierten von g(t) und y(t), so haben wir 


Lois) =[2 4 +20 p + 8(p)) Slej) + 6; [4 +17 p +h(p)] B(@) ; (10) 


das heiBt: die Beziehung zwischen den Transformierten der Spannungen und der Dehnungen ~ 


ist, bis auf den Wert der Koeffizienten, gleich der Beziehung (2) zwischen den Spannungen und 
den Dehungen in einem elastischen idealen Kérper. Man versteht also leicht, da® auch bei 
Giltigkeit der Gleichung (9) die Gleichung (3) und folglich die Gleichung (4) giltig bleiben, 
womit also die Gleichung (4), wie behauptet, auf den elastischen dissipativen Systemen erwei- 
tert ist. 
(Eingegangen am 1, Juni 1953.) 
Anschrift des Verfassers: Professor D. Graffi, Bologna, Istituto Matematico dell’ Universita. 


1 D. Graffi, Ann. mat. pura appl. (4) 17 (1939) $.173; Mem. Acc. Scienze Bologna (10) 4 (1946/47) S.103. 

* Diesen Satz hat schon Lord Rayleigh, Theory of Sound Nr. 108, ausgesprochen; er setzt aber voraus, 
daB beide Krafte und die Verschiebungen harmonisch von der Zeit abhangen. Hier wird der Satz, mit 
einem anderen Verfahren, auf mit der Zeit beliebig variable Krafte erweitert. 
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Beitrag zur Berechnung schiefwinkliger Platten. 
Von J. Krettner. 


1. Einleitung. Das Problem der schiefen oder parallelogrammférmigen Platte, das in der 
Baustatik, namentlich im Briickenbau eine gewisse Bedeutung hat, ist im wesentlichen erst 
nach dem Kriege, vor allem im Ausland oder von auslindischen Autoren tin Angriff genommen 
worden. Zwei Wege wurden vornehmlich eingeschlagen. Entweder legt man rechtwinklige 
Koordinaten, also auch die gewéhnliche Differentialgleichung fiir die Plattenbiegung zugrunde, 
um dann allerdings infolge der schiefen Randbedingungen cine wesentliche Erschwerung der 
Rechnung in Kauf zu nehmen, wie die Arbeiten von Anzelius? und Rongved? zeigen, wenn 
man nicht andere, z. B. elliptische Funktionen 
heranziehen will, welche in orthogonalen Koor- 
dinaten eine schiefe Periodizitaét aufweisen. 
Andererseits kann man nach dem Vorbild von 
Lardy und Favre von vornherein schiefwink- 
lige, dem Problem angepaBte Koordinaten ver- 
wenden, wobei aber dann in der auf schiefwink- 
lige Koordinaten umgerechneten Differential- 
gleichung gewisse Zusatzglieder auftreten, die 
natiirlich auch wieder rechnerische Schwierig- 
keiten mit sich bringen. Hier soll der Lardy- 
Favresche Weg gewahlt werden, in Ankniipfung 
an meine allgemeinen Ausfiihrungen iiber die 
Darstellung der Tensoren in allgemeinen Koor- 
dinaten *, Durch eine nicht ganz leichte Um- Abb. 1, Zusammenhang zwischen rechtwinkligem und 
formung der bei der Lisung der inhomogenen schiefwinkligem System. 
Differentialgleichung sich ergebenden Koeffi- 
zienten gelingt es — das diirfte zunachst das wesentliche Ergebnis der vorliegenden Arbeit 
sein — mit Hilfe der Ubertragung der bei rechtwinkligen Platten verwendeten Methode 5, die 
auftretenden Doppelreihen fiir die Durchbiegung in einfache, gut konvergierende Reihen zu 
transformieren. In der vorliegenden Mitteilung soll zunachst der Fall der momentenfrei aufge- 
legten, schiefwinkligen, vollbelasteten Platte zugrunde gelegt werden. 


2. Die Differentialgleichung der Biegeflache. Bekanntlich lautet die Plattengleichung in 
rechtwinkligen Koordinaten, wenn B die Biegesteifigkeit der Platte bedeutet 


AAw ==. (1) 


Mittels der Transformationsgleichungen (Abb. 1) 


=f —Vte Pp, Y= 


cos |p 
ergibt sich : zs 
1 f02w 2w ; w 

fee S ' es 3 
Aw = cos? alse 1 Oy? eRe Ox 55) (3) 

und o _ s 

1 [otw 04w Be w : w w 

AAw = 099s Sor F2(1 +2 sin Pasa — 4800 aes eral (4) 


1 P, Lardy, Schweiz. Bauztg. (1949) S, 207 u. 419; H. Favre, Schweiz. Bauztg. (1942) S. 35, ol, 60; 
H. Favre, Bulletin Techn. de la Suisse Romane (1946) S. 321; J. Nielsen, Beitrag zur Berechnung schiefw. 
Platten. Akademie der Techn. Wissenschaften, 

2 A, Anzelius, Bauing. 20 (1939) S. 478—482. 

3 Rongveld, Uber das Problem schiefw. Platten, Diss. Berlin 1945. 

4 J. Krettner, Osterr. Ing.-Arch. 7 (1953) S. 11; J. Krettner, Ing.-Arch. 21 (1953) S. 112—118. 

6 §. Iguchi, Eine Lésung fiir die Berechnung der biegsamen rechteckigen Platte. Berlin 1933; 
W. Miiller, Ing.-Arch. 21 (1952) S. 278. 
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Setzt man s, = sin y und c, = cos g, so erhalt man aus (1) die partielle Differentialgleichung 
vierter Ordnung 


O4w O4w d4w tw \ 
ee ter irre ly tm taaeat ©) 


Man kann nun das Integral dieser Plattengleichung gemiB w = wy + w, aus dem Integral wy 
der vollstindigen Gleichung A4w = p/B und aus der Lisung w, der homogenen Gleichung _ 
AAw = 0 zusammensetzen. Das partikulare Integral der inhomogenen Gleichung erfaBt dabei | 
die einwirkende Oberflachenbelastung der Platte, wahrend das allgemeine Integral der homo- 
genen Gleichung den Randbedingungen des Problems angepaBt werden muB. 


3. Partikulares Integral der vollstandigen Gleichung. Aus dem Verhalten der Platte ist 
ersichtlich, da8 fiir die Funktion der Biegeflache stets wy (a — x, b — y) = wy (x, y) sein muB. 
Gegeniiber der Rechteckplatte, wo doppelte Symmetrie besteht, ist demnach ein erweiterter 
Ansatz notwendig, der in der Form 


W(x, ¥) =D) D (tna cos 7 4 cos y + dn sin" sin 7S (m,7e 1 3 oe) ey 


angesetzt werden kann. Fiir den Druck, der wie im Falle der Rechteckplatte doppelte Symme- 
trie aufweist, ist fiir m,n —1,3,5,... der Ansatz 


; Bis ak IU 
ye > > Joplin “= x sin — y (7) 
zu wahlen. 


Verwendet man zur Abkirzung 
= 2 m\2]2 2 /nm\2 
gal RA earch 
=i ma\ (na\ | (/ma\2 nm\2 
baa MENS) 5): o 


He Bo Co 
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so ergeben die Ansitze (6) und (7) durch Einsetzen in die Differentialgleichung (5) fiir die 
Koeffizienten c,,, und d,,,, die beiden Bestimmungsgleichungen 


@ Cnn +6 dan =0, 


c 


biews +4 dia a 


n 


ast e o a moGlamels | 
mn G2 _— $2 mn (5) a) eal) al t 
eee gg (EN ee 
ied Goo (cis 


Im Falle der Rechteckplatte, fiir die s) = 0 und cy = 1 ist, erhalt man 


mit der Liésung 


(9) 


eer 


(10) 
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4, Umformung der Koeffizienten in dem Integral Wo der Biegefliiche. Das wichtige Problem 
besteht nun darin, die errechneten Koeffizienten von Wy so umzuformen, daB die vorliegenden 
Doppelreihen in einfache Fourierreihen entwickelt werden kénnen. Zunichst ist 


em aera ya ea oe 
a — B 2\a+b ==)" | 
11 
ont (14) 
ey Seer ry aad 
Setzt man } = s, + ic, und 3 = s,—ic,, so ergibt sich fiir @ + 6 und a — b nach einer 
leichten Rechnung 
Eig us) Fs Wire at Re 
(+e) C+ 2) 


(=) oF (ee ey ay 


Damit erhalten wir fiir die Quotienten dieser Ausdriicke 


ve ee 


1 
Gea 9 


a 


a 


7 a 
cotton 
hee 

\ 


io. 1 1 1 ®) (12) 
a+b e2) | ("= , cae max jae ee og nn|[ma Be é 
=) o —oy | [2% + 0% saree ys a 
1 1] 1 1 2 
ab nm\? mx na|? | [ma na? ma na||mo nz| {° (3) 
C—— sided , v Foy, J J | gv 
("F) wo — oy = 87 ae | 072" o 
Fithren wir noch die Abkiirzungen ein 
© = ai O = TID ST (65 Hie) i = Wi — Cy) = nt — On | 
(14) 
o = aie =i 8 = = (8) — § ¢) i = = (is + ¢)) =ynt+6,, 
wobei 
A ==, Yn = Same OIF = ey (15) 
gesetzt wurde, so ergibt sich bei geeigneter Zusammenfassung aus (11), (12) und (13) 
a at ee Oh db: o —m re (16) 
ee (hj —@)? |(m2 x? + @?)? | (m2 a? + 2)? (m2 7? + wm?) (m2? x? + ?)}’ 
b want mM 7 W 4 mm = ms:(o+@) : (17) 
Bp  (d6—oP fe mm + w)2 (m2 a? + &2)2 (m2 x2 + w?) - (m2 x? + 2) 
Die Koeffizienten der Doppelreihe nehmen damit folgende Werte an: 
2a* O5Ct o o o+o (18) 
Cnn = 7 2 1 2 722 272 2 2 2 772) |? 
te (hb —o)? n ce me? += wp)? * (m2 x2 + 7)? (m2 ae? + w*) (m? a? + ’”) 
d att a ow? L ie 20°o 
Pie aeeoale gt) (m2 x? + ~@?)2 © (m am) (m? n% + 2) 9 (mat) (m? a? + w?) (m2? + A?) (19) 
m ot ™m 3 rc ma 7 | 
(mm? a? -—- o7)? (m? 7? + 62)2 | (m2? x? + w?) (m? x? + 65?) 
Wir setzen nun noch : ie 
= Vee a 4 Po a zi) 
ee ¢ re — m3 B 
und fiihren folgende Funktionen ein: 
cosmaé _— cosmamé a cos ma & 
fi () = (m2 me + 2)? ’ fr (E) ao (m? 72 + &?)? ? fs() =S (m? 7? + w?) (m2? 2? + ob”) 
mit 
F,(é) =— 2 Ki [of,(é) + of) — @ +) fs(€)] (20) 
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und 
sin ma é sin ma é sinmaé 
81(5) =( (m x) (m? me Ew)? ? 8() (Gm x) (m2 a? + a)?” &3() = (mz) ) (m2 x? + w?) (m? 72+ 5?) ” 
nmazé (ma)sinmaé& (m x) sin m x & 
aie) = ae 656) =e ane + 86) = Ga at Gata OT 
mit 


Gn(E) = K & [w* g,(€) + &? g(€) — 2 & w g(£) — 84 (6) — 85 (6) + 2 &6()1- (21) 


Damit geht die Gleichung fiir wy itber in 


y(n) = >) (re ee 2a) (22) 


n 


5. Wichtige mathematische Formeln zur Umwandlung der Doppelreihe. Zur Umgestaltung 
der auftretenden Doppelreihen benétigen wir verschiedene Summierungen von Fourier- 
Reihen, die mit Hilfe der bei A. Kneser! angegebenen Formel 

cosm(x—®) 1 , wOojad SS 
>) m? +. a? ~~ o@ 2aGinaz te eeatb apie (23) 
méglich sind. Setzt man in dieser Reihenbeziehung an Stelle von 7 — @ bzw. a die Werte 


7 — =m — bzw. = , so erhalt man die dem Problem angepaBte Form 


cosmaé 1 Coj w (1 — &) 1 tie 
_ m+ wz 2 | wo Sin w =F | Coals tte Eye Ge 
Mit der Beziehung 

i Sop 1 1 35 
(mt a Lo) (mPa EH eer ae oe aa as 

erhalt man aus (24) die fiir uns wichtige Formel 
= cosmaé oo! il Cofa(l—-g) Gof d(l1—é)y] - 1 26 
(m? a2 + w?) (m? a? + o?) 2 \62— w? o Sin w oO: Sing | wart” oe) 


Ersetzt man in dieser Beziehung die Variable € durch 1 — &, so ist cos mz (1 — £) = (— 1)™ 
cos mz &, und Formel (26) geht tiber in 


a, (—1)™ cos ma E =+| ] ee = Goj HE\__ 1 : (27) 
(m? 72 + w®) (m2? a? + 2) 2 |? — w? \w Gih w 6 Gin & (w+ &)? 


o Sin & 


Durch Subtraktion der Reihenformel (27) von (26) erhalt man eine der Formel (26) ent- 
sprechende fiir ungerade m giiltige Bezichung: 


SS cosmaé ro [1 —(—1)™]cosmaz& 
_ (m? a? + @?) (m2 2 + A?) 2 m? 3? + w?) (m? a? + 2) 


28 
1 ae Sin a = =| Sa 
Sgt =o) - mae Gn 13551) 
io) Goji > © So} > 
Differentiiert man diese Formel nach der Variablen é so ergibt sich 
awe tl 
we (mz) sinmaé& _ 1 Gojo(S—s) Gof (5 —¢) 

(m2 a? + w2) (m2 a? + 2) 4 (6? + oe? — er 29 
2 ea Goi 


(m = 1,3,5,...). 


1 A, Kneser, Die Inte 


gralgleichungen und ihre Anwendungen in der mathem, Physik, S. 148, Braun- 
schweig 1911. ; 
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Ersetzt man in der Formel (24) die Variable € durch 1 — &, so geht diese Formel in 
(—1)"cosmaé 1 / Gojwé 1 
= m? 72 +. 2 +(5 Gino | 4 (30) 


tiber. Aus Formel (24) und (30) ergibt sich durch Subtraktion fiir ungerade m 


cosmaE TS {1 —(—1)"] cos maé 1 Sin w (5 —<) 


< (ma? +o) neat sot ee (ME Ly Oy c,) a(S) 


o®j>  (0<é<)). 


Daraus erhalt man unmittelbar durch Differentiation nach (7) 


> cosmae 2 Goiz ee Se) eine (ye 5|| mecha 2 
( 


™m a at 16 w* Coj? > eo) 


Nochmalige Differentiation nach é ergibt 


@ 1 
P 2 £60} + Ginw (5 —é) + Ginwé 
(mz) sinma& 2 2 
(m? x? + op? = 5 (eee Pee Pas ee ore (33) 
16 w Gof? > 


Integriert man Formel (23) auf beiden Seiten nach der Variablen # von @ bis z, so ergibt sich 
eine weitere fiir uns wichtige Reihenbeziehung 


SSeEeu “Sa, Se Cs ae ee a (34) 


m(m? + a?) 2 a? Ginaz | 


Ersetzen wir in dieser Beziehung wieder 2 — # durch a & und az durch , so erhalten wir 


__ sinmaé ] Gin w(1— &) 


(m x) (m? a? Fo) ore es uae (m = 1, 2,3,4,...). (35) 


Mittels der Beziehung (25) erhalt man damit 
sinmaé _1{__1 [Gino—s _ Sin s—~#) + wart (36) 
2h (m? x? + w)(m? a? +62) 2 aaa a? Gin w oO Sin d 
Differentiation von Forme] (35) nach w ergibt 
sin ma & 
2 (m a) (m? 7? +- a”)? 


; ON 
-) ate ae Gin w(1 — &) + w[& Gin w Coj w(1 — €) — Sin w€] o 
= 4 


Sin? w 


2(1 ‘|. 


Ersetzt man nun noch in Formel (37) und (36) € durch 1—é&, und beachtet, daB 


sin mz (1 — &) = — (— 1)” sin mz £ ist, so gehen die Formeln tiber in 
(—1)"sin maé a ee Sees ae) é (38) 
y (m2) (m? 7? + w?) (m2 a? + 5?) ees: o—o' |e’ Sinn A Sind (w+)? } ” 


(—1)"sinma& _ a : ae Gin w Gin w E+ w [Coj w Sin w E— E Gin w Coj w €] —2¢} : (39) 


(m m1) (m? x? + w®)? Sin? w 


Subtraktion der Formel (38) von (36) bzw. (39) von (37) ergibt wieder die entsprechenden 
Formeln fiir ungerade m 


1 yf ‘ 


x sin ma é& ed) a 1 
i (mz) (m? 22 + w)(m?n?-+ 6?) 4 6? — w? 0? Gof 2 r Hg 


52 
sin ma & 
pare mt) (m2 x? + w?)? 
L ] 1 
ota 1 |, 2 Co} ag Gino| 5 —é| + 2 Coj ol -—é)|+0 Ginwé 


4 eh 


(6) 
ae 
4 Coj 5} 
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(41) 


(m = 1,355.4. <4) 


Damit sind die zur Umwandlung der Doppelreihen notwendigen Reihenbeziehungen hergeleitet. 


6. Die Umwandlung der Doppelreihe. Fiihrt man nun in die GréBen F,(€) und G,(€) von 
(20) und (21) die ermittelten Reihenbeziehungen ein, so erhalt man nach langerer Zwischen- 


rechnung und geeigneter Zusammenfassung 


eee 1 4 o ] a 
b Coj-¥ 6of » (> —£) + w Go} 6oj d (F—4) 


Rif) = = 
8 We A o 
| (w+ 0) Coj | bof > 
i220 aS ie | | 2¢ ei ‘les ) 
1 & Gojiz Sin wo\y —é) +o vj Sind \y —é 
l 8 p t 
(o 6)? oj $ Cot 
Oy 1 o 1 ; (22) 
l & Sof > Gino (5 —£) —o6oj > Sind — | 
4 1 
(6 —o) (6: w) Co} F 6oj 5 
: o oD 7 
1 & Sof? gy Sof oF +o Cop? Z Coj o 
16 6\? 
(oe) (Gof 2 05 =) 
fa ‘q , 
oO Coj> Sin w (= 4 +o Coj> Gin ee = 
C () K € 2 2 2 2 
n Tis 4 0 
| ( w) Go} Coj 
0 1 t ] 
1 8 Gof Cojo = =, + 0 Gof F Coj 6 = = 
ei j 
(0+ &) Gof F Gof 
1 OGo2-S Sino § + oGoPS Sin dé 
pat ae (43) 
(+d) (coj F Gof 2) [ 
p 1 1 
: é? Gof Cojo Came — 0 oj F 6j Bea) 
4 2 (ad? — w?) (a0) 
(m+ 0) Cojz Coj > 
bie ae. sar peau | 
1 SF 6oiels 1) wf FOES —2) | og oy 
2 (a? — w?) Co ey 2) 4(m+ o)? 
Z 2 
Dabei ist nach (14) und (15) 
© +@® =27,1, 
o-—-wo =26,, 
O:°o=— ak 
y) 
By : WON ‘ 
db? = (8 —y3) +27 i= (7) [(o8 — 52) + 2s) ey i], (44) 
2 2 
wo” = (dn — Yn) — 2Yn On t oa a [(¢6 — 86) — 2 s9 eg t], 
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Wir fithren nun noch die folgenden Ausdriicke ein: 


M, = Gof 4, & cos yn (1 — &) + Cof 6,(1 — 2) cos », é, 
M, = Sin On é cos Yn (1 — &) — Sin On (1 — é) COS Yn ee | 


M; = (1 +cos Vn Coj On) Coj On € cos y, & + sin yn Stn 6, sin y, € Gin 6, E, 32) 
— M, =(1 +cosy, Coj 6,) Gin 6, cos y, € + sin Yn * Gin 6, sin y, é Coj 6, é, 
N, = Gin 6, é sin y, (1 — &) + Gin 6, (1 — &) sin y, é, 
= Woj 6, & sin y, (1 — &) — of 6, (1 — €) sin y, €, (46) 


Nz = (1 +cosy, €of 6,) Gin 6, &- sin y, € — sin Yn Sint 6, Col 6, &- cos y, €, 

N, = (1 +cos yp, Cof 6,) Cof 6,€+ sin y, € — sin y, Gin 6, Gin 6, + cos Vnés 
Durch Einsetzen dieser Ausdriicke in (42), (43) und (22) erhalt man schlieBlich nach einer 
langeren Zwischenrechnung fiir die Funktion w, die einfache Fourierreihe 


Fn Ma) 
a o.Co23 (SoM, + ¢) Nj). 


ie ear 
Py b4 co | Le M; — cy N3) 
Goji = ee. + cos 8 Sa 


wy(§, 7) = eB =| 


: ; | cos 7 
+ [2 8 N, + sy(3 — 2 s?) N,] Pen yan es 
nd (oj = 7 co +608 > x) 


sus “he N.) (47) 
co a Cote SON a 
| Pea Co (Co M, — sy Nz) + 2 ¢§ (Goj = oe cos Ts) 


mn A 
Go ee Z * ¢9 + cos 5 Sa 


5 
S068 = = 7] | 


— [2 8 M, — s,(3 — 2s?) N,] (n =1,3,5,7°**). 


| n> (oj = & + cos 7. s)| 
Wird der Anfangspunkt des Koordinatensystems in den Punkt (a/2, 0) verschoben, so hat man 
in den Formeln (45), (46) und (47) lediglich an Stelle von € den Wert & + & zu setzen. 

Kommt schlieBlich der Anfangspunkt des Koordinatensystems in den Mittelpunkt der Platte 
zu liegen, so ist noch 7/A durch s+ y/A za ersetzen und zu beachten, daB damit 


aes n+1 
‘ nee?) = Tt Sy oat ee ee or LL 
sin nas + a ( 12 cos yn und cos nals i) =(— 1)? sin 70 (48) 


wird. 


7. Der Spezialfall der rechteckigen Platte. Im Falle der rechteckigen Platte ist p = 0 und 
damit s, = 0 und c, = 1, sowie alle Nx = 0. Aus Formeln (46) ergibt sich ferner 


M, = Gof .6(1 — £) + Gof 6, € =2 Cof Gof 6, pee ‘). 


M, = Gin 6, — Gin 6, (1 — 6) = — 260 F Gin 4, (3 —§). 
= (1 +Gof 6,) Cof 5,€ = = 2Cof? = Cof 4, €, 
M, = — (1 + Goj 6,) Gin 6, & =— 2602+ 6 Gin 6, &. 


Aus (47) erhalten wir damit im Falle der Rechteckplatte 
wy(€, 7) = 

1 1 2 5 Oa 

TE N2E Goi 5G Si n (5 — ) + Gin” g| + 460155 i Gof > (5-2) | sin = 7] 


jl Ss 
5B un 
= ay 


(nee leas Sy Tees) = 129) 
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Wendet man, die unter (41) angegebene Reihenformel auf den Fall der Rechteckplatte an, so 
geht (49) in die wohlbekannte Formel von Navier iiber: 


mun 
sin mz & sin = 7 1 


16 po a4 
wo(5,1) = 7B Paes (m+ n) (m? 2? + w?)?” a 


Verschiebt man den Koordinatenanfangspunkt in den Punkt (a/2, 0) der Platte, so ist wieder & 


durch = + &zu ersetzen und man erhalt in Ubereinstimmung mit Iguchi 


mun un gi) IL, ai un UN >. wn . on 
ae 2 Coj 97 Coj Zé fF an Sin 57 Gof FF —F Gol a7 sin Fal] sin 7 
wo (é, 1) = sB 1— wun | ns 
m 2 oj? OM 
(51) 
Wird der Anfangspunkt noch in den Mittelpunkt der Platte verlegt, so ist in der Formel (51) 


n—1l1 


noch sin =n durch (—1) 2 cos ae zu ersetzen. 


8. Lésung der homogenen Differentialgleichung 14w = 0. Legt man den Koordinaten- 
anfangspunkt in die untere linke Ecke der Platte, so kann das Integral w, der homogenen 
Gleichung bei Beriicksichtigung der Polarsymmetrie zum Mittelpunkt in der Form angegeben 
werden: 


w, (7) = 
= 3 format 2)+9fr—3] fontmna(o—2) ras 4) Smal 3) 


+sin ma|(§ Seer %o(n _ 5 {comin matea(n — 7) +eom(7 — 7) 60 mace —3)] 
{ 
\ 


+2 [mess —a)+ (3) dy» Gof mx e4(6—F) + dan (§—7) Sitt max o(E = 
+sinna]9(§ 3) +(n— Z)| {den SitrnxeE— 3) + dan (§—] Cof jneg(§—F \. 


Die Koeffizienten c;x und d;x lassen sich hier, wie in ahnlichen Problemen (eingespannte Recht- 
ecksplatte usw.) zwar nicht in endlicher Form darstellen, sind aber, wie Lardy fiir den Fall der 
vollbelasteten schiefwinkligen Platte eingehend gezeigt hat aus den Randbedingungen mittels 
Interationsverfahren gut bestimmbar. 

In einer weiteren Arbeit soll auf andere Belastungsfalle und auf die Bestimmung der Mo- 
mente naher eingegangen werden. 


(Eingegangen am 24, Juni 1953.) 
Anschrift des Verfassers: Dr. J. Krettner, Miinchen 13, Gentzstr. 5 III]. 
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Zur Stabilitat der periodischen, erzwungenen Rollschwingungen 
eines Schiffes. 


Von O. Grim, 


Die periodischen erzwungenen Schwingungen eines einfachen Schwingers mit nicht gerader 
Kennlinie, die durch die Lésungen der Differentialgleichung 


% + f(x) = P sin (wt) (1) 
beschrieben werden, sind nicht immer stabil. Fiir bestimmte Wertepaare P, w gibt es drei ver- 
schiedene Lésungen dieser Differentialgleichung, von denen eine Lésung instabil ist und daher 
nicht verwirklicht werden kann. 

Die bekannten, bisher durchgefithrten Untersuchungen der Stabilitat dieser Schwingungen! 
gelten aber nur fiir kleine Abweichungen der Federcharakteristik f(x) von dem linearen Verlauf. 
Es gibt indessen Faille, in denen die Abweichungen sehr gro8 sind. Ein solcher Fall liegt z. B. 
bei den Rollschwingungen eines Schiffes, fiir die die Bewegungsgleichung 


0% +EK—) =M (2) 


lautet, vor. Das Stabilitatsmoment Gh, das bestrebt ist, das gekrangte Schiff wieder aufzu- 
richten, ist nur fiir kleine Winkel » proportional zu diesem Winkel. Es gibt aber auch Roll- 
schwingungen mit Amplituden, fiir die das lineare Gesetz fiir das Riickstellmoment langst nicht 
mehr gilt. Und da die GréBe des Stabilitatsmomentes so bemessen werden muf, da das Schiff 
niemals in die Gefahr des Kenterns kommen kann, so ist der ganze Bereich bis zu dem stati- 
schen Kenterwinkel g,, fiir den das Stabilitatsmoment zu Null wird, von praktischer Bedeu- 
tung (Abb. 1). 

Nun geniigt das bisher bekannte Ergebnis der erwaihnten Untersuchungen, namlich da8 in 
dem Fall, daB es drei verschiedene Lisungen der Bewegungsgleichung gibt, die beiden Lé- 
sungen mit dem kleinsten und dem gréBten 
Amplitudenwert von @ stabil sind und die dritte 
Lésung instabil ist, fiir das Beispiel des Schiffes 


Glug) 


Abb. 1. Abb. 2. 


nicht: Wenn z. B. die Kurven fir die Lésungen der Bewegungsgleichung (ohne Dampfung) des 
Schiffes in dem (gy, T)-Diagramm Abb. 2 bis zu der Periode T = ov gefiihrt werden, miissen 
die Kurven fiir T = 00 in die drei Punkte einmiinden, die auch als statische Gleichgewichts- 
lagen unmittelbar aus dem Diagramm fiir das Stabilititsmoment gewonnen werden konnen 
(Abb. 1). Wenn das Ergebnis der Untersuchungen der Stabilitat der erzwungenen Schwin- 
gungen auf diesen Fall ausgedehnt werden diirfte, miiBte der jeweils oberste und unterste 
Ast der in Abb. 2 gezeichneten Kurve stabil sein. Die Betrachtung des Verlaufs des Stabili- 
tatsmomentes in Abb. | zeigt aber sofort, daB bei statischer Beanspruchung, also fiir T = oo, 
nur die Gleichgewichtslage (a) stabil sein kann. Daraus muB geschlossen werden, da8B auch 
der oberste Ast in dem Diagramm Abb. 2 nicht nur stabile Bewegungszustande, sondern von 
einem noch unbekannten Punkt an auch instabile Bewegungszustinde umfafit. Es ist offen- 


1 Siehe z. B. J. J. Stoker, Non-linear vibrations. S. 114 u. 213 ff., New York 1950. 
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sichtlich, daB fiir das Beispiel des Schiffes der vollstindigen Lésung des Problems eine prak- 
tische Bedeutung zukommt. 

Die bekannten Stabilitatsuntersuchungen wurden auf dem folgenden Weg gefiihrt: Die 
periodischen, erzwungenen Schwingungen des Schiffes sind durch die Lésungen der Bewegungs- 
gleichung (2) gegeben. Wenn nun eine kleine Abweichung p von der periodischen Lésung 
irgendwie zustandekommt, gilt hierfiir 


OP +4) +Ghy +y) =M (3) 
oder, wenn man hiervon (2) abzieht 
Oy +Clk(p +p) — Ky)] = 0. (4) 


Diese Differentialgleichung fiir », in die man fiir 7 die aus (2) resultierende periodische Lésung, 
deren Stabilitat untersucht werden soll, einzusetzen hat, ist durch periodisch veranderliche 
Koeffizienten gekennzeichnet. Fir diese Differentialgleichung existieren stabile und instabile 
Bereiche, in denen die Lisung y entweder mit der Zeit abklingt und verschwindet oder unbe- 
schrankt anwachst. Durch die Feststellung, ob fiir eine bestimmte Lésung die Lésung der 
Differential gleichung (4) in einen stabilen oder instabilen Bereich fallt, ist die Frage nach der 
Stabilitat der Lisung  entschieden. 

Es entsteht also nun die Frage, ob dieser Weg in dem Beispiel des Schiffes auch bis zu so 
grofen Winkeln q, fiir die die Hebelarmkurve h(y) den maximalen Wert iiberschritten hat und 
sich wieder der Null nahert, beschritten werden kann und ein brauchbares Ergebnis liefert. 
Diese Frage kann ohne Schwierigkeiten beantwortet werden, wie folgt. 

Fiir diese Untersuchung wird eine einfache Hebelarmkurve angenommen, die keinen Wende- 
punkt besitzt und die mathematisch durch 


% Hy) = oY (2 ae (5) 
ausgedriickt werden kann. AufSerdem wird zunachst angenommen, dafs die erzwungene Roll- 
schwingung genau sinusformig verlauft: 

~Y =f sin (wt), 
d. h. also, daB das erregende Moment durch den aus (2) hierbei bestimmten Wert 


M a x 3 m2\ Org. . 
ow, ( wo? |, ji we = sin (wt) + @ =e i sin (3%) 
Pr Pre 


4 


gegeben ist. Ob diese Bewegung @ sin (wt) stabil und daher verwirklichbar ist, folgt aus der 
Gleichung (4), die in diesem Fal] die Form 


Res , 3 , 3G 
yy +o] r 9 | Pp 


annimmt. Das ist eine Mathieusche Differentialgleichung. Wenn diese Gleichung in die Nor- 
malform 


-cos(2@ 7) = (6) 


py +y(A +ycos x) =0 (7) 
gebracht wird, sind die Parameter 4 und y 
Ir Gaye 8 GE 3 w2 @ 
veer snee te jp POPE 
to (1 2 a Pat ee ees p2 (8) 


Fir die Mathieusche Differentialgleichung enthalt die bekannte Ince-Struttsche Karte die 
Grenzkurven y = f(A), die die stabilen Bereiche der Lésung von den instabilen Bereichen 
trennen. Da fiir jedes Wertepaar /, y aus (8) ein Wertepaar y/o, P/~, bestimmt werden kann 
kénnen die Grenzkurven der Ince-Struttschen Karte in das Diagramm (w/o, @/@x) ba 
(T/T, y/g.) tbertragen werden (Abb. 3). Dieses Diagramm zeigt mere hee und anechams 
lich, fiir welche Wertepaare w, y die erzwungene periodische Schwingung @ sin (wt) stabil oder 
instabilist. Das Diagramm zeigt, daB tatsachlich der oberste Ast der Kurve fiir die erzwungene 
Schwingung (Abb. 2) nicht fiir alle Perioden stabil ist. Das Diagramm 3 enthalt auch die Kad 


; Se 
Pe et a so Fe ; : . ° 
oO” = Wo ( ease we 7 5 fiir die das erregende Moment kein Glied mit der Kreisfrequenz (a) enthalt. 
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Die obere der beiden von T/T, =1 ausgehenden Grenzkurven (b) 
Winkel y der Kurve w? = o (1 — = 
PR 

Grenzkurve und die Feststellung, daB die In- 
stabilitat der erzwungenen Schwingung schon 
bei Winkeln unterhalb des Kenterwinkels gp, 

_ beginnen kann. | 


schmiegt sich fiir kleine 


eng an. Bemerkenswert ist die mit (a) bezeichnete 


Das Diagramm Abb. 3 gilt streng nur fiir 
die rein sinusférmige erzwungene Schwingung. 
Es wird vermutlich fiir ein rein sinusférmiges 
erregendes Moment nicht viel anders aussehen. 
Um zu einem Urteil hieriiber zu kommen, wird 
zunachst als erster Schritt in die Bewegungs- 

- gleichung (2) der Ansatz 


P =, sin (wt) + G2 sin (Bt) (9) 
eingefiihrt und das Verhaltnis der Amplituden 
Q2/P1 so gewahlt, da fiir das erregende Moment 


das Glied mit der Kreisfrequenz 3 @ wegfallt. : Abb. 3. 
Daraus folgt das Amplitudenverhaltnis 


G2 _ 1 GF I 
Pr 4 Vi, 9 a it 2 3 Pi 
we 2 QE 
und fiir die Untersuchung der Stabilitat dieser Bewegung an Stelle von (6) die Differential- 
gleichung 


(10) 


- 2 3 Qi 3 Pi Pz Pi P 
+ pw} + (1 its 2@t) +3! cos(4mt)| = 0. 
yry@ Paeeele oca9 a, cos(2@ t) 4 pa cos (4m t) (11) 


Diese Rechnungen zeigen, daB der Einflu8 des Gliedes @, sin (3et) fiir groBe Frequenzen w bzw. 
groBe Verhaltnisse w*/w% sowie fiir kleine Winkel ~ bzw. kleine Verhaltnisse (,/q, sehr klein sein 
mu, Daraus folgt, daB fiir ein rein sinusférmiges Moment die Grenzkurve (a) in dem (g, T)- 
Diagramm bei kleinen Perioden ebenso verlaufen mu wie in dem Diagramm Abb. 3; das 
Gleiche muB fiir die Grenzkurven (6) und (c) fiir sehr kleine Winkel gelten. 

Es sind aber noch weitergehendere Folgerungen méglich: Wenn in (2) das Moment M = 0 
gesetzt wird, so liefert diese Gleichung die freie Schwingung, die natiirlich periodisch verlaufen 
mu. Durch Differentiation dieser Gleichung erhalt man 

Of +E pW) =0. (12) 
Diese Differential gleichung ist aber identisch mit (4), wenn ~ durch y ersetzt wird, wenn (4) auf 
die freie Schwingung angewandt und yw als sehr k!ein angesehen wird. Da nun in (12) die aus der 
freien Schwingung folgende Lisung @ ein periodischer und damit beschrankter Wert ist, mu} 
auch yp in (4) fiir diesen Fall ein periodischer und beschrankter Wert sein, wobei dieser Wert p 
zwar ebenso wie ~ von der Zeit abhangig sein muB, aber als sehr kleine GréBe angesehen werden 
kann. Das beweist, daB, gleichgiiltig wie die Hebelarmkurve h(y) verlauft, die Kurve in dem 
(yp, T)-Diagramm, die die freie Schwingung bezeichnet und der sich die beiden oberen Aste der 
Amplitudenkurve fiir die erzwungene Schwingung nahern, gleichzeitig Grenzkurve zwischen 
einem stabilen und einem instabilen Bereich sein mu&. Auf einem ahnlichen Wege kann ein 
allgemein giiltiges Gesetz fiir die Lage der Grenzkurve (c) gefunden werden. An der Stelle, an 
der die Amplitudenkurve in dem (gy, T)-Diagramm eine vertikale Tangente hat, mu8 sowohl 
die Bewegungsgleichung (2) als auch die Gleichung (3) eine rein periodische Lésung haben. 
Denn man kann in (3) die Unbekannte (py +-y) auch als auf der Amplitudenkurve liegende eng 
benachbarte Lisung von (2) auffassen. Daraus folgt, daB, gleichgiiltig wie die Hebelarmkurve 
_ verlauft, eine Grenzkurve durch die Punkte, fiir die die Amplitudenkurven lotrechte Tangenten 
_ haben, gehen muf. 
Es ist also sicher, da® die Grenzkurve (6) fir ein rein sinusférmiges Moment bei groBen 
- Winkeln tiefer liegt als in Diagramm Abb. 3, Und es ist damit nur erwiesen, daB, falls drei 
Lisungen der Bewegungsgleichung (2) méglich sind, die Lésung mit der mittleren Amplitude 
- fiir alle Perioden instabil ist, bzw. daGB ein ahnliches Gesetz auch gilt, wenn mehr als drei Lé- 
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sungen miglich sind. Es mu8 nun aber weiter gefolgert werden, daB die Grenzkurve (a) des 
Diagramms Abb. 3 auch fiir das rein sinusférmige erregende Moment existiert. Es ist ja auch 
schon bewiesen, da& sie fiir sehr kleine Perioden T ebenso liegen muB wie in dem Diagramm 
Abb. 3, und es liegt nahe, ihre Lage so anzunehmen, daB sie sich der nun etwas tiefer — nam- 
lich auf der Skelettlinie — liegenden Grenzkurve (6) ebenso nahert wie in dem Diagramm 
Abb. 3 (Abb. 4). Es erscheint daher sicher, da® auch fiir ein rein sinusférmiges Moment der 


Ss 


id. 
05 SS 05 
£instabil 
i 
é 
ha 
[7 
{] 
/ af 
I 
0 F 2 0 
Abb. 4. Abb. 5. 


obere Ast der Amplitudenkurve fiir die erzwungene Schwingung nicht bei allen Perioden stabile 
Bewegungszustande erfaBt, und insbesondere, da der instabile Bereich schon unterhalb des 
statischen Kenterwinkels g, einsetzt (Abb. 5). Damit erscheint das Problem als geliést. Es 
wurden zwar nicht die genauen Lagen der Grenzkurven bestimmt, aber die Genauigkeit, mit der 
sie gezeichnet werden kénnen, erscheint fiir das Beispiel des Schiffes ausreichend. 

Die Lésung des Problems ist auch fiir eine beliebige andere Hebelarmkurve leicht zu finden. 
Zumeist ist an Stelle von (5) der Ansatz 


hy) = ag + by — cy 


passender. Fiir die Lage der Grenzkurven (b) und (c) gelten die oben hierfiir festgestellten Ge- 
setze, und fiir die Grenzkurve (a) kann der Wert fiir T = 0 durch Annahme einer sinusfér- 


Abb. 6. 


migen Schwingung @ sin (wt) bestimmt werden. An Stelle der Differentialgleichung (4) er: 
halt man 


; Mes eee 302, eS co 
Oy+Cy aby bp — Zep —(z bp" — Fe!) cos 21) —F cpt cos (4001) = 0, 


und fitr @ = oo erhalt man die Lage der Grenzkurve (a) hiermit aus 


a +S bp —2cpt=0. 


as 
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Dieser Winkel ist, solange a positiv bleibt, immer kleiner als der statische Kenterwinkel Py: 
fiir den 


a + boy; — co, = 0 
gilt (Abb. 6). Es beginnt also auch in diesem Fall der instabile Bereich unterhalb des statischen 
Kenterwinkels @,. 

Kinschrankend mu erwahnt werden, da in der fiir die Stabilitat der Bewegung maB- 
gebenden, nicht linearen Differentialgleichung (4) bei den weiteren Untersuchungen nur die 
linearen Teile verwendet wurden. Diese Untersuchungen gelten daher streng genommen nur 
fiir sehr kleine Abweichungen von dem untersuchten Bewegungszustand. Es ist jedoch wahr- 


 scheinlich, daB die gewonnenen Erkenntnisse durch diese Einschrankung nicht geschmilert 


werden. 


(Eingegangen am 26. Juni 1953.) 


Anschrift des Verfassers: Dipl.-Ing. O. Grim, Hamburg 1, Berlinertor 21. 
Hamburgische Schiffbauversuchsanstalt 
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Zur Theorie der Vierpilzplatte. 
Von W. Miiller. 


1. Einleitung. Die in einigen fritheren Mitteilungen gegebene Darstellung + einiger Falle 
von rechteckigen Platten soll im folgenden durch die Behandlung einer Platte erganzt werden, 
die von vier paarweise symmetrisch zu den Mittellinien gelegenen Einzelkraften belastet baw. 
gestiitzt wird. Wir wollen die so gekennzeichneten Falle wegen der Verwandtschaft mit dem 
Pilzdeckenfall kurz unter der Bezeichnung Vierpilzplatte zusammenfassen *, Um die Unter- 
suchung mdglichst natiirlich und systematisch durchzufiihren, werden wir dabei von der 
Annahme rechteckiger Druck- bzw. Stittzflachen und von der Darstellung der Durchbiegung 
durch Doppelreihen ausgehen, die nach V. Lewe * mit Hilfe der Fourierschen Methode ohne 
Schwierigkeiten aufgestellt werden kénnen. Zur Weiterbehandlung dieser Reihen, also vor 
allem zum Zwecke der Umwandlung in einfache Reihen sollen zunichst einige allgemeine 
Reihenentwicklungen betrachtet werden, die auch fiir viele andere Falle aus der Plattentheorie 
von Nutzen sind. 


2. Einige fiir die Plattentheorie wichtige Reihenformeln. a) Wir gehen aus von zwei bei 
A. Kneser 4 vorkommenden, bereits an anderer Stelle herangezogenen Formeln. 


a) Zunachst haben wir 


3 cosny _ 2 Goi w(x —y) d (0<y<22) (1) 


n? + q@? 20 Ginw =a 202 


i et LRP BP See 


worin @ eine von n unabhangige Konstante bedeutet. Differentiieren wir nach w, so entsteht 
die weitere Formel 


cosny  a[wnlojay + Cojo(x—y) Ginoxz+oy Ginwy Ginw (a —y)] 1 9 
> wa 43 Gin? wa ~ 2@4° (2) 


Zur Bestatigung fithren wir den Grenziibergang w — 0 durch, wozu die Benutzung der drei 
ersten Glieder in den Entwickelungen 


Gine alter tpytt- |, Cojx=1+ > P45 m4... (3) 


erforderlich ist, weil sich die Glieder mit w? und wtim Zahler aufheben. Es entsteht dann die 
bekannte Fouriersche Formel 


cosy aD ey Lo Deen I 3 rs; 
oder mit y = 0 
cosna@ 4/1 Lie, UP den Lin 
vas na (34 So oe 6) - ) 


Durch Differentiation der Formel (2) nach y resultiert dagegen 


>) nsinny _ aly Gin a at Gof w(x —y) — ax Ginw y)} 
(n? + w?)? 4a Siw wa : (5) 


1 W. Miller, Ing.-Arch. 20 (1952), S. 278; 21 (1953), S. 63: Osterr. Inc.- 
2 F. Télke, Ing.-Arch, 5 (1934), S. 187. nai amma soe ec: 
: V. Lewe, Bauingenie ur 3 (1922), S. 314 u. 344; 4 (1923), S. 453. 
A. Kneser, Die Integralgleichungen und ihre Anwendungen in der math. Physik, $.157 u.167, Braun- 


Met ce vgl. auch S. Iguchi, Kine Lésung fiir die Berechnung biegsamer Platten, S.15 u. 46f., 
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Fiir den Grenziibergang w — 0 braucht man die beiden ersten Glieder von (3) und erhalt dann 


Oy at = hy (a7) 2a —) | 


3 
n n 
oder 


sin n 7 1 (9) 
> Staee 1 tg @—9). | 


n 


b) Wenn nur iiber die ungerade Zahlenfolge summiert wird, geht man aus von der Grund- 
formel 


: ‘x 
cos ny a Sin w (3 = ») 
n2 + @? =o wn (0 << y <a zt) ; (7) 


w= 1, 8) 5s. 4 Coj = 


Differentiation nach w liefert alsdann 


© GR F % . 
cos ny 2 a Go} 97 E y Coj @ ie = ») + Sinw 2 — ») —w 2 Cojwy 
(n? + w?)? << ness 3 (8) 
n=1, 3,5 16 w? Gof? —- 


Weitere Differentiation nach y gibt die Formel 


nsinny 


Quy Coj = sin o(3 —>) +2 CGinwy 
(n? + w?)? on 


(9) 
Tes s15 16 w Cop = 


Fiir w@ — 0 muB man in der Reihe (9) nur die ersten Glieder der Entwicklungen (3), in der 
Reihe (8) dagegen die beiden ersten Glieder benutzen. Dann ergibt sich durch eine leichte 
Rechnung 


4 7 4 
a ae 69? | 40%) = (1— 2 oe) (1 + 2¢— 20%), (10) 
n=1,3, 5 
sn nne — = 9(1—o). (11) 
n=1, 3,5 


c) Subtrahiert man nun zwei entsprechende, fiir die gesamte Zahlenfolge (n = 1, 2,3...) 
und fiir die Folge der ungeraden Zahlen (1, 3,5...) giiltigen Reihen, so ergeben sich weitere 
spater zu verwendende Entwicklungen, bei denen n die gerade Zahlenfolge (n = 2, 4, 6 usw.) 
durchlauft. So hat man 


> cosny ] 
n=2,4, oe es ap | 
oe te i oie | | U 12 
oaGohar +2 eins" i o(5 —7) +207 Gino eine(y—7)| : (12) 
= 7 3 ea ° 
o) 4 
16 w? Gin? = | = | 
: x [27 Gin %2 Goi o(% —7)— 2 Gino 7| 
Ss eee 2 2 (0<y<z). (13) 
n=2, 4, aan ee, 16 0 Sin? = 
Im besonderen wird 
pin One 14 
Sp ap CANO) (laa 20) 3 (14) 
n=2, 4,6... 
cosnao 4/1 5 Sea 15 
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I. Die aufliegende Vierpilzplatte mit der Randbedingung M =0 


3. Die Doppelreihe fiir die Durchbiegung der momentenfrei aufgelegten rechteckigen Vier- 
pilzplatte. Wir setzen eine rechteckige Platte mit den Seitenlangen a und 6 voraus, die auf 
vier symmetrisch zu den Mittellinien und parallel 
zum Grundrechteck gelegenen Rechteckflachen mit 
den Seiten 2% a und 2 f a gleichmaBig belastet und 
an den Randern so aufgelegt ist, daB w und die Mo- 
mentensumme 

E h® 

M =— N Aw (Sareea (16) 
verschwinden (Abb. 1). Die Koordinatenachsen x, y 
sollen mit zwei Rechteckseiten zusammenfallen und 
der Anfangspunkt in die linke untere Ecke des Recht- 
ecks gelegt werden: die Mittelpunkte der Lastflachen 
mégen die Koordinaten %, ¥o3 @ — %q Yo3 Xp. o — Yo» 
a— %; b — y, haben. Machen wir fiir die Belastungs- 
funktion den Ansatz 


P(%, ¥) = af (*) aly), (17) 
wo f (x) und g(y) als Sinusreihen 


f(z) = > 4, sin , 


Abb. 1. Rechteckplatte mit vier rechteckigen Last- ee GUS 
oder Stiitzflachen. 8(y) a B, ae b (18) 
n 


darstellbar sind, so ergeben sich die Koeffizienten A,, und B, nach der Fourierschen Methode. 
Wenn beriicksichtigt wird, da®B f und g innerhalb der Druckflache gleich Kins und auBerhalb 


gleich Null sind, so erhalt man mit Einfithrung von & =—, =2, 4 ale 
a a a 
a 
2 : co —x)— 
2S 6 oats [ (=) sin ma = de =| cos" epee ke [eos m= ee 
a a m potas m ch PS SE 
4 : : 
= —— (1 — cos mz) sinmaasin mz &, ; (19) 
B =—4 (1~— cos nz) sinnz® sinna ® 
ine Pathos A A 


Man sieht daher, daB die Koeffizienten mit geradem Zeiger m oder n verschwinden, oder es wird 


ee : on. : 
A, =—— sin ma asin ma &,, B, =—-sin na sin nx (m,n =1,3,5..). (20) 
Man erhAlt also fiir die Druckfunktion die Doppelreihe 
64 q Le : : : : 
p(xy) = a ese yin ma & sin m7 & sin mz & sin n A sin n nt sin 1 at (21) 


, (msn = 1,3,5.4.) 
Benutzt mau ferner die Differential gleichung fiir die Durchbiegung w 


™ n 


NAA w= ps 
so ergibt sich durch Koeffizientenvergleich fiir die elastische Flache die Doppelreihe 
— 64qb4 sin maa sin mx Sysin ma Esin na sinn asin n a1 
~ Nae > m n(n? + m2 22)? (22) 


RS Fa Beas BAS en 
die in dem m der Folge vorausgesetzten Fall konzentrierter Punktbelastung in 


16 P33 a sin ma é,sin ma Esinn asin nat 


Ena (m2 m® a2) (23) 


n 
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tibergeht. Ohne uns auf die vielen in (22) und (23) enthaltenden Sonderfalle einzulassen, 
bemerken wir nur, daB, wenn Uberdeckungen der Druckflachen vermieden werden sollen, 
folgende Bedingungen erfiillt werden miisser: 


1 B 1 1 1 
0<a<—, a a a + & Sz: B+ Sz° 
Wenn «= &=6/A=n,)/A=1/4 wird, soerhalt manden Navierschen Fall der vollbelasteten Platte. 


4, Darstellung von w durch eine einfache Reihe. Mit Hilfe der in Ziff. 2 gegebenen Reihe (8) 
kann man nun ohne Schwierigkeit etwa die auf den Zeiger n sich bezichende y-Reihe summieren, 
wenn man zunachst unter der Voraussetzung 7 < 7, fiir das Produkt der sin-Funktionen die 
halbe Differenz zweier cos-Funktionen einsetzt, also 

sin noe sin n sae =z (cs ne — cos nas 
schreibt und @ durch m, ferner @ durch (y, — 7)/A baw. (7 + )/A ersetzt. Es ergibt sich 
dann im Gebiete I (0< 7 < ) folgende Reihendarstellung von w: 


- 12 ae sinmaésinmaé [mai Cinman, Sinman 
fy — ee 
Yeas, 5 m3 oj mae Coj m 7 
: ; A 
+2mnan Cinmay Gin ma — no) (24) 


+ 2Goj ma(5 — np) (Sin may — man Cof ay : 


Die Durchbiegung wyim Gebiet II (nj) Sy SA 
— 7) folgt aus dem Ausdruck (24) durch Ver- 
tauschung von 7 und 7, wahrend wy im Ge- 
biet IIT (A — yn) Sn < A) aus wy durch Erset- 
zung von 7 durch 4 — 7 hervorgeht (vgl. Abb. 2). 
An den Ubergangsstellen verhalten sich, wie 
man leicht zeigen kann, w und die erste Ablei- 
tung von w stetig. 


a 


7 
0 


Abb, 3. Hiuhenlinien der Biegeflaiche einer aufgelegten Vier- 
lastplatte (A = 1,5, & = 1/3, 1) = 1/2). 


Abb. 2. Gebietseinteilung fir die Berechnung der 
Vierpilzplatte mit Punktlasten oder Punktsttitzen. 


Natiirlich ist es auch mdglich, die Reihe (24) auf dem Wege der Uberlagerung aus der friither 
behandelten Einpilzplatte herzuleiten, wenn man als Lastpunkte die Punkte x» Yo 3 @ — %q »¥o ; 
b— y, 3 a— %),6— Yo zugrundelegt. 
‘4 Te Nas 3 eben wir die Hohenlinien der Biegeflache fiir den Fall einer Platte aufgezeichnet, 
deren mafgebende Konstanten die Werte & = 1/3, NM =1/2, A = 1,5 haben. Man erkennt hier 
wieder, daB die Lastpunkte fiir die elastische Fliache keine singularen Stellen sind. 
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Wenn man den Navierschen Fall der gleichmafig belasteten Platte mit dem vorliegenden 
Fall verbindet und gleichzeitig die Punktlasten umkehrt, also als Stiitzkrafte einfiihrt, so 
erhailt man die momentenfrei aufliegende belastete und gestiitzte Vierpilzplatte. Die Durch- 
biegung im Navierschen Fall 1aBt sich aber folgendermafen als cinfache Reihe darstellen: 


= 74 Ji}y Ol sinmaé , A 
1 ai errr an z] 


m=1, 3,5... m®> Co) m ‘eos 


mah Ginma 
2 


(25) 


7 : A 
7 +man Sin mo(n =} : 
Co} ™ ey. J 
Bei gleichzeitiger Punktstiitzung erhalt man also die Durchbiegung im Gebiet I in der Form 
: (2w,.)r =e Ut. 

5. Berechnung der Momente fiir die Vierpilzplatte. Dic Biegungsmomente m, und m, und 

das Drillungsmoment m,,, lassen sich, wie bekannt, in der Form darstellen 
Ge 0) ae @ ew ge) ew 

= == eaalle oa eee (ee 26 
wenn man unter » die peas Konstante versteht. Sehr einfach gestaltet sich die Berech- 
nung der mit 1 + y dividierten Momentensumme. Es ergibt sich der Ausdruck 


n in ma E Co} m Mees about. 
M _ My a ity oe 4 P a 4 bass 2 i: vas ‘ : "\2 iG % is be (27) 
ee =, i eyatereert m Cos M4 ° 


in dem allerdings zunichst die Singularitaten verdeckt erscheinen. Es laBt sich dann zeigen, 
daB die Biegungsmomente auch in folgender Weise aus dem Moment ableitbar sind 1: 


1l+y 1 oM oM oM 
m, an 9 M 9 ql \(n on 0 any at d is i) | 


(28) 


4 aM 
Um das Moment M durch #-Funktionen darzustellen, ie wir die komplexe und die kon- 


jugiert komplexe Variable € = &+in,f =&— i7 ein und gehen aus von dem Fall der 
Kinlastplatte, fiir den folgender Ausdruck gefunden wurde ?: 


: Ce oa) 
= Hie! : —, 
ti tonne Cec re Ce ER ot 
6,(—+ 0,( 5 
wenn unter i, der Realteil des Logarithmus des #-Quotienten verstanden wird. Nach der 


Vorschrift dieses Ausdrucks findet man dann auch die von den drei anderen Kraften her- 
rithrenden Momente, wenn man statt €,) nach der Reihe folgende Koordinaten einsetzt: 


Rae a CMe how kta an 4 


Es wird dann 


is 1 —C Cael sp & 
raat l= 
M, = 5— Reln = ro 
1 ome €é—1+4, 
: 2 i 2 
a : fe hg 
‘ gs a= aaa 
M, == Reln : ; (30) 
2 x pe ey 
i 2 1 
p(Cti=b— id i este ) 


1 
2a Helm eg en ale i arceng memes 


1 Vel, i ‘Nadai, Elastische Platten, S. 86, Berlin 1925; K. Grein, Pilzdecken, S. 37, Berlin 1948, 
* Vel. W. Miiller, Ing.-Arch. 21 (1953), S. 64. Rie 
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Wendet man nun folgende, aus der Theorie der 3-Funktionen bekannten Transformations- 
formeln } an: 


(2 £5) = + Alo). 


LA ; a 


Ady Fig) = beet 0), 


1 A : 
feti (vk — 14) = bar ome a(n), 
so heben sich, wie man sieht, die Glieder auBerhalb der Logarithmen heraus, und die Momente 
stellen sich dann iibersichtlich dar als Logarithmen von Quotienten der vier Funktionen 


(32) 


M, a HR, In 


0, (555), 8), 
Im Ganzen 1aBt sich das Moment bei der Vierpilzplatte in die Form 
M = 
$+h\_ (E—S — ( —E be —E 
a5) 5) (4) 55) 58) 58) 8) 58 


of GE) on(ESe) 0, (EE%) 0, (E58) 0, (3%) 0 ESS) 0,655) 465%) 
(33) 


se 
Pd ae 


HR, In 


setzen. Durch Anwendung der Beziehung 
: il 


o(v) By(v) (rv) Fs(v) = 5— A(2 v) B, (0) (34) 
kann man den Ausdruck noch etwas umgestalten. 
Um die Absolutbetrage zahlenmaBig zu berechnen, setzen wir 


E+ =n, E—f =m, Nt+M=%s> N—-N = 82> 


De aig? To + 15S, DeTiaeetts5 Ye + Us, 
eee me 8 ge 
Dann wird : . 
Cee aie Ne a 
pee te pe Es phe Ga ety a U5 


Fiir die Absolutbetrage der 9-Funktionen benutzt man folgende, frither abgeleiteten Naherungs- 
formeln ? : 
|9,(v) |? ~ 4 gt? (Coj zs — cos xr) & —g—2qcosar Coja 7 : 


| B,(v)|? ~ 4 gi/? (Coj ws + cos mr) a —g+t+2qcosar Co] x : , (35) 


| P(v) |? ~ 1— 4qceosarCofas + 24? (cos2ar+ Coj2zs), 
\9,(v)|2 ~ 1+ 4qceosarCojas + 24? (cos2ar-+ Coj2xs), 


1Vgl. G. H. Halphen, Traité des fonctions elliptiques, Bd. I, S. 254, Paris 1886—91; ferner F’. Tricomi, 


Elliptische Funktionen, S. 156, Leipzig 1948. 
2 W. Miiller, Ing.-Arch. 21 (1953), 5S. 65, 66. 
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wo noch die Werte v, v, vz v, und die entsprechenden Indizes der GréBen r und s einzusetzen sind 
und in unserem Falle (A = 1,5) die auftretenden Konstanten folgende Werte annehmen: 
gq =e—74 = 0,00898, 4q'/2 =0,3791, 4q =0,0359, 24? = 0,0001614. 

Mit Hilfe dieser Formeln haben wir nun die Absolutbetrage der vier 0-Funktionen fir vierzig 
Argumente (€ = 0, 0,1, 0,2 bis 0,5i7u. 0 bis 0,75) und danach die Momentenbeiwerte cy = 
M 2a/P berechnet. Um den Vergleich mit der Einlastplatte, deren Lastpunkt im Mittel punkt 
liegt, durchzufiihren, haben wir die mit 4 dividierten c,,-Werte mit den entsprechenden 
¢,7 Werten der Einpilzplatte in den Tabellen 1 und 2 zusammengestellt. 


Tabelle 1. cm-Werte der momentenfrei aufgelegten Vierpilzplatte. 


U/] 

a 0,1 0,2 | 0,3 0,4 0,5 0,6 | 0,7 0,75 
0,1 Ont 0,355 0,537 0,755 0,903 0,996 - 0,895 0,888 
0,2 0,315 0,601 1,038 oO 2 1,832 1,746 1,534 1,450 
0,3 0.555 0,946 1,550 2,477 3,233 3,008 2,495 2,459 
0,4 0 586 OD 1 788 2 665 3,536 3,159 2,897 2,824 
0,5 0,599 1,052 1,689 2,342 2,805 2,888 2°857 2,847 

Tabelle 2. 
Cm- Werte der in der Mitte belasteten Einlastplatte zum Vergleich mit der gleichbelasteten Vierlastplatte. 
U7 | 

mae 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,75 
0,1 0,146 0,304 0,481 0,737 0,898 | 15102 22 1,238 
0,2 0,250 0,586 0,934 1,348 1,810 2,262 2,580 2,614 
0,3 0,398 0,826 1,313 2,110 2,678 4,480 4,276 4,380 
0,4 0,464 0,836 1,585 2,558 3,431 4,823 6,774 Tn202 
0,5 0.490 1,034 1,681 2,816 3,749 5,702 10,061 oo 


Die Abb. 4 und 5 geben ferner die Veranschaulichung der Momentenverteilung, bzw. die 
Héhenlinien der Momentenflachen bei der Vierlastplatte und der Einlastplatte. Die Zuord- 
nung der mit Zahlen versehenen M-Linien zu den Werten der Momentenbeiwerte ist derart, 
daB der Linie n dem Werte c, = n- 0,3 entspricht. 


r 


Gy 


L 


Abb. 4. Héhenlinien der Momentenfliche einer aufliegenden Abb. 5. Héhenlinien d a i 
I BAY, er Momentenfliche d 
Vierlastplatte (A = 1,5, & = 1/3, m = 1/2). (A= 1,5, & = 1/2, ‘m= 1/9). area 


, 
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II. Die querkraftfrei eingespannte Vierpilzplatte. 
6. Die Lewesche Ausgangsreihe, Bei querkraftfriex Kinspannung, die z. B. den Bedingungen 


dw 0 [0?w 
39> delge + 2 — Gl = ee) 


an den Randern xp = 0 und x = aentspricht, kénnen in den Reihenentwicklungen nur die Funk- 
tionen cosnz¢ und cos nw n/A auftreten. Wenn wir cine gleichmaBige Belastung p, je Flachen- 
einheit und vier Pilzstiitzen mit rechteckigen Stiitzflachen bei der Stiitzkraft q voraussetzen, 
so kénnen wir die Belastungsfunktion p(x y) den Ansatz machen 


P(* ¥) = Po — 4f(%) aly) 
mit (37) 
f(*) == Ag+ >” Am cos mzc~, ay) => Bot > Bi cosnn. 


Dann erhalten wir nach der Fourierschen Methode folgende Werte der Koeffizienten: 


a b 
2 Pa ifr B 
A=Z[fe)dz=8a, B=F | sy)dy =87, 
0 0 
2 ‘i x 2 5 x|%+a0 x |a—%+aK 
An == | f(%) cos m To dx ge a [sin es ae asm marl 
0 


' 4 4 
=— (1+ cos m Jt) sin m7 a cosma &y = —— [1 + (— 1)™]sin maacosmaz &, 
m7 


is Dane es Mo 
B, eka 1)"] sinna 7 COs N a. 
Man sieht daher, daB in diesem Falle die Glieder mit ungeraden m bzw. n verschwinden, und 
daB fiir gerade m oder n die Koeffizienten die Werte 


ce Sis B No a, 5 
As, == sin mm a cos mmf, Bn a erg, CMe tee (m,n = 2,4,6...) (38) 


erhalten. Damit wird aber 


ff) =4a+—. 2 = sin mao. cos mx &, cos ma & ; 


m=2,4, 6... 
Baise B 5" ee ce” Lan are” 
gn) ca Fh ee : for n A A ho 
und die Lastfunktion erhalt die Form 
1 ee “ 
p(én) =— 2 Po(za > ai mizacosmma &cosmz é 
™m 
7 1 or B No wis 39 
+ =, D) = sin na 7 cos na) cos na (39) 


; B n 
-- wD sin m7 % cos m7 Ecos ma § sin n1—cos naa cos na] : 
wenn wir als Gleichgewichtshedingung 
l6gapa*=—p,ab oder 16qga0B =pyd 


beriicksichtigen. Wenn die Stiitzflachen unendlich klein werden, so geht die Belastungs- 
5* 
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funktion itiber in 


n 


p (&m) =—2p(>) cos mz &) cos ma & + >" cos n 1" cos nm 


m 


+ 2 >">’ cos mz & cos mm cos nm! cos na] (nn = 254, 07) (40) 


Macht man fiir w einen entsprechenden Ansatz mit unbestimmten Koeffizienten, so erhalt 
man durch Einsetzen in die Differentialgleichung und Vergleich beider Seiten fiir w eine Summe 
zweier einfacher und einer Doppelreihe. Nach Hinzufiigung einer Konstanten hat man also 


EDEN Meces ne coe 
W = Wy Nat cos mx £,cos ma E+ >’ cos nv? cos 7 
m™m 


n 


1] UI 
cos ma &,cos mx & cos nt —" cos n 1 (41) 


A A =9.4.6...)5 
se OS OS (n2 -- m? 42)? (m,n ‘ ) 


= Wy + Wy + Wy + Ws - 


7. Umformung der Reihen. Die Doppelreihe kénnen wir zunachst auf eine einfache Reihe 
reduzieren, indem wir etwa die auf den Index n und die Variable 7 sich beziehende Reihe 
summieren mit Hilfe der unter (8) angegebenen Formel. Setzen wir namlich unter der Voraus- 
setzung 7 < % 


No nel No Non 
cos NT Cos MI —- — z (cos ma 7 + cosnz 7 3 


so erhalten wir folgende Gleichung: 


J 1 
ee . az m HA Coj mx rp Co man 
Pa mend g., 8 (m A)® Sint m ne Sin m na 
+ Gof max(F— Ig) (Co mary — mae Sin man) Gey 
i 1 
+ mam Coj man Sint ma (-5 m) (mai 


Das letzte von 7 unabhangige Glied auf der rechten Seite fiihrt nun, wie man erkennt, wieder 
zu der Reihe — w, , die sich daher mit der ersten Reihe unserer Summe aufhebt. Damit bleibt 
aber auBer der Konstanten nur die Summe einer einfachen, nur » enthaltenden und einer beide 
Variablen enthaltenden einfachen Reihe iibrig. Wenden wir aber auf die erste Reihe w, die 
Formel (15) an, so ergibt sich 


cos nv cos nw —- 


A A 


ibe) -(9) +28) +202) 6-0) <7] 


und daher_insgesamt fiir die Durchbiegung, wenn wir noch statt m die Zahl 2n einfiihren, 


n=2,4,6... 
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wobei dann n die ganze Zahlenreihe 1, 2, 3 usw. durchl auft, 


4 
wy = yt FPG [nt + A? 9? — 6 9 (A= M0) 9? + 8A — M0)? 


16 N (nz)? Ginnaa — Ginn wA 


__ Po ab Ss Serene | ea 
piel PRs Us 


+ Coj nx (A— 2 mp) (Cof 2nxy —2nay Gin 2 nz) 
+2na% Gin Bn (5 — no) Gof 2m 29) 


Aus der im Gebiet I giiltigen Funktion wy; ergeben sich die Ausdriicke fiir wy im Gebiet II, 
bzw. der Ausdruck fiir wy im Gebiet III nach der oben angegebenen Regel durch Vertau- 
schung von 7 und 7, bzw. durch Ersetzung von 7 durch A — n. 


Abb. 6 gibt eine Veranschaulichung der Héhenlinien der nach (44) berechneten Biegeflache 
fiir den Fall A = 1,5, & = 1/3, yn) = 1/2. Die Stiitzpunkte treten als singulire Stellen fiir die 
Anschauung nicht hervor; erst in den zweiten Ableitungen von w macht sich die Singularitat 
bemerkbar. 


- Es verdient hervorgehoben zu werden, daB w, ein partikulares Integral unserer Gleichung 
darstellt, namlich die Durchbiegung eines Plattenstreifens von der Breite b, der in der x-Rich- 
tung unendlich lang ist. Dieselbe Form erhalt man, wenn man nach der elementaren Balken- 
theorie die elastische Linie eines gleichmaBbig belasteten Balkens von der Lange 6b berechnet, 
der durch zwei, um y, von beiden Enden abstehenden Auflagern gestiitzt und an den Enden 
(durch die angrenzenden Teile) so eingespannt wird, da die Tangenten hier horizontal bleiben 
(vgl. Abb. 7). Die beiderseitigen Kinspannungsmomente erhalt man auch aus der allgemeinen 
Momentensumme, fiir die sich aus (44) folgender Ausdruck herleitet: 


Abb. 6. Héhenlinien der Biegeflache der querkraftfrei eingespannten Vierpilzplatte (4 = 1,5, & = 1/3, = 1/2). 
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Abb. 7. Plattenstreifen von der Breite b, gleichmaBig belastet, symmetrisch gestiitzt und an den Enden eingespannt 
(Partikularlésung w,). 


My => pat : = [a(2 4% — 3 B29 = 208) — 


eet patton Be 


n Ginn ah 


Wales, 0 as 


(45) 


entsprechend fiir die Gebiete II und III nach der oben ausgesprochenen Regel. Setzt man im 
ersten-Teil auf der rechten Seite 7 = 0,s0 kommt als Einspannungsmoment der Wert 


if eee, 1 
My = Po (20 — mb — 5 7) = [0 (6 — 90) — FH (46) 


Eine gewisse Bestatigung der Rechnung gibt die Herleitung des Falles & = 1/A = 1/2, bei 
dem die vier Stiitzen in der Mitte des Rechtecks zusammenfallen. Wir erhalten dann 


b\2 3|2 pab (— 1)» cos 2na& 
Tay aes (sl "ie eciens — (xn) Ginna 


x [Cof2nay (1+ naAGoj tg nwd)—2nay Gin2nayJ, + (47) 


oA hoe £ SK (—1)"cos2na& Coj 2na7 
M=> pe? ae gale B= 248) = >, n Ginnaa 


n 


Das ist der an anderer Stelle behandelte Fall einer durchlaufenden Pilzdecke mit einem regel- 
maBigen Gitter von Stiitzkraften. Die Formeln gehen in die frither gegebenen Formeln iiber, 
wenn man a und b durch 2 aund 2 bersetzt. 1. 


8. Darstellung des Momentes durch 3:-Funktionen. Um den Ausdruck fiir das Moment 
umzuformen, gehen wir aus von der identischen Gleichung 


cos2nz&cos2nz&Coj2nz7 Coj 2 na(S —') =F Be eos 2 na(E—% +45] 
+ cos 2 na(¢ +f, — ct cos 2 na(¢—~%— i) + cos 2 na(é—fy+ iz)| : 


Benutzen wir ferner die Reihenbeziehung fiir den Logarithmus der $)-Funktion 


See = In JT (1—q?") — In 8(2), se 


—7 Ginn awa 
so erhalten wir fiir das Moment die Darstellung 
M =+ ae Any — 29? — 272 =a | 
+] Helm (0 — 65 + iz) Bo(E +f — ¥-7) O66 —& — iz) a(o ++i 2) 
Nun besteht aber die weitere Relation 


Reln d(vF Tid) =Ingq74T wey + Reln lr) =F AIF amt Tela Hr), (60) 


1 Vgl. W, Miiller, Ing.sArch. 20 (1952), S. 2874. - 
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wenn WIr jetzt v — v, + 1 v,setzen. Damit geht der Ausdruck fiir das Moment iiber in 


1 afl I 
a tla (1 — a) — 2 + 08) 
il = i 51 

re R_ In Co (4 == Cy) OG ++ oo) H(C —%) B(o + ca | | a 


wo jetzt die GréBen 7 und y symmetrisch auftreten. Mit den Werten A = 1,5 wird g = 0 00898 
(7/4) n§ = 09,5236. Daher haben wir, wenn wir alle Konstanten zusammenfassen : { 


1 pa®fl - ae 
Mm 5 Oe FG Meln AC —2) C+.) CB) HC +E) 


— 5 80,6159] = 7 2“, | 


(52) 


Abb. 8. Héhenlinien der Momentenflache der querkraftfrei eingespannten Vierpilzplatte (A = 1,5, & = 1/3, m) = 1/2). 


Wir haben Momentenzahlen fiir einige Werte von £ undy in dem vierten Teil der Platte berechnet 
und in der Tabelle.3 zusammengestellt. Ferner haben wir in Abb.8 die Héhenlinien der Momen- 
tenfliche gezeichnet. Dabei ist zu beachten, da® lings der Linie Null ein Ubergang von posi- 


Tabelle 3. ¢m-Werte der querkrafifrei eingespannten Vierpilzplaite. 


x 0 0,1 0,2 0,3 0,4 0,5 | 0,6 0,7 | 0,75 

Bl ac a J ee | 
0 0,284 0,267 0,233 Onl? 0,077 0,014 —0),009 — 9,055 —(0,069 
O51 0,279 0,261 0,225 0,135 0,049 —(0,021 —),040 —0,077 —(0),087 
0,2 0,268 0,247 0,203 0,088 —(0,041 | —0,161 —(,139 | —0,145 —0,141 
0,3 0,255 0,231 0,174 0,035 —0,165 == (52000 | —Onalia —(,216 
OFS 0,251 OR 27 0,171 0,023 —(,196 | — w —,299 | —0,225 —(0,227 
0,4 0,244 0,219 0,162 0,012 —0,185 | —0,393 —0,287 | —0,234 —0,321 
0,5 0,241 0,215 0,157 0,013 —0,152 | —0,273 | —0,255 —0,235 —0,232 
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tiven zu negativen Momenten stattfindet, entsprechend dem Wechsel zwischen konvexer 
und konkaver Kriimmung, der Flache. Die Zuordnung der M-Linien zu den Werten von c,, 
ist derart, daB die Linie mit der Nummer n der Wert n, 0,07095 zugehért. In den Eckpunkten 
der Platte (§ = 7 = 0) haben wir die gréBten Werte c,, = 0,2838, in der Mitte der Platte 
(€ = 0,5; 7 = 0,75) wird c, = — 0,232. Die vier Stittzpunkte sind singulare (Wirbel-)punkte 
im Felde der M-Linien, in denen das Moment negativ unendlich wird. 


9. Schlu8bemerkung. Die in der vorstehenden Arbeit angewendeten Grundsatze geben 
den Zugang zu mehreren Erweiterungen, auf die wir nur kurz hinweisen wollen. Zunachst kann 
der Fall der endlichen, rechteckigen Stiitzflachen ohne groBe Schwierigkeiten von dem anfangs 
gegebenen Ansatz aus weiter entwickelt werden, wobei es dann erforderlich ist, bei Behandlung 
der 7-Reihen auf die Verschiedenheit z. B. der w-Funktionen in fiinf verschiedenen Streifenge- 
bieten Riicksicht zunehmen. Ferner liegt es nahe, die Untersuchung auch auf andere Falle der 
Mehrpilzplatte, etwa auf die Neunpilzplatte auszudehnen, wozu alle Voraussetzungen gegeben 
sind. Der Verf. behalt sich vor, auf einige dieser Probleme spater zuriickzukommen. 
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